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principally  by  L.  M.  Linson,  B.  Kivel,  S.  Powers,  and  L,  F.  Cianciolo. 
Volume  II  contains  Chapters  1  and  2  prepared  principally  by  G.  M.  Wey]  and 
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ABSTRACT 


Several  aspects  of  the  analysis  of  barium  ion  clouds  are  presented 
including  ion  cloud  modeling,  comparison  of  radar  and  optical  data,  and 
correlation  of  data  with  theory.  A  quantitative  model  has  been  developed 
from  which  various  properties  of  barium  ion  clouds,  primarily  their  size, 
time  history  of  the  peak  electron  concentration,  and  height-integrated 
conductivity  can  be  estimated.  These  estimates  are  in  excellent  agreemtnt 
with  the  observations  of  Secede  ion  clouds.  The  modeling  of  the  radiation 
transport  in  the  ion  cloud  has  been  confirmed  by  a  Monte  Carlo  calculation. 
The  optical  analysis  of  Spruce  at  R  +  14  min  has  been  corrected  for  effects 
of  cloud  geometry  and  the  results  for  ion  density  are  found  self-consistent 
and  in  agreement  with  radar  measurements.  A  summary  review  is  given  of 
the  values  and  limitations  of  photographic  data.  The  current  status  of 
theoretical  understanding  of  the  dynamics  of  barium  ion  clouds  is  reviewed. 
Particular  attention  is  given  to  their  motion,  deformation,  and  the  properties 
of  striations  including  onset  time,  scale  size,  and  dissipation. 
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technical  summary 

VOLUME  I 

Chapter  I  describes  a  quantitative  model  from  which  the  properties 
of  ion  clouds  can  be  estimated.  Of  particular  importance  are  the  ion  cloud 
size  perpendicular  to  the  magnetic  field,  the  peak  electron  content  on 
magnetic  field  lines,  the  time  history  of  the  peak  electron  concentration, 
and  the  height-integrated  Pedersen  conductivity  of  the  ion  cloud.  These 
properties  are  given  as  u  function  of  the  altitude  of  release  (140  km  to 
300  km)  and  chemical  payload  weight  in  the  range  0.  1  kg  to  400  kg  and 
are  found  to  be  in  excellent  agreement  with  observations  of  Secede  ion 
clouds.  The  effects  of  the  neutral  wind,  the  ambient  electric  field,  the 
ionospheric  conductivity,  the  latitude  of  release,  and  the  time  of  day  and 
year  are  also  discussed.  The  results  of  this  modeling  are  useful  for  future 
test  planning. 

Chapter  2  reports  on  a  radiation  transport  study  that  confirms 
the  optical  modeling  used  in  Chapter  3,  Radiation  diffusion  in  a  slab  of 
purely  scattering  atoms  with  only  one  transition  is  calculated  using  a 
Monte  Carlo  method.  Partial  redistribution  because  of  the  Doppler  effect 
is  included.  The  scattering  is  assumed  to  be  isotropic  and  incoherent. 
P.esults  are  obtained  for  1)  the  amount  of  diffuse  reflection  and  transmission 
as  a  function  of  optical  depth,  and  2)  the  angular  distribution  and  spectrum 
of  the  reflected  component.  The  results  are  found  to 
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obtained  from  a  differential  equation  for  radiation  diffusion  without  redis¬ 
tribution.  Tl. j  Jight  reflected  is  characterized  by  an  aijedo  close  to  one 
at  the  line  center  when  the  optical  depth  is  large  compared  to  one.  The 
reflected  intensity  increases  with  the  amount  oi  scattering  material 
corresponding  to  an  increase  in  an  effective  line  width.  The  total 
amount  of  diffusely  transmitted  radiation  is  insensitive  to  the  amount 
of  scattering  material  and  has  two  spectral  bands  centered  in  the  line 
wings  where  the  optical  depth  is  about  two.  The  angular  distribution  of 
reflected  radiation  Is  Lambertian. 

Chapter  3  contains  a  correlation  of  Spruce  data.  Optical  analysis 

for  event  Spruce  at  R  +  14  min  has  been  improved  to  take  account  of  the 

angles  between  the  cloud  face  and  the  solar  illumination  direction  and  the 

lines  of  sight  from  Tyndall  and  Barin.  The  correction  for  the  dependence 

of  these  angles  leads  to  consistent  ion  densities  and  better  agreement  with 

the  radar  determined  density.  The  peak  ion  column  density  perpendicular 

12  -2 

to  the  cloud  face  is  determined  to  be  ~3  x  10  cm  .  Comparison  of  the 
extent  of  the  ion  cloud  as  observed  by  the  Thomson  scatter  radar  measure¬ 
ment  with  the  longer  extent  observed  optically  leads  to  an  estimate  of  cloud 
width  of  about  2  km  at  a  distance  of  about  5  km  west  of  the  trailing  edge. 

Chapter  4  reviews  the  use  of  photographic  data  for  estimating 
ion  density  and  ion  inventory.  Methods  for  analysis  of  optically  thin  and 
thick  clouds  are  considered.  An  effective  line  width  method  is  preferred 
for  thick  clouds  in  order  not  to  underestimate  ion  column  density  by  using 
the  linear  relation  valid  for  optically  thin  clouds.  In  either  case  a  semi- 
empirical  theory  is  required  because  of  the  many  unknowns.  The  unknowns 
are  grouped  as  problems  of  1)  data  reduction  and  2)  interpretation  of 
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reduced  data.  Data  reduction  problems,  which  include  film  reciprocity, 
dependence  of  sensitivity  on  film  temperature,  and  atmospheric  trans¬ 
mission  may  be  alleviated  using  known  st£  r  intensities  as  standards. 

Unknowns  in  interpretation  include  the  solar  spectral  flux  and  radiation 
transport  details.  The  good  results  of  the  data  correlation  in  Spruce  at 
R  i  14  min  lead  to  the  conclusion  that  despite  unknowns  a  scn'i-empirical 
analysis  of  photographic  data  can  be  used  to  measure  barium  cloud  ion 
column  density. 

Chapter  5  consists  of  a  review  of  our  current  understanding  of 
barium  ion  cloud  dynamics.  It  is  found  that  the  most  important  parameter 
that  determines  the  behavior  of  ion  clouds  is  the  ratio  of  the  height-integrated 
Pedersen  conductivity  of  the  ion  cloud  to  that  of  the  ambient  ionosphere.  A 
theoretical  treatment  of  the  motion  of  ion  clouds  shows  that  low  conductivity 
clouds  tend  to  move  with  the  local  E  x  B/B  velocity  while  high  conductivity 
clouds  tend  to  move  with  the  local  neutral  wind.  Quantitative  agreement 
between  theory  and  observations  is  achieved  by  assuming  that  the  conjugate 
ionosphere  in  the  southern  hemisphere  is  contributing  to  the  ambient 
conductivity.  A  simple  model  is  shown  to  explain  several  aspects  of  the 
deformations  of  ion  clouds  including  their  steepening  on  the  backside, 
elongation  in  the  direction  of  motion,  and  narrowing  transverse  to  the  direction 
of  motion.  The  properties  of  striations  including  their  onset  time,  scale 
size,  and  dissipation  are  treated.  It  is  found  that  the  normalized  onset  time 
is  a  minimum  for  clouds  of  intermediate  conductivity.  The  simple  two- 
dimensional  model  predicts  the  scale  size  of  striations  when  the  effects  of 
the  cloud  polarization  field,  conductivity,  and  finite  size  are  included  as  well 
as  electron  collisions.  The  dissipation  of  striations  is  controlled  by  a  combi¬ 
nation  of  diffusion  and  convection.  Areas  where  more  theoretical  work  is  needed 


arc  also  d  is  cussed. 
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CHAPTER  1 

QUANTITATIVE  ION  CLOUD  MODELING 


J .  QUANTITATIVE  ION  CLOUD  MODELING 


Lewis  M.  Linson 


A.  INTRODUCTION 

We  will  present  quantitatively  the  variation  of  ion-cloud  properties 
such  as  si/.c  and  density,  as  a  function  of  the  altitude  of  release  (from  140 
to  300  km),  and  payload  in  the  range  from  0.  1  kg  to  400  kg.  Knowledge  of 
these  cloud  properties  over  this  range  of  parameters  will  be  useful  for 
future  test  planning  and  is  necessary  as  an  input  for  theoretical  work  that 
describes  the  dynamics  and  distortions  of  barium  ion  clouds.  The  plasma 
aspects  of  barium  ion  cloud  behavior  are  described  in  Chapter  5.  The 
model  described  below  has  been  chosen  to  give  optimum  agreement  with 
the  best  available  optical  and  radar  data.  A  correlation  of  the  results  of 
this  modeling  with  the  data  obtained  during  the  Secede  I  and  II  tost  series 
conducted  at  mid- latitudes  will  be  discussed. 

In  the  following  three  sections  we  will  describe  the  ion  cloud  pro¬ 
perties  at  early,  intermediate,  and  late-times  after  release.  Early-time 

behavior  refers  to  less  than  a  minute  after  release.  Intermediate-time 
behavior,  which  is  dominated  by  ambipolar  diffusion  parallel  to  the  magnetic 
field,  typically  will  last  from  one  to  15  minutes  after  release.  Long-time 
behavior  typically  occurs  after  some  tens  of  minutes.  In  Section  E  we 
will  also  discuss  the  time  variation  of  the  electron  density  as  obtained 
from  the  above  modeling  for  several  typical  ion  clouds.  In  Section  E  we 
compare  the  results  of  the  modeling  with  optical  and  radar  data.  The 
height- integrated  Pedersen  conductivity  of  the  ion  cloud  will  also  be  pre¬ 
sented  as  a  function  of  altitude  and  payload.  In  Section  F  we  will  describe 


the  influence  of  several  ambient  parameters  whose  variation  has  not  been 
included  in  the  modeling.  In  particular,  the  effects  of  the  neutral  wind, 
the  ambient  electric  field,  the  ionospheric  conductivity,  the  latitude  of  the 
release,  and  the  time  of  day  and  year  will  be  discussed.  We  will  also 
indicate  the  importance  of  the  results  of  this  modeling  for  determining  the 
behavior  of  ion  clouds  as  they  are  influenced  by  the  ambient  electric  field. 
The  limitations  of  the  modeling  will  also  be  discussed. 

B.  EARLY- TIME  BEHAVIOR 


By  the  early-time  behavior  of  the  ion  cloud  we  mean  the  time  during 
which  the  ion  cloud  is  formed  by  the  photoionization  of  barium  atoms.  The 
deposition  of  barium  ions  from  the  barium  vapor  typically  takes  place  in 
less  than  a  minute.  The  single  most  important  parameter  that  describes 
the  early-time  behavior  is  the  collision  time,  r  ,  of  barium  with  the  ambient 
atmospheric  constituents.  The  mobility,  /it,  of  barium  ions  in  nitrogen  gas 
has  been  measured,  *  from  which  we  may  obtain  the  barium  ion-nitrogen 
collision  time, 


mBa  8.  6  x  1(T  , 

T  =  -  s  - — - seconds 

6  N 

n 

where  mga  *s  the  atomic  mass  of  barium  and  e  is  the  magnitude  of  the 
electron  charge.  The  collision  cross  section  has  been  assumed  to  vary 
inversely  as  the  square  root  of  the  temperature,  and  is  a  weighted 
atomospheric  concentration  in  cm  ^ , 


(1) 


N  =  N_t  +  N0  +  0. 8N^. 
n  N2  °2  ° 


(2) 


It  has  been  assumed  that  the  collision  cross  section  for  oxygen  molecules 
and  atoms  is  the  same  as  it  is  for  nitrogen  molecules.  The  fact  that  the 
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lighter  oxygen  atoms  are  less  efficient,  in  stopping  a  heavier  barium  ion 

]  h 

is  expressed  by  the  appropriate  Langevin  factor,  (1  +  mt,  m  )  '  x 

u  cl  N  2 

(1  +  na j-j  ^  h1q)  ^  ^  ~0.  8,  where  m and  mQ  are  the  masses  of  the 

nitrogen  molecule  and  oxygen  atom  respectively.  We  also  assume  that 

barium  atoms  have  the  same  collision  cross  section  with  air  and  hence  T 

given  by  Eq.  (1)  will  represent  the  barium  atom-neutral  air  collision  time. 

The  CIRA  1965  atmospheric  model  6  at  1800  hours  has  been  used 

to  obtain  the  neutral  temperature,  mass  density,  and  N2>  0?,  and  O 

densities  as  a  function  of  altitude.  The  model  6  is  appropriate  for  the 

solar  flux  conditions  that  existed  during  the  second  half  of  January  1971. 

We  have  assumed  that  the  ion  temperature,  is  the  same  as  the  neutral 

temperature,  Tn,  since  the  ions  and  neutrals  are  in  close  thermal  contact. 

Our  model  for  the  altitude  dependence  of  the  electron-to- ion  temperature 

2 

ratio  was  influenced  by  the  recent  paper  by  Wand  .  We  have  taken  the 
electron-to- ion  temperature  ratio,  Te/T^,  as  varying  linearly  from  one 
at  1  30  km  tc  a  value  of  1 . 6  at  21  0  km  and  then  decreasing  linearly  to  1 .  2 
at  300  km.  The  barium- neutral  and  barium- ion  ambipolar  diffusion 
coefficients,  Dn  and  ,  respectively,  are  given  by 

k  T  /  T  \ 


D 


m 


B 


~  T;  DH  =  Dn  1  +  tH 
a  \  n  / 


(3) 


and  are  shown  in  Fig.  1  along  with  the  ion-neutral  collision  time,  T,  as  a 
function  of  altitude,  k  is  Boltzmann's  constant. 

We  also  need  to  know  the  rate  at  which  barium  atoms  are  photo- 


ionized  and  oxidized  by  the  ambient  air.  Haerendel  of  the  Max- Planck 
Institute  has  consistently  reported  photoionization  times  of  order  19  seconds 
while  Best  of  AFCRL  reports  a  time  closer  to  30  seconds.  In  our  model 


ALTITUDE  (km) 


wo  will  choose  24  seconds  as  being  within  20%  of  the  correct  value.  We 
shall  also  Lake  the  oxidation  rate  constant  of  barium  by  atmospheric 
to  be  7.  5  x  10"  ^  cm^/s  cc.  The  fraction,  f,  of  available  barium  atoms 
t Hat  are  ionized  is  equal  to  the  ratio  of  the  photoionization  rate  to  the 
total  loss  rate  of  barium  atoms  due  to  both  photoionization  and  oxidation, 

(  =  - i - r? - ■ 

1  +  1.8x10  V  Nn 

U2 


The  time  constant,  t.,  is 
’  i 


24 


li  24  f  1  +  1.8  x  10" 9  N 


seconds 


O. 


(5) 


and  is  shown  as  a  function  of  altitude  in  Fig.  2. 

We  first  give  a  qualitative  description  of  how  the  ion  cloud  is  formed 
before  proceeding  with  a  quantitative  model  that  seems  to  fit  the  data  very 
well.  We  shall  be  concerned  here  with  the  chemical  mixture  that  was  used 
during  the  Secede  II  test  series.  The  barium  thermite  was  based  on  the 
reaction  of  2.  5  moles  of  barium  per  mole  of  cupric  oxide  according  to  the 
formula. 


2.  5  Ba  +  CuO  — •  Ba  O  +  Cu  +  1.5  Ba.  (6) 

In  addition,  1.  8%  of  the  thermite  weight  was  barium  azide.  As  a  result, 
roughly  half  of  the  chemical  payload  weight  was  available  as  excess  barium, 
amounting  to  ~10  D  barium  atoms  per  48-kg  barium  release.  Ground 
tests  have  shown  that  most  of  the  gaseous-iiquid  two-phase  flow  is  vented 
from  the  canister  in  a  few  tenths  of  a  second.  Optical  data  at  high  altitudes 
indicate  that  the  exit  velocity  is  typically  of  the  order  of  1  km/sec.  Hence, 
a  mixture  of  liquid  Ba  O,  Cu,  and  Ba  and  their  vapors  rapidly  expands 
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until  the  gaseous  products  reach  equilibrium  with  the  ambient  atmospheric 
density.  Since  barium  is  the  most  volatile  product,  most  of  the  vapor 
should  consist  of  barium  atoms.  Studies  of  the  barium  ejection  process 
ha-'re  led  to  the  conclusion  that  most  of  the  barium  atoms  that  are  produced 
are  formed  by  evaporation  from  hot  liquid  barium  droplets  ejected  from 
the  canister.  The  cooled  liquid  droplets  that  do  not  vaporize  form  the 
particulate  ring  whose  center  of  mass  follows  the  rocket  trajectory. 

After  the  rapid  expansion  by  this  barium  vapor  has  ceased,  the  barium 
atoms  begin  diffusing  through  the  ambient  atmosphere.  At  the  same  time, 
solar  radiation  begins  exciting  ground- state  barium  atoms  into  metastable 
states.  The  barium  atoms  in  the  excited  states  are  then  lost  by  the  com¬ 
peting  processes  of  ionization  and  oxidation  by  molecular  oxygen,  with  the 
time  constant  E  given  by  Eq.  (5).  Once  the  ions  have  been  created  they  are 
tied  to  the  earth1  s  magnetic  field  lines  and  can  no  longer  diffuse  perpen¬ 
dicular  to  the  magnetic  field.  Hence  the  scale  size  of  the  ion  cloud 
transverse  to  the  magnetic  field  should  bear  some  relation  to  the  trans¬ 
verse  scale  size  of  the  neutral  cloud  at  early  times.  Also  within  a  few 
time  constants  E,  the  barium  atoms  either  become  barium  ions  or  have 
been  oxidized  to  become  BaO. 

We  shall  be  concerned  here  only  with  the  dense  ion  cloud  formed 
from  the  initial  ionization  of  barium  atoms.  We  shall  not  be  concerned 
with  the  less  dense  ionization  that  is  formed  by  the  photodissociation  of 
BaO  and  subsequent  ionization  of  the  barium  atom.  This  second  process 
lakes  place  on  a  time  scale  at  least  two  orders  of  magnitude  slower  than 
the  primary  process  and  gives  rise  to  densities  two  orders  of  magnitude 
less  than  the  peak  densities  of  concern  to  us  here. 
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Now  iet  us  discuss  a  quantitative  model  for  the  formation  of  the 
main  ion  cloud,  the  formation  of  which  depends  on  the  early-time  behavior 
of  the  neutral  vapor.  Initially,  the  barium  vapor  at  high  density  compared 
to  the  ambient  can  be  taken  to  fill  a  sphere  whose  outer  edge  expands 
radially  at  1  km/s,  a  much  larger  velocity  than  a  barium  atom  thermal 
velocity  at  ambient  temperature.  At  a  later  time,  after  the  barium  atoms 
have  lost  most  of  their  directed,  radial  velocity  and  have  expanded  so  that 
their  concentration  is  small  compared  to  the  ambient  concentration,  the 
distribution  of  barium  atoms  can  be  taken  to  be  a  spherical  gaussian 
expanding  by  diffusion  through  the  ambient  air  at  less  than  a  thermal 
velocity.  There  exists  no  model  that  incorporates  the  physical  processes 
responsible  for  the  transition  between  this  initial  and  later-time  behavior. 
In  order  to  develop  a  quantitative  description  of  this  transition,  we  shall 
adopt  the  following  model.  We  shall  assume  that  the  initial  mass  of  vapor 
uniformly  fills  a  sphere  of  radius  R(t)  and  that  the  initial  momentum  of  the 
vapor  mass  is  shared  with  the  mass  of  ambient  air  contained  inside  this 
radius.  The  assumption  that  the  velocity  within  the  sphere  increases 
linearly  with  radius,  r,  from  zero  at  the  center  to  dR/dt  at  r  =  R  satisfies 
flux  conservation.  We  define  an  equal  mass  radius,  aQ,  as 

a  =  (3M/47rp)1/3  (7) 

where  M  is  the  mass  of  the  vapor  and  p  is  the  mass  density  of  the  ambient 
atmosphere.  The  equation  expressing  the  conservation  of  momentum 
described  above  can  be  written 


1 


(8) 
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giving  the  time  behavior  of  the  radius  of  the  vapor  bubble.  When  the  dR/dl 


has  decreased  to  a  typical  barium  thermal  velocity  of  0.  3  km/s,  the  vapor 

bubble  has  a  radius  R  =  1.  27  a  and  it  takes  1.  92  a  seconds  to  reach  this 

o  o 

radius.  Until  R  reaches  a  at  t  =  1.  25  a  seconds,  the  barium  vapor  is 

o  o  1 

denser  than  the  surrounding  atmosphere  and  will  certainly  be  optically  thick 

to  the  metastable  exciting  solar  radiation.  Hence,  we  shall  assume  that  the 

formation  of  barium  ions  begins  only  1.  25  a^  seconds  after  release. 

Equation  (8)  underestimates  the  radius  of  the  neutral  cloud  as  a 

function  of  time.  A  single  barium  atom  traveling  at  1  km/s  will  have  a 

mean  free  path  of  r  km  where  r  is  the  collision  time  given  in  Eq.  (1)  and 

plotted  in  Fig.  1.  A  high-  velocity  barium  atom  will  lose  only  l/3  of  its 

velocity  in  a  single  collision  with  a  typical  air  molecule  which  is  five  times 

lighter.  Hence  we  might  expect  the  mean  veolcity,  V,  of  a  cluster  of 

barium  atoms  to  behave  exponentially  with  time  as  exp  (-t/3r  )  during  the 

0.  67  aQ  seconds  that  the  barium  vapor  is  less  dense  than  the  surrounding  air. 

Thus  we  assume  that  the  radius  of  the  neutral  barium  vapor  is  larger  than 

1.  27  a  ,y 
o  1 


( 

/ 


,0.  67  a 


V(t)  dt  =  ir  [  1  -  exp(-0.  67  a  /3r)]  km 


(9) 


After  1.  92  seconds  after  release  we  assume  that  the  neutral  vapor  con¬ 
tinues  to  increase  in  size  by  ordinary  diffusion  giving  rise  to  a  typical 
gaussian  radius  (defined  as  the  e-folding  radius  for  the  density)  as  a  function 
of  time  of 

A/1 


->  (t)  = 

n 


.  33  +  37  ^  1  -  exp  £-0.  67 a  J 3t  J 


+  4D  (t  -  1. 92a  ) 


) 


km.  (10) 
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While  the  neutral  cloud  is  expanding  according  to  the  model 
described  above,  ions  are  being  formed  and  deposited  on  magnetic  field 
lines.  We  assume  that  the  scale  size  of  the  ion  cloud  transverse  to  the 
magnetic  field  is  given  by  the  gaussian  radius  of  the  neutral  cloud  at  two 
ionization  time  constants  after  ionization  is  assumed  to  have  begun,  or 

al  =  an  ^ 1  * 

whe  re 

t,  =  2t.  +  1.25a  seconds.  (12) 

1  1  o 

For  our  model  we  have  assumed  that  the  mass  of  the  vapor  pro¬ 
duced  in  the  barium  release  is  l/5  of  the  total  chemical  payload  weight. 

The  equal  mass  radius  defined  by  Eq.  (7)  sets  the  initial  scale  size  for 
the  neutral  cloud  and  the  initial  time  scale  of  the  expansion  and  is  shown 
in  Fig.  3  as  a  function  of  altitude  and  payload.  The  ion  cloud  scale  size 
transverse  to  the  magnetic  field  given  by  Eqs.  (10),  (11),  and  (12)  is  shown 

in  Fig.  4  as  a  function  of  altitude  and  payload. 

23 

We  now  assume  that  6.  25  x  10  P  barium  atoms  are  produced 
where  Pis  the  chemical  payload  in  kilograms.  This  number  corresponds 
to  a  30%  vaporization  efficiency  for  the  barium  release,  which  is  consid¬ 
erably  higher  than  has  been  obtained  in  laboratory  tests  but  is  not  incon¬ 
sistent  with  some  estimates  of  the  total  ion  inventory  in  large  barium 
releases  as  we  shall  discuss  below.  The  mass  chosen  above  for  the  vapor 
corresponds  to  40%  more  than  the  mass  of  the  barium  atoms  alone.  This 


choice-  appears  to  load  to  good  agreement  of  the  model  with  observations, 
although  there  is  good  reason  to  believe  that  the  vapor  consists  almost 


i 


entirely  of  barium  atoms.  Only  the  fraction  f  given  by  Eq.  (4)  of  the  barium 

atoms  produced  will  form  ions.  An  estimate  of  a  peak  electron  density,  n  , 

23 

can  be  made  assuming  that  f  x  6.  25  x  10  x  P  ions  are  distributed  throughout 
a  sphere  of  gaussian  radius  km: 


n 


P 


f  x  6.  25  x  108  P  -3 

—37273  Cm 

al 


(13) 


The  value  of  n^  is  shown  as  a  function  of  altitude  and  payload  in  Fig.  5. 

Since  ions  cannot  diffuse  across  the  magnetic  field  lines,  the 
number  of  ions  on  each  magnetic  field  line  remains  constant.  Hence  we 
can  easily  estimate  the  peak  electron  content  deposited  on  a  magnetic  field 
line  by  integrating  the  gaussian  distribution  parallel  to  the  magnetic  field 
resulting  in 


N, 


f  x  6.25  x  10 


13 


Ta.2 


P  -2 
—  cm 


(14) 


The  maximum  number  of  electrons  cn  a  field  line  given  by  Eq.  (14)  is 
shown  as  a  function  of  altitude  and  payload  in  Fig.  6. 

C.  INTERMEDIATE  TIME  BEHAVIOR 

In  the  preceding  section  we  have  developed  a  model  that  suggests 
that  at  t  -  tj  s  after  release  given  by  Eq.  (12),  the  ion  cloud  is  a  sphere 
with  a  gaussian  radius  given  by  Eq.  (11)  and  with  a  peak  electron  density 
given  by  Eq.  (13).  Subsequently,  the  electron  density  decays  as  the  ioni¬ 
zation  spreads  out  parallel  to  the  magnetic  field  by  ambipolar  diffusion. 
Hence  th».  gaussian  length  parallel  to  the  magnetic  field  is  given  by 


-13- 


r  2  I1/2 

iL  (t)  =  +  4  D  (t  -  t^)J  km. 


The  electron  density  at  one  minute  after  release  is  then  given  by 


Nixl°  -3 

n  (60)  = - cm 

S*  iL  (60) 


(£  is  given  in  kilometers)  and  is  shown  as  a  function  of  altitude  and  payload 

- 1  /2 

in  Fig.  7.  The  electron  density  begins  to  decrease  approximately  as  t  ' 
at  times  long  compared  to 


which  is  shown  as  a  function  of  altitude  and  payload  in  Fig.  8 
D.  LATE -TIME  BEHAVIOR 

The  atmosphere  was  assumed  to  be  of  constant  density  in  the  pre¬ 
vious  section  since  the  size  of  the  ion  cloud  was  small  compared  to  a 
neutral  density  scale  height.  However,  the  large  value  of  the  ambipolar 
diffusion  coefficient  causes  the  ion  cloud  to  grow  to  an  appreciable  fraction 
of  the  neutral  atmosphere  scale  height  in  the  direction  parallel  to  the 
magnetic  field.  The  distribution  of  ions  parallel  to  the  magnetic  field  will 

become  distorted  from  a  gaussian  shape  and  assume  a  new  quasi-equilibrium 

3 

shape  as  obtained  by  Bannister  and  Davis  and  illustrated  for  the  case  of 

4 

barium  ion  clouds  by  Pendyala  , 

n(z>  =  nmax  exP  {-it  -  rr  [exp('TT-> -1]}  (18> 
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Fig.  8  The  Time  at  Which  Ambipolar  Diffusion  Parallel  to  the  Magnetic 
Field  Begins  to  Cause  a  Decrease  as  t"1/2  in  the  Peak  Electron 
Density  of  the  Ion  Cloud 
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where  x.  is  measured  from  the  peak  of  the  electron  density  and  H.  is  the 
barium-ion  scale  height,  D|j  /gy .  On  the  topside  at  large  positive  z,  the 
density  varies  as  exp  (-z/lL)  while  near  the  peak  the  electron  density  has  a 
gaussian  shape 


n(2)  '  "max^t-^A2]’  2  «  Hn 


with  a  gaussian  length  of 


iA  -  4ZH.Hn  . 


This  gaussian  length  is  shown  as  a  function  of  altitude  in  Fig.  9. 

The  profile  given  by  Eq.  (18)  is  in  quasi-equilibrium  and  indicates 
that  the  peak  ion  density  eventually  decreases  to  a  constant  value  at  a 
given  height.  This  asymptotic  value  of  the  density  in  the  ion  cloud  is  given 


N.  x  10 

"A  = - - -  <2 

tA 

and  is  shown  as  a  function  of  altitude  and  payload  in  Fig.  10.  Before  this 

equilibrium  distribution  is  reached,  the  peak  electron  density  has  been 

decreasing  inversely  as  the  gaussian  scale  length  £  (t)  given  by  Eq.  (15) 

-  1  /2 

or  roughly  as  t  .  An  estimate  of  the  time,  t^,  at  which  the  electron 

density  stops  decreasing  and  approaches  its  asymptotic  value  given  by 

1  /2 

Eq.  (21)  is  made  by  setting  £  (t^)  equal  to  (I-TH^)  '  resulting  in 


lA  "  tl  + 


H.  H  -a. 

in  1 


which  is  shown  as  a  function  of  altitude  and  payload  in  Fig.  11. 


D  3423 


20 


30 


j£A(km) 


40 


Fig.  9 


The  Gaussian  Length  Along  the  Magnetic  Field  for  a  Barium 
Cloud  in  Diffusive  Equilibrium 
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In  the  previous  discussion,  it  has  been  assumed  that  the  ion  cloud 
has  remained  at  the  altitude  of  release.  In  fact,  barium-ion  clouds  have 
descended  in  altitude  with  typical  vertical  velocities  of  20  m/s.  This 
change  in  altitude  c  -n  come  about  from  a  combination  of  factors.  The 
principal  causes  are  gravitational  drag  on  the  ions,  components  of  the 
neutral  wind,  or  V^,  in  the  downward,  or  northward,  directions, 
respectively,  and  a  westward  component  of  the  electric  field,  E  .  These 
factors  lead  to  a  downward  drift  velocity  of  the  ion  cloud  of 

?  E 

ViD  =  (VD  +  8T  )  (1  "  &  COS  9d]  +  &  (VN  Sin  0d  +  IT }  COS  dd  (23) 

where  0^  is  the  dip  angle  of  the  magnetic  field  and  Oci^cl  is  the  coupling 
K  parameter,  which  depends  on  the  conductivity  of  the  ion  cloud  as  described 

in  Chapter  5D.  The  resulting  change  in  altitude  has  important  effects  on 
the  time  dependence  of  the  peak  electron  density.  The  first  effect  is  a 
decrease  in  the  rate  at  which  f  (t)  increases  due  to  smaller  D||  at  lower 
altitudes.  Thus  the  peak  electron  density  can  be  expected  to  decrease  more 

-i/z 

slowly  than  t  if  the  downward  drift  of  the  ion  cloud  is  of  the  order  of 

H  /t.  ~  4  gr  [  obtained  by  using  Eq.  (22)]  or  larger.  This  criterion  becomes 
n  A 

easier  to  satisfy  at  lower  altitudes.  The  second  effect  is  an  increase  in  the 
asymptotic  value  of  the  peak  electron  density,  n^,  since  ft ^  also  decreases 
at  lower  altitudes. 

We  are  now  in  a  position  to  describe  the  time  variation  of  the  peak 
electron  density  in  an  ion  cloud.  We  assume  that  the  ions  are  continuously 
created  with  a  time  constant  t.  given  by  Eq.  (5)  and  deposited  in  a  sphere  of 
gaussian  radius  a^  given  by  Eq.  (11)  for  times  less  than  t^  given  by  Eq.  (12). 
At  times  greater  than  t^,  the  electron  density  begins  to  decrease  inversely 
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as  the  gaussian  length  £  (t)  given  by  Eq.  (15)  increases.  The  electron 
density  decreases  until  times  t^  given  by  Eq.  (21)  at  which  time  the  density 
begins  to  approach  its  asymptotic  value  n^  given  by  Eq.  (21).  Thus  we 
have 


n(t)  =  ^ 


[]  ‘  exp(-(t-1.25ao)/ti)]  ; 


Ni 

JTF(t] 


[l  -  exp  (-{t  -  1. 25  aQ)/ti)]  ; 


n . 


t  <  1.25a 

o 

1.25  a  <  t  <  t, 
o  I 

tj<  t<  t^ 

,>2tA 


(24) 


The  time  histories  of  the  peak  electron  concentration  obtained  in 
this  manner  for  several  typical  cases  of  interest  are  shown  in  Figs.  12  and 
13.  Figure  12  shows  the  time  history  for  48  kg  payloads  at  altitudes  of  150, 
190  and  250  km,  while  Fig.  13  shows  a  comparison  of  several  different 
payloads  at  around  190  ;  m  altitude.  Note  that  the  peak  electron  density 
given  by  expression  (24)  is  typically  approximately  10%  less  than  the  peak 
given  by  Eq.  (13)  due  to  the  fact  that  all  of  the  barium  atoms  have  actually 
not  been  photoionized  by  t  =  2t^.  It  should  be  pointed  out  that  it  is  possible 
for  the  peak  electron  density  to  increase  siowly  at  late  times  due  to  a  slow 
change  in  altitude.  This  possibility  is  indicated  by  the  shaded  widths  of 
the  lines  at  late-time  corresponding  to  a  50-km  descent  in  altitude  for  the 
high  altitude  releases.  In  general  it  has  been  observed  that  the  location  of 
the  peak  in  the  electron  density  does  tend  to  decrease  in  altitude  with  time. 
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Since  the  asymptotic  value  of  the  density,  Eq.  (21),  varies  inversely  as 
the  asymptotic  gaussian  length  £  ^  given  by  Eq.  (20),  the  late-time  electron 
density  can  slowly  increase  as  the  cloud  descends  in  altitude.  A  more 
rapid  descent  in  altitude  will  cause  the  electron  density  to  decrease  more 
slowly  than  shewn  in  Figs.  12  and  13  as  discussed  above.  One  must  not 
base  the  late-time  behavior  of  the  peak  electron  density  on  the  altitude  of 
release  if  it  indeed  is  descending  in  altitude. 

E.  COMPARISON  WITH  DATA 


We  are  now  in  a  position  to  evaluate  the  results  given  by  this  model. 
Although  the  model  has  a  reasonable  physical  basis,  the  values  of  the 
parameters  used  were  chosen  to  give  agreement  with  optical  and  radar 
data.  The  first  point  to  make  is  that  the  value  of  the  barium-neutral 
collision  time,  ,  ,  given  by  Eq.  (1)  is  in  agreement  with  the  optical  and 
radar  measurements  of  the  time  variation  of  the  gaussian  radius  of  the 
neutral  cloud  and  the  time  variation  of  the  gaussian  length  of  the  ion  cloud 
parallel  to  the  magnetic  field.  In  addition  to  the  physical  modeling,  this 
parameter  is  the  single  most  important  parameter  in  obtaining  the  results 
that  have  been  given  above. 

The  optical  data  have  also  been  used  to  obtain  an  estimate  of  the 
early-time  expansion  of  the  neutral  cloud  and  estimates  of  the  most  likely 
scale  size  of  the  ion  cloud  transverse  to  the  magnetic  field.  This  estimate 
is  difficult  to  make  accurately  because  the  ion  cloud  is  optically  thick  to 
its  resonance  radiation,  which  tends  to  result  in  Loo  large  an  estimate. 
However,  by  recognizing  that  one  is  primarily  concerned  with  the  scale 
size  of  the  high  density  portion  of  the  ion  cloud,  one  could  sometimes  obtain 
the  best  estimate  of  the  appropriate  scale  size  many  minutes  after  release 
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when  the  dense  core  has  moved  away  from  the  denser  portions  of  the  sur¬ 
rounding  ion  bridge.  This  determination  at  these  later  times  can  be  compli¬ 
cated  by  the  fact  that  the  ion  cloud  has  elongated  considerably,  parallel  to 
the  magnetic  field.  Hence,  we  have  had  to  use  our  judgment  in  obtaining  the 
best  estimates  with  which  we  wish  to  match  the  model. 

The  HF  long  pulse  radar  data  provide  a  lower  bouro  on  the  peak 
electron  density.  At  early  times  when  the  ion  cloud  has  not  been  distorted 
very  much  from  a  gaussian  shape,  this  measurement  probably  provides  a 
fairly  reliable  estimate  of  the  peak  ion  density.  However,  due  to  the 
changing  aspect  of  the  ion  cloud,  the  additional  ionization  deposited  in  the 
ion  bridge,  and  the  structuring  that  takes  place  in  the  ion  cloud  at  later 
times,  the  HF  data  provide  a  less  reliable  estimate  at  later  times.  Hence, 
we  have  tended  to  place  a  larger  reliability  on  the  times  of  the  dropout  of 
the  highest  frequencies  tended  to  occur  in  the  time  range  t^  <  t  <  t^  where 
the  peak  density  given  by  Eq.  (24)  is  felt  to  be  most  reliable. 

Table  1  presents  a  comparison  of  the  data  with  the  values  of  several 
parameters  obtained  from  the  model.  The  ion-cloud  radius  is  estimated 
from  optical  data,  and  the  ion-cloud  densities  are  given  for  each  event  at 
the  time  after  release  indicated.  We  see  that  the  model  estimate  of  the 
scale  size  of  the  ion  cloud  transverse  to  the  magnetic  field,  a^,  is  in 
excellent  agreement  with  estimates  obtained  by  the  optical  data.  The  value 
of  the  peak  electron  density  obtained  from  the  HF  data  for  the  48-kg  payload 
at  190  km  represents  an  average  of  the  Dogwood,  Spruce  and  Quince  data  at 
100  s  after  release.  The  model  gives  excellent  estimates  for  these  nominal 
releases.  In  particular  the  estimate  for  the  electron  content  agrees  with 
the  value  produced  by  Event  Quince  as  measured  by  Faraday  rotation. 

Events  Dardabasi  I  and  Redwood  I  are  also  modeled  very  accurately. 
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TABLE  1.  CORRELATION  OF  MODEL  WITH 


in  the  case  of  Olive  the  returns  from  the  36  MHz  radar,  corres- 

7  -  3 

ponding  to  a  density  of  1.  6  x  10  cm  ,  dropped  out  at  660  seconds  after 

_  1/2 

release.  By  extrapolating  backward  using  a  t  '  decay  law,  one  obtains 

7  -  3 

an  estimate  of  3.  1  x  10  cm  at  three  minutes  after  release.  At  the  same 

time  the  incoherent  backscatter  radar  was  measuring  a  peak  concentration 

of  the  order  of  5  x  10  cm  .  This  latter  radar  technique  has  consistently 

tended  to  give  values  of  peak  electron  density  somewhat  higher  than  those 

indicated  by  the  HF  radar.  We  also  note  that  Olive  II  and  Olive  III  have  not 

been  included  in  the  tabulated  model  estimate.  Each  of  these  releases  can 

be  expected  to  be  similar  to  Dardabasi  at  3  minutes  after  release  and  could 

7  -  3 

thus  contribute  an  additional  1.  1  x  10  cm  to  the  indicated  value  of  the 

peak  density  for  Olive.  Since  these  two  events  did  not  occur  at  the  same 

7  -3 

altitude  as  the  main  ion  cloud,  this  latter  estimate  of  4.  5  x  10  cm 
represents  an  upper  bound  to  the  estimated  peak  density  of  Olive.  On  the 
other  hand,  the  total  electron  content  in  the  three  Olive  releases  is  one- 
third  larger  than  the  result  given  in  Table  1  for  Olive  I  alone. 

The  model  predicts  a  high  value  of  the  electron  density  for  Nutmeg  I. 
This  value  is  not  inconsistent  with  the  data.  The  36-MHz  signal  dropped 
out  at  160  seconds  after  release  and  the  next  two  lower  frequencies  dropped 

_  i  /2 

out  considerably  faster  than  at  '  law  would  predict.  However,  optical 
records  show  that  Nutmeg  I  was  definitely  striated  before  105  seconds  after 
release  and  previous  experience  has  shown  that  the  HF  radar  tends  to  under¬ 
estimate  the  ion-cloud  electron  density  when  the  cloud  is  striated.  We  sug¬ 
gest  that  the  value  given  by  the  model  is  an  accurate  estimate  of  the  electron 
density  in  Nutmeg  I. 
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The  time  histories  for  the  peak  electron  density  for  the  sample 
cloud  shown  in  Figs.  12  and  13  also  are  quite  reasonable  and  do  not  conflict 
with  any  radar  data.  In  particular  the  asymptotic  values  of  the  density 
shown  in  Figs.  10,  12,  and  13  provide  an  explanation  for  the  fact  that  HF 
returns  from  ion  clouds  last  for  very  long  times  after  release.  This  result 
indicates  that  the  scattering  at  late  time  is  indeed  from  electron  densities 
that  exceed  the  critical  density.  It  should  be  remembered  that  the  asymptotic 
values  given  in  Fig.  10  are  for  a  release  whose  peak  remains  at  the  release 
altitude.  As  the  ion  cloud  descends,  the  inventory  on  a  field  line  remains 
constant  and  the  peak  density  increases  as  H  ^  given  by  Eq.  (20)  decreases. 
The  asymptotic  value  can  also  be  obtained  by  cross  -  plotting  along  constant 
content  contours  with  the  aid  of  Fig.  6. 

The  most  important  parameter  that  is  necessary  for  describing  the 
behavior  of  barium  ion  clouds  once  they  have  been  produced  is  the  height- 
integrated  Pedersen  conductivity  given  by 


L. 
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Ba 


/  — 
J  1  +  (Hr) 


,  x  J  3.2  x  10'11  x  or  „ 

?  n  (z)  dz  -  - = -  N,  mho 

1  +  (Qrr  *■ 


(25) 


-  2  —  1 

where  N  „  is  in  cm  and  =  eB/mn  =  35  s  for  a  magnetic  field  of 
H  '  Ba  ° 

l/2  gauss.  The  variabion  of  the  cloud  conductivity  wi .  .ltitude  and  payload 
is  shown  in  P'ig.  14.  We  shall  make  use  of  these  estimates  for  the  ion-cloud 
conductivity  when  we  compare  the  theory  of  ion-cloud  dynamics  with  data 
obtained  during  the  Secede  test  series. 


F.  DISCUSSION 


There  are  a  number  of  effects  whose  variations  have  not  been 
included  in  the  above  modeling  and  that  have  an  influence  on  the  properties 
of  the  ion  cloud.  These  effects  include  the  ambient  neutral  wind  and 
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Fig.  14  The  Height-Integrated  Pedersen  Conductivity  for  a  Barium  Ion 
Cloud  as  a  Function  of  Altitude  and  Payload 


electric  field,  the  optical  thickness  of  the  neutral  cloud,  the  cloud  size  and 
inventory,  the  latitude  of  the  release,  the  time  of  day  and  the  season  of  the 
year  of  the  release,  and  the  effects  of  barium  cloud  dynamics.  We  shall 
discuss  each  of  these  briefly,  indicating  the  assumptions  and  limitations  of 
the  model  and  how  it  may  be  extended  to  include  some  of  these  effects. 

The  ions  have  been  assumed  to  remain  at  rest  after  the  barium  atoms 
have  been  ionized.  If  there  is  a  large  wind  of  the  ambient  atmosphere  at 
the  altitude  of  the  release,  or  if  there  is  a  large  ambient  ionospheric 
electric  field,  there  will  be  a  relative  velocity  between  the  neutral  barium 
atoms  which  are  being  ionized  and  the  barium  ions  thus  created.  Hence  the 
ionization  will  be  distributed  over  a  larger  area  in  the  plane  transverse  to 
the  magnetic  field.  As  a  result  the  electron  density  will  be  less  than  if  no 
neutral  wind  existed  as  assumed  in  the  model.  For  instance,  a  neutral 
wind  of  100  m/s  at  an  altitude  of  190  km  would  displace  the  neutral  cloud 
by  2.  5  km  during  the  25-s  ion-creation  time  appropriate  to  that  altitude 
given  by  Eq.  (12).  A  simple  calculation  indicates  that  the  volume  in  which 
the  ionization  is  now  distributed  has  been  increased  by  approximately  50%: 


Volume  ~ 


,  V  t 

1+4-2- 

4  ai 


~  1.5 


This  effect  probably  accounts  for  the  fact  that  the  electron  density  in  Event 
Plum  I  which  had  a  large  neutral  wind  present  was  approximately  2/3  that 
of  the  Spruce  and  Quince  ion  clouds.  During  Secede  III  the  electron  density 
of  ion  clouds  released  during  high  ambient  neutral  winds  and  electric  fields 
was  significantly  reduced  from  the  values  given  by  this  model. 
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We  have  not  treated  the  influence  that  the  thickness  of  the  neutral 
cloud  to  the  metastable-exciting  solar  radiation  has  on  the  production  of 
ions  in  the  ion  cloud.  Many  of  the  barium  atoms  on  the  dark  side  of  the 
neutral  cloud  should  be  shielded  from  this  radiation  until  the  excited  atoms 
have  been  consumed  by  photoionization  and  oxidation.  This  effect  should 
cause  a  lengthening  in  the  time  in  which  the  bulk  of  the  ionization  is  created. 
As  a  result,  the  electron  density  in  the  ion  cloud  should  be  somewhat  less 
than  that  given  by  the  model,  since  the  neutral  cloud  would  have  diffused 
to  a  larger  size  and  any  ambient  winds  would  have  further  displaced  the 
neutral  cloud.  In  spite  of  this  influence,  this  model  has  been  seen  to 
give  excellent  agreement  with  observations  of  barium  ion  clouds.  Perhaps 
the  fact  that  we  have  assumed  that  the  size  of  the  ion  cloud  is  determined 
by  the  size  of  the  neutral  cloud  at  two  ionization  time  constants  rather  than 
one  after  an  equal  mass  radius  has  been  reached,  has  partially  compensated 
for  the  optical  depth  effect.  This  effect  should  also  cause  a  more  non- 
uniform  deposition  of  the  ions  in  the  ion  cloud,  which  is  probably  of  more 
importance  for  the  initial  condition  to  be  assumed  in  developing  the  theory 
of  barium  ion  clouds  than  it  is  here,  in  giving  typical  values  for  key 
parameters. 

We  note  that  in  order  for  our  modeling  to  give  agreement  with  the 
radar  data  we  have  assumed  a  30%  vaporization  efficiency  for  the  production 


of  barium  atoms.  This  value  is  somewhat  higher  than  the  normally  accepted 
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reported  vaporization  efficiencies  approaching  45%  in  the  cases  of  Events 
Kumquat  and  Lime  when  the  fraction  of  barium  atoms  lost  to  oxidation  has 
been  accounted  for.  The  30%  figure  is  in  agreement  with  the  inventories 


reported  for  Events  Sapsucker  and  Titmouse.  By  using  beacon  data 
Prenatt^  has  developed  an  ion-cloud  model  for  Plum  that  corresponds  to  a 
20%  vaporization  efficiency.  In  view  of  the  excellent  agreement  that  this 
model  has  produced  with  the  optical  and  radar  data  we  feel  that  this  30% 
figure  is  realistic  for  the  chemical  mixtures  used  during  Secede  II. 

If  the  ion  inventory  had  corresponded  to  a  5%  to  10%  vaporization 
efficiency,  the  size  of  the  ion  cloud  transverse  to  the  magnetic  field  would 
have  to  be  smaller  than  the  value  given  by  a^  in  Eq.  (11)  in  order  that  the 
electron  densities  remain  as  high  as  reported  by  radar.  A  minimum 


possible  size  that  might  be  considered  reasonable  is 


a  =  (a  2  +  4D  t.)1/2 
m  o  n  i 


corresponding  j  the  size  of  the  neutral  cloud  at  one  ionization  time  constant 

after  the  vapor  has  reached  an  equal-mass  radius.  This  size  was  found  to 

be  between  0.  6  and  0.  7  times  a given  by  Eq.  (11)  throughout  the  altitude 

and  payload  range  covered  by  the  model.  Hence  there  is  a  possibility  that 

an  8%  vaporization  efficiency  could  be  assumed  if  the  transverse  scale  size 

for  the  Spruce  ion  cloud  was  2.  5  km  instead  of  3.  75  km.  We  note  that 

Prenatt  has  estimated  a  gaussian  radius  of  2.  8  km  for  Event  Plum  I,  and 

8 

the  incoherent  backscatter  radar  results  at  early  times  in  Event  Spruce 
indicate  a  gaussian  radius  of  2.  5  km  transverse  to  the  magnetic  field.  It 
should  be  noted  that  the  same  experiment  reports  an  effective  initial 
gaussian  length  of  4.  3  km  parallel  to  the  magnetic  field.  At  22  seconds  after 
release,  these  two  sizes  should  have  been  equal.  However,  it  is  not  likely 
that  the  Spruce  ion  cloud  was  gaussian- shaped  at  that  time  due  to  optical 
thickness  effects.  These  small  values  for  the  transverse  scale  size  of  the 
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ion  cloud  appear  to  conflict  with  the  optical  data.  It  is  fair  to  suggest  that 
the  scale  size  of  the  ion  cloud  transverse  to  the  magnetic  field  and  the  total 
ion  inventory  are  somewhat  unresolved  questions.  In  spite  of  this  uncer¬ 
tainty  it  appears  that  the  electron  densities  produced  by  the  above  model 
are  in  excellent  agreement  with  the  data. 

The  modeling  that  has  been  done  has  tacitly  assumed  midlatitude 

locations  with  magnetic  dip  angles  larger  than  45°  but  southward  of  the 

auroral  zone.  One  parameter  that  changes  as  a  function  of  altitude  is 

the  downward  drift  due  to  gravity  as  given  in  Eq.  (23).  Above  200  km, 

this  drift  is  significant  for  large  barium  ion  clouds.  The  latitude  of  the 

release  will  also  have  an  effect  on  the  limiting  value  c.f  the  density  n. 

that  is  reached  by  the  ion  cloud.  In  particular  the  value  given  by  Eq.  (21) 

should  be  decreased  by  the  factor  sin  6 d  ,  since  the  ion  cloud  can  diffuse 

parallel  to  the  magnetic  field  for  a  longer  time  before  it  has  diffused  a 

scale  height  in  eltitude.  As  far  as  high  latitude  is  concerned,  it  should  be 

noted  that  this  model  should  be  applied  with  caution  due  to  the  large  neutral 

winds  and  ambient  electric  fields  that  can  exist  in  and  near  the  auroral  zone 

and  to  the  large  temporal  and  spatial  variability  of  both  of  these  parameters. 
q 

Lins  on  has  shown  that  the  conductivity  of  the  conjugate  ionosphere 
in  the  southern  hemisphere  influences  the  motion,  distortion,  and  striation 
characteristics  of  barium  ion  clouds.  Although  the  conductivity  of  the  local 
ionosphere  will  be  more  typical  of  nighttime  values  than  of  daytime  values, 
the  conductivity  of  the  conjugate  ionosphere  can  depend  on  both  the  time  (.  f 
day  { i .  c.  ,  whether  it  is  a  morning  or  eve  "ling  twilight  relearn)  and  the 
season  of  the  year.  Hence  the  behavior  of  barium  ion  clouds  does  depend  on 
whether  the  conjugate  ionosphere  is  in  sunlight  or  in  darkness.  For  instance, 
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during  the  morning  releases  of  Secede  I  in  May,  the  conjugate  ionosphere 
was  in  darkness  while  during  the  evening  releases  of  Secede  II  in  January, 
the  conjugate  ionosphere  was  in  sunlight.  Hence  in  planning  future  tests 
both  the  time  of  day  and  the  season  of  the  year  need  to  be  considered  if 
one  is  attempting  to  make  predictions  of  the  motion  of  barium  icn  clouds 
or  the  onset  time  for  the  formation  of  striations. 

The  modeling  described  in  this  section  has  not  included  any  of 
the  dynamics  of  barium  ion  clouds  that  are  discussed  in  Chapter  5  that 
summarizes  the  relevant  theory.  In  particular  the  motion,  distortion, 
and  stria tion  behavior  of  barium  ion  clouds  depends  on  the  ratio  of  the  ion- 
cloud  height-integrated  Pedersen  conductivity  to  the  height-integrated 
Pedersen  conductivity  of  the  ionosphere  in  both  hemispheres.  In  order  to 
model  these  plasma  effects  the  conductivity  of  the  ion  cloud  can  be  esti¬ 
mated  with  the  aid  of  Fig.  14.  Once  this  value  is  known,  an  estimate  can 
be  made  of  the  expected  ionospheric  conductivity.  It  is  important  to 
remember  to  include  the  contribution  to  this  conductivity  from  the  iono¬ 
spheres  in  both  hemispheres.  In  particular,  near  the  magnetic  equator 
the  magnetic  field  lines  on  which  the  ion  clouds  are  created  never  leave  the 
ionosphere  and  an  accurate  estimate  carried  out  over  the  entire  length  of 
the  field  line  must  be  made.  Once  the  value  of  this  ratio  of  conductivities 

is  known,  then  an  estimate  can  be  made  for  the  motion  of  barium  ion 

9 

clouds. 

Once  barium  ion  clouds  have  striated  so  that  very  small  scale 
lengths  perpendicular  to  the  magnetic  field  have  been  formed,  then 
diffusion  perpendicular  to  the  magnetic  field,  which  has  been  neglected  in 
the  above  modeling  due  to  the  large  scale  sizes  involved,  becomes 
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important  and  is  the  mechanism  by  which  the  ionization  is  eventually 
dissipated.  The  time  scale  for  this  dissipation  is  also  discussed  in 
Chapter  5  describing  the  theoretical  results. 
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PARTIAL  REDISTRIBUTION  IN  RADIATION  TRANSPORT 


2.  PARTIAL  REDISTRIBUTION  IN  RADIATION  TRANSPORT 

B.  Kivel 


A.  INTRODUCTION 

Artificial  clouds  produced  by  chemical  releases  in  the  vxpper  atmos- 

12 

phere  are  of  current  interest.  ’  The  clouds  are  produced  at  twilight  so 
that  they  may  be  observed  in  the  scattered  sunlight  while  the  sky  is  dark. 
Photographic  information  allows  determination  of  cloud  position  and  spatial 
extent.  It  is  also  possible  to  estimate  the  column  densities  of  the  scatter¬ 
ing  macerial  using  calibrated  intensity  data.  However,  because  of  optical 
depth  effects,  the  relation  between  the  intensity  and  the  amount  of  material 
is  not  linear.  The  purpose  of  this  paper  is  to  assist  in  understanding  the 
relation  between  column  density  and  scattered  intensity  using  simple  models. 

We  consider  the  problem  of  scattering  by  a  slab  of  atoms  with  only 
one  resonant  transition  that  scatters  isotropically.  Of  primary  interest  is 
the  effect  of  the  partial  redistribution  of  radiation  energy  which  is  caused 
by  the  Doppler  effect.  Collisions  are  neglected,  which  is  reasonable  for 
the  high  altitudes  o*’  interest.  For  simplicity  the  effect  of  polarization  will 
also  be  omitted.  We  approach  the  problem  using  two  models.  The  first  is 
an  analytic  treatment  which  does  not  allow  for  redistribution  of  photon 
energy  in  the  scattering  event  and  the  second  is  a  Monte  Carlo  calculation 
which  allows  for  the  correct  redistribution  caused  by  the  Doppler  effect. 

In  many  radiation  transport  problems,  especially  where  collisions 

3  4 

dominate,  it  is  customary  to  use  complete  redistribution.  ’  No  redistri¬ 
bution  is  also  used,  in  particular  it  has  been  used  in  analysis  of  the  atmos- 
pneric  sodium  in  the  airglow.  The  importance  of  redistribution  is  that  it 
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allows  photons  which  start  near  the  center  of  the  line  and  have  a  small  mean 
free  path  to  change  in  energy  to  regions  where  the  mean- free- path  is  longer 
and  it  is  easier  to  escape  the  cloud.  An  analytic  solution  taking  into  account 

the  effect  of  partial  redistribution  has  been  used  for  the  interpretation  of 

/ 

spectra  lrom  an  optically  thick  nebula.  Comparison  of  the  Monte  Carlo 
with  the  analytic  results  for  no  redistribution  allows  one  to  estimate  the 
importance  of  the  redistribution  effect  in  the  upper  atmosphere  scattering 
produced  by  artificial  clouds. 

In  Section  B  we  gi're  the  analytic  solution  for  diffusion  of  monochro¬ 
matic  energy  photons  with  constant  mean- free- path  through  a  one- dimensional 
slab.  This  solution  is  well  known  and  hence  forms  the  basis  for  testing  the 
Monte  Carlo  results.  In  Section  C  we  give  the  physical  basis  and  in  Sec¬ 
tion  D  the  mathematical  details,  for  the  Monte  Carlo  calculations.  Section 
E  contains  a  comparison  of  the  results  for  both  methods  and,  in  particular, 
compaies  the  reflection  and  transmission  coefficients.  A  discussion  of  the 
results  is  given  in  Section  F. 
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B.  ANALYTIC  MODEL 


In  this  section  the  problem  of  resonant  isotropic  scattering  neglect¬ 
ing  the  effect  of  polarization  by  scattering  atoms  with  a  Maxwell- Boltzmann 
distribution  of  velocities  in  a  one- dimensional  slab  is  treated.  We  solve  the 
one-dimension  diffusion  differential  equation.  This  is  a  well  known  elemen¬ 
tary  problem  and  a  basic  discussion  of  it  may  be  found  in  standard  texts 

7 

such  as  Morse  and  Feshbach. 

We  treat  the  problem  of  an  incident  flux  normal  to  one  face  of  the 
cloud  and  vary  the  photon  energy  and  slab  thickness.  Our  primary  interest 
will  be  in  regions  in  which  the  cloud  is  optically  thick  since  the  optically 
thin  solution  which  relates  intencIty  linearly  to  the  number  of  scatterers  is 
trivial.  It  will  be  seen  from  the  results  that  in  the  optically  thick  regions 
the  radiation  is  essentially  reflected  from  the  illuminated  face.  However, 
in  the  line  wings  as  the  cloud  becomes  less  thick  a  fraction  of  the  light  is 
diffusely  transmitted  through  the  slab.  In  the  still  less  thick  regions  a 
part  is  transmitted  without  scatter  and  of  the  fraction  which  is  scattered, 
about  half  is  reflected  and  half  transmitted. 

We  start  with  a  one  dimensional  diffusion  equation  for  the  radiation 


flux  ( J ^  )  at  frequency  v  ;  i.  e.  , 


9  p 
r  v 

*  £ 


(i) 


c 


whe  re 


X 


is  the  velocity  of  light 

is  the  mean-free-path  at  frequency  v 


where  the  subscript  o  stands  for  the  solar  illuminated  side  of  the  cloud, 
and 

<5> 

where  the  subscript  r  stands  for  the  values  on  the  dark  side  of  the  slab 
of  optical  depth  normal  to  the  plane  face  of  the  cloud.  These  boundary- 
conditions  are  derived  approximately  in  Morse  and  Feshbach.  They  are 
required  to  give  an  average  of  zero  for  the  diffusing  flux  from  outside  into 
the  slab.  The  result  is  justified  by  more  rigorous  solutions  of  the  corres¬ 
ponding  integral  equation 
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The  solution  of  the  differential  equation  (Eq.  (3)),  is 


Pv 


3—z  , 

—I  e  I  a  z  +  b 
c  v 


(6) 


where  the  two  unknown  constants  a  and  b  are  fixed  by  the  boundary  condi¬ 
tions  (Eqs.  (4)  and  (5)).  In  this  way  we  find 


31 


a  =  - 


„  /1-7-  0.3e  T\ 

\  1.4  +  T  / 


and 


+  0.  7  a 


(7) 

(8) 


As  an  example  consider  the  case  with  t  =  2.  3.  A  plot  of  the 
diffusing  radiation  density  p  normalized  by  the  solar  flux  I  ^  divided  by 
the  velocity  of  light  (p  /(I^/c))  as  a  function  of  distance  in  the  slab  | 
normalized  by  the  radiation  mean-free-path  X  is  shown  in  Fig.  1.  The 
density  peaks  at  a  distance  close  to  optical  depth  of  1  (£  «  X  )•  There  is 
more  radiation  on  the  solar  illuminated  side  (|  =  0)  than  on  the  dark  side 
(£,  =  2.  3  X  ).  Also  the  gradient  of  density  is  larger  on  the  solar  illuminated 
side  and  this  implies  more  radiation  is  re-emitted  from  the  solar  side 
than  from  the  back.  The  figure  also  indicates  that  the  boundary  conditions 
are  satisfied,  since  the  curve  of  p  extrapolates  to  zero  at  a  distance 
A  %  ~  0.  7X  from  the  boundary  (i.  e.  ,  at  |  =  3.  0  f  and  at  |  =  -  0.  7X  ). 

At  higher  t  the  values  of  R  =p/(l/c)  are  larger  for  the  same 

z  and  the  peak  occurs  at  larger  z.  For  example,  at  r  =  4.  6  at  z  =  0, 

R  =  1.  5;  the  maximum  occurs  at  z  =  lne(-3l/ca),  i.  e.  ,  at  z  =  1.28  with 

R  -  2.  6.  On  the  other  hand  the  value  of  R  at  the  back  face  is  lower; 
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Fig.  1  Diffusion  solution  in  a  slab  of  optical  depth  t  =  2.  3.  This  figure 
illustrates  that  the  radiation  density  peaks  near  optical  depth  of 
about  1  and  that  the  gradient  of  density  at  each  boundary  extra¬ 
polates  to  0  at  0.  7  of  a  mean- free-path  away  from  the  slab. 
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i.  e.  ,  at  z  =  4.  6  R  =  0.  58.  These  differences  provide  the  basis  for  the 
larger  reflection  coefficient  for  the  slab  of  higher  optical  depth. 

We  can  now  express  the  emitted  radiation  resulting  from  diffusion 
outward  v. hich  is  in  proportion  to  the  gradient  of  the  radiation  density  at  the 
boundary  (see  Eq.(l)).  However,  because  of  the  boundary  conditions 
(Eqs.  (4)  and  (5))  the  re-emission  is  in  proportion  to  the  density  of  radia- 
tion  at  the  boundary.  Taking  account  of  the  fact  that  e  of  the  incident 
light  is  transmitted  without  scattering  we  have 


and 


(1-e 


-r 


(1-e 


P 


(9) 

(10) 


where,  as  before,  the  subscript  o  stands  for  solar  illuminated  side  and 
T  for  the  opposite  side.  For  r  -  2.  3,  illustrated  in  Fig.  1,  we  have 
I Q/lv  -  0.  55  and  I^/l^  =  0.  35  and  the  remaining  10%  of  the  incident  flux 
is  transmitted  unscattered. 

We  can  now  expre  •> ~  the  intensities  of  re-emission  analytically  as 
a  function  of  r  .  Defining  relative  intensities  for  reflected  (IR )  and  trans¬ 
mitted  (IR)  components,  we  have 


iR  -  I  =1  a 


and 


iT  - 


T  ~T 

--  =  a  —  e 


(ID 


(12) 


where  (compare  Eq.  (7)) 


a  = 


1.7-  0.  ^  e 
1.4  +  r 


-  T 


(13) 


-47- 


so  that  the  total 


From  Eqs.  (11)  and  (12)  we  obtain  i^+  irp  =  1  -  e 
re-emitted  light  is  equal  to  the  total  scattered  light  and  energy  is  conserved. 
In  the  limit  of  very  small  optical  depth  (r  —  0)  we  have  ip  =  i  =  r  /2 
which  is  the  expected  result  for  a  thin  slab.  Here  the  amount  of  light 
scattered  is  in  proportion  to  the  amount  of  material  and  for  isotropic  re¬ 
radiation  l/2  of  the  light  comes  out  of  each  face  of  the  slab.  In  the  limit 
of  large  optical  depth  (r  ■—  °°)  we  have  ip  —  1  and  i  j  ■—  0  so  that,  again 
as  expected  for  a  semi- infinite  non- absorbing  scattering  medium,  all  of 
the  light  is  reradiated  from  the  illuminated  face. 

These  results  can  now  be  applied  to  a  scattering  medium  with  a 
single  scattering  line  with  Gaussian  profile.  In  Fig.  2  we  give  the  rela¬ 
tive  intensity  (or  line  profile)  for  an  incident  flux  independent  of  frequency, 
as  a  function  of  frequency  normalized  to  the  Doppler  half-width  for  an 
optical  depth  at  the  line  center  of  100.  The  optical  depth  as  a  function  of 
wavelength  is 

2 

T  =  100  e  X  (14) 


where 


x  =  Av/A  Vp 


the  ratio  of  the  frequency  from  the  line  center  (A  v  )  to  the  Doppler  half 
width  (A  v  ).  The  curve  in  Fig.  2  labeled  FRONT  corresponds  to  ip 
and  BACK  to  i  j.  In  this  case,  the  light  reflected  from  the  solar  illumi¬ 
nated  face  is  about  5  times  larger  than  that  scattered  from  the  dark  side. 
The  radiation  from  the  back  face  comes  mainly  from  regions  in  the  line 
wings  where  the  optical  depth  is  about  2.  In  the  far  wings  where  the 
scattering  slab  is  optically  thin  the  intensities  are  the  same  front  and  back. 
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RELATIVE  INTENSITY 


At  the  line  center  all  light  is  reradiated  from  the  front  and  the  total  inten¬ 
sity  scattered  from  the  front  is  determined  by  the  effective  line  width  (i.  e. 
the  wavelength  region  where  t  >1). 

We  have  assumed  throughout  that  there  is  no  redistribution  of  radia 
tion  in  wavelength  because  of  the  absence  of  collisions.  This  assumption 
is  not  true  because  of  the  Doppler  effect,  where  the  recoiling  atom  may 
gain  or  lose  energy  at  the  expense  of  the  photon.  The  importance  of  this 
effect  is  explored  in  the  following  sections  using  Monte  Carlo  calculations. 
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c.  PHYSICAL  BASIS  FOR  MONTE  CARLO  CALCULATION 

The  problem  to  be  solved  is  the  scattering  of  radiation  by  a  slab 
of  low  density  gas  with  only  one  transition  and  no  line  broadening  except 
the  Doppler  effect. 

Each  scattering  event  is  incoherent  because  there  is  a  large  phase 
change  as  a  result  of  the  scattering.  ^  Even  radiation,  which  deviates  from 
the  transition  energy  of  an  atom  at  rest  by  an  amount  large  compared  to 
the  Doppler  width,  is  in  resonance  since  the  difference  in  energy  is  com¬ 
pensated  for  by  a  change  in  the  atom's  kinetic  energy.  As  a  result  we  can 
1)  treat  the  scattered  light  as  being  emitted  isotropically  from  the  scatter- 
ing  atom,  and  2)  neglect  the  interference  effects  of  the  several  scattering 
centers. 

The  correct  problem  allows  the  photon  to  diffuse,  both  in  space 
and  energy.  The  photon  can  be  absorbed  only  by  those  atoms  that  have 
the  specific  projected  velocity  along  the  photon  direction  required  by  the 
Doppler  effect.  The  velocity  of  the  atom,  of  course,  can  be  appreciably 
larger  since  it  is  only  the  projected  component  that  is  fixed.  The  re¬ 
emitted  photon  can  be  shifted  to  the  red  or  the  blue  independent  of  the  initial 
incident  photon  energy  and  depending  only  on  the  total  velocity  of  the  scatter¬ 
ing  atom  and  the  direction  of  emission  relative  to  the  atom's  velocity 
vector.  We  assume  that  the  reradiation  is  isotropic  and  pick  the  emission 
direction  at  random  and  independent  of  the  velocity  of  the  scattering  atom. 
The  energy  of  the  emitted  photon  is  determined  by  the  Doppler  effect,  using 
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the  projected  velocity  of  the  scatterer  on  the  direction  of  emission.  The 
mean-free-path  of  the  emitted  radiation  is  determined  by  the  photon  energy. 
For  a  Maxwell- Boltzmann  velocity  distribution  the  mean-free-path  has  a 
Gaussian  dependence  on  energy.  Since  the  number  of  scatterers  is  pro- 
portional  to  exp  -  (v,|  M/2kT),  where  v^  is  the  component  along  the 
direction  of  the  photon  flight  path  and  the  Doppler  shift  is  given  by  the 
relation  Ae  =  vjj  v  / c. 

Because  we  are  interested  in  a  uniform  slab  which  is  infinite  in 
area  we  are  not  concerned  with  the  motion  of  the  photons  parallel  to  the 
face  of  the  slab.  It  is  only  necessary  to  calculate  the  projected  motion  nor¬ 
mal  to  the  face.  This  calculation  is  done  by  using  the  polar  angle  between  the 
photon  direction  and  the  direction  normal  to  the  face.  The  length  of  motion 
between  scattering  is  determined  at  random  assuming  an  exponential  atten¬ 
uation  with  distance  and  a  mean-free-path  determined  by  the  atomic  cross 
section  and  the  Maxwell- Boltzmann  distribution  of  scatterers.  The  pro¬ 
jected  length  of  the  distance  between  scattering  events  is  used  to  deter¬ 
mine  the  successive  positions  relative  to  the  face  of  the  cloud.  The  photon 
is  tracked  from  initial  entrance  until  it  emerges  from  either  face  of 
the  slab.  If  the  photon  reemerges  from  the  face  it  entered,  then  it  is 
reflected.  Reemergence  from  the  rear  face  after  scattering  is  diffuse 
transmission.  It  is  also  possible,  even  in  optically  thick  cases,  that  a 
photon  traverses  the  cloud  without  scattering;  we  call  this  unscattered 
transmis  sion. 
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D.  MATHEMATICAL  DETAILS  OF  MONTE  CARLO  METHOD 


There  are  three  variables  which  determine  the  photon  state; 

1)  its  location  in  the  slab,  given  by  the  normal  distance  z,  which  goes 
from  0  at  the  face  of  initial  incidence,  to  a  variable  upper  limit  ^  ; 

2)  the  polar  angle,  0  ,  fixing  the  direction  of  motion  with  respect  to  the 
z  axis;  and  3)  the  photon  energy,  given  by  the  non-dimensional  ratio,  x, 
of  the  deviation  from  the  line  center  divided  by  the  Doppler  half-width, 

A  ~  „  =  v  v  /c,  where  (7  is  the  wave  number  value  for  the  transition, 

vq  =  sj  2kT/M  a  nominal  particle  speed  at  temperature  T,  for  mass  M, 

with  c  =  3  x  10^  cm/s  and  k  =  1.  38  x  10  ^  erg/°K.  A  typical  set  of 

initial  conditions  for  these  parameters  are  xq  =  0;  i.  e.  ,  at  the  line  center 

0  =0,  incident  normal  to  the  face,  and  z  =  0,  at  the  incident  face  of  the 

o  o 

scattering  region. 

The  mean-free-path  for  a  Maxwell- Boltzmann  distribution  of 
x2  / 

particle  speeds  is  \  -  e  /(nQ)  where  n  is  the  particle  density  and 
Q  is  the  peak  Doppler  cross  section 

O  =  nTtt  rgf/  A~  D  (15) 


whe  re 

r  =2.8x10  ^  cm 
e 

f  =  transition  number 


For  Ba  ion  clouds  at  X  =  4934A,  T  =  900°K  and  f  =  0.  35,  Q  ^  10  *^crr?. 

For  n  =  3  x  10  cm  ,  l/nQ  l/ 3  km.  For  simplicity  and  generality  we 
have  normalized  the  distance  scale  z  by  dividing  by  the  mean- free- path  at 
the  center  of  the  scattering  line  l/nC.  The  probability  of  going  a  distance  s  and 
being  scattered  in  an  interval  ds  is  dp  =  e  ds.  We  determine  s  by 
picking  p  at  random  from  a  uniform  distribution  between  0  and  1  and 
setting  s  =  lne(l/p).  In  this  way  s  covers  the  range  from  0  to  » 

— -  g 

and  has  a  probability  of  occurrence  in  proportion  to  e 

The  distance  the  photon  moves  perpendicular  to  the  slab  face  is 

2 

then  determined  for  s,x  and  the  polar  angle  according  to  A  z  =  se  cos  0  . 

The  photon  goes  forward  for  0  <  90°  .  The  distance  moved  is 

x2 

longer  in  the  line  wings,  increasing  as  e 

The  velocity  vector  of  the  scattering  atom  is  determined  so  that 

the  component  (v^  )  in  the  direction  of  the  incident  photon  is  given  by  the 

Doppler  effect;  i.  e.  ,  v^  =  c  A~/~  =  vqx.  The  perpendicular  component 

(Vj.  )  is  picked  at  random  consistent  with  a  Maxwell- Boltzmann  distribu- 

2  2 

tion;  defining  u  =  v^  /vq  we  set  u  =  ln^fl/p)  where  p  is  a  random 
number  from  a  uniform  distribution  between  0  and  1.  The  polar  angle 
between  the  initial  photon  direction  and  the  velocity  vector  is  then 
vjj  =  tan  *  f°r  positive  x,  and  ijj  =  180°  -  tan  *  —  for  negative 

x.  The  azimuthal  angle  of  the  velocity  vector  around  the  direction  of  the 
incident  photon  is  <6  ,  picked  at  random  between  0  and  360°  (rf>  =  360  p 
where  p  is  uniform  random  0  to  1). 

We  will  need  the  polar  angle  of  the  velocity  vector  v/ith  respect 
to  the  z  axis,  Q.  This  is  determined  by  the  law  of  cosines  in  spherical 
trigonometry 
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cos  9 


cos  ip  cos  @  +  sin  41  sin  0  cos  d>. 


(16) 


The  angle  of  emission  is  taken  ?s  isotropic.  We  define  the  polar 
angle  with  respect  to  the  velocity  vector  as  0  and  the  azimuthal  angle  as 
$.  The  element  of  solid  angle  is  sin  0  d©  d  $  =  dp  d  $  where 

p  s  cos  0.  We  pick  p  at  random  between  +1  and  - 1  and  $  between 
0  and  360°  (p  =  2p- 1  and  $  =  360  p  and  p  is  uniform  random  0  to  1). 

The  Doppler  shift  of  the  emitted  photon  is  determined  by  the 
velocity  component  along  the  direction  of  emission.  The  shift  i-~  given  by 

I  2  — 

the  product  of  the  scattering  atoms  velocity  wu  +  x  and  cos  0. 

The  polar  angle  of  the  emitted  photon  with  respect  to  the  z  axis 
(3)  is  determined  by  the  law  of  cosines; 

cos  0  =  cos  0  cos  9  +  sin  ©  sin  9  cos  $  .  (17) 

With  the  values  for  x  and  0  the  distance  to  the  next  scattering 
event  is  calculated  as  described  above.  The  photon  is  followed  in  this 
way  until  it  emerges  from  either  the  front  or  rear  face  of  the  slab. 


E. 


RESULTS 


Calculations  have  been  made  for  eleven  incident  photon  energies 

x  =  0,  ±  0.  5,  ±  1 ,  ±  1.  5,  ±  2  and  ±  2.  5.  All  cases  are  for  a  normal 
o 

angle  of  incidence,  0  =0,  and  all  start  at  z  =  0.  We  have  followed  the 

°  o  o 

trajectories  of  100  photons  for  each  of  these  cases.  The  calculation  was 

done  on  an  IBM- 360  series  145  using  FORTRAN  IV  and  OSMFT;  the  1100 

cases  took  less  than  10  minutes  of  computing  time.  The  results  are 

analyzed  to  give  the  fraction  of  photons  reflected  and  transmitted.  The 

dependence  on  the  initial  energy  of  incidence  and  on  the  slab  thickness  is 

reported  in  this  section.  The  energy  spectrum  of  the  reflected  photons  is 

given  for  a  uniform  spectral  flux  of  incident  photons.  We  also  correlate 

the  angular  dependence  of  the  reflected  photons. 

Comparison  of  the  Monte  Carlo  calculations  for  reflection  and 

diffuse  transmission  with  the  diffusion  solution  without  redistribution  fe>r  depth 

'S'  =  ICO  as  a  function  of  x  are  shown  in  Fig.  3.  For  x  <  1  the  optical 

depth  is  >  30  and  95%  or  more  of  the  incident  light  is  reflected  according  to 

2 

both  methods  of  calculation.  At  x  =  1.  5,  r  =  £e  w  10  and  according  to  both 
methods  the  reflected  part  has  dropped  to  about  85%,  the  rest  being  trans¬ 
mitted.  At  x  ~  2,  t  ~  2,  about  2C%  is  transmitted  unscattered  and  of  the  rest 
50%  is  reflected  and  30%  diffuse  transmitted  according  to  Monte  Carlo. 

Here  both  reflection  and  transmission  differ  from  the  analytic  solution; 

Monte  Carle  gives  slightly  more  reflection  and  less  transmission.  At 
x  >  2.  5  more  than  80%  is  unscattered  and  the  scattered  radiation  is  approxi¬ 
mately  half  reflected  and  half  scattered. 
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Coefficients  of  reflection  and  transmission  as  a  function  of  the 
normalized  incident  energy.  The  solid  lines  are  the  analytic 
solutions  without  redistribution.  The  data  symbols  are  obtained 
by  the  Monte  Carlo  method  with  the  correct  redistribution.  The 
results  arc  for  an  optical  depth  ^  =  100  and  an  incident  flux 
normal  to  the  face  of  the  scattering  slab  (9Q  =  0).  The  reflec¬ 
tion  coefficient  is  close  to  1  at  the  line  center,  drops  to  about 
l/2  at  twice  the  Doppler  half  width  where  the  optical  depth  is 
about  2  and  the  transmission  coefficient  peaks  in  the  region 
where  the  optical  depth  is  about  2. 


In  conclusion,  Doppler  redistribution  Hoes  not  allow  signifi¬ 


cant  transmission  where  there  is  large  optical  depth.  Where  the  optical 
depth  is  about  2  in  the  line  wing, the  Doppler  redistribution  estimated  from 
the  Monte  Carlo  result  implies  somewhat  less  transmission  (not  more) 
than  given  by  the  analytic  method.  However,  the  total  reflection  for  a 
uniform  input  spectrum  is  close  to  70%  in  both  cases.  The  total  reflected 
intensity  increases  with  increasing  scattering  material  because  the  increas¬ 
ing  optical  depth  leads  to  a  larger  spectral  range  where  the  optical  depth 
exceeds  2;  i.  e.  ,  the  effective  line  width  increases  with  an  increase  in  the 
column  density  of  scattering  atoms  in  the  slab. 

1 .  Dependence  on  Slab  Thickness 

Information  has  also  been  obtained  on  the  effect  of  slab  thickness. 
The  results  are  shown  in  Fig.  4.  The  analytic  results  for  the  relative 
intensity  as  a  function  of  the  cloud  thickness  in  terms  of  the  mean-free- 
path  at  the  line  center  is  shown  by  the  solid  lines.  Monte  Carlo  calculations 
for  finite  slab  thicknesses  were  also  carried  out  and  the  results  are  shown 
by  the  data  symbols.  Within  the  statistical  uncertainty  of  the  finite  samples 
the  Monte  Carlo  results  agree  with  the  analytic  method  which  has  no  re¬ 
distribution  of  photon  energy  as  a  result  of  scattering.  These  results  are 
for  an  incident  of  photon  with  energy  at  the  Doppler  half  width,  i.  e.  , 

=  1.  Thus  we  again  find  that  the  Monte  Carlo  calculation  agrees  with 
the  method  that  neglects  the  effect  of  redistribution. 

The  transmitted  component  peaks  where  the  optical  depth  is  about 
_x2 

2;  i.  e.  T  ~  e  '  ~  6  e  .  This  is  consistent  with  the  result  in  Fig.  3. 

As  a  result  the  total  amount  of  light  diffusely  transmitted  is  not  sensitive 
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to  the  optical  depth  and  hence  insensitive  to  the  amount  of  scattering 
material.  Of  course,  the  spectral  separation  of  the  bands  increases  with 
increasing  optical  depth. 

2.  Reflected  Spectrum 

The  final  energy  (x)  of  the  reflected  photons  has  been  recorded  and 
it  is  possible  to  estimate  the  spectral  profile  for  a  uniform  incident  spectral 
flux.  For  the  uniform  incident  spectrum  the  emitted  line  is  symmetric 
about  the  central  frequency.  However,  it  should  be  noted  that  the  spectrum 
depends  on  the  distribution  of  the  incident  flux.  Photons  near  the  line  center 
are  emitted  with  the  profile  approximately  like  that  given  by  the  uniform 
distribution,  while  photons  incident  at  the  line  wings  are  reradiated  in  the 
wings  where  comparable  amounts  of  red  and  blue  photons  emerge.  The  dis¬ 
tribution  for  emerging  photons  gives  slightly  more  reradiation  on  that  side 
of  the  line  center  at  which  the  incident  flux  started.  The  difference,  how¬ 
ever,  is  not  statistically  proven  by  the  limited  number  of  calculations  carried 
out  in  this  study.  For  example,  taking  all  positive  photons  incident  (i.  e.  , 

x  =  +0.  5  to  +  2.  5)  the  difference  between  the  final  x  >  0  and  x  <  0  is  less 
o 

than  5%  of  the  total  number  of  events.  A  similar  result  occurs  for  the  incident 
negative  energy  cases  (xq  <  0). 

The  total  number  of  events  are  equally  divided  between  final  x 
values  greater  than  and  less  than  0.  As  a  result  we  summarize  the  results 
using  the  absolute  value  of  x.  The  final  energies  are  divided  into  half 
integer  regions  of  x  space.  The  results  are  given  in  Table  1.  The  first 
column  gives  the  final  energy  range  of  the  reflected  photons.  The  next 
column  gives  the  fraction  expected  in  the  energy  range.  The  third  column 
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is  the  number  of  the  769  observed  reflected  photons  expected  to  have  the 
final  energy  in  the  given  range  and  the  last  column  is  the  number  observed 
in  the  Monte  Carlo  results. 

The  curve  of  the  expected  spectrum  given  by  the  diffusion  differen¬ 
tial  equation  without  redistribution  is  the  same  as  that  for  the  reflected 
component  given  in  Fig.  3.  This  follows  because  no  redistribution  requires 
that  the  final  energy  equals  the  incident  energy  (x  -  x  ). 


TABLE  1 

REFLECTED  SPECTRUM 


x  Range 

Fraction 

Analytic 

Analytic 

Number  Expected 

Monte  Carlo 
Number  Observed 

0  -  0.  5 

.256 

197 

170 

0.  5  -  1.0 

.  252 

194 

183 

1.  0  -  1.  5 

.237 

182 

157 

1.  5  -  2.  0 

.  175 

135 

164 

2.  0  -  2.  5 

.  066 

51 

85 

2.  5  -  3.  0 

.  012 

9 

10 

The  Monte  Carlo  results  give  769  reflected  photons  from  a  uniform 
distribution  of  1100  incident  photons  in  the  range  xq  =  -2.  5  to  +2.  5,  that 
is,  there  are  70%  reflected.  A  numerical  integration  of  the  analytic  result 
shows  that  the  reflection  in  the  same  range  from  xq=0  to  2.  75  is  also  approxi¬ 
mately  70%.  That  both  methods  give  the  same  average  reflection  is  seen  in 
Fig.  3  from  the  agreement  between  the  analytic  and  Monte  Carlo  results.  Be¬ 
cause  of  this  agreement  we  normalize  the  area  of  the  Monte  Carlo  histogram 
for  the  emitted  spectrum  to  the  area  under  the  analytic  curve  for  the  reflec¬ 
tion  coefficient  in  the  region  x^O  to  2.  75.  This  normalization  is  accomplished 
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by  dividing  the  number  of  events  (N  )  in  each  half  integer  range  of  x  by 
202.  Similarly,  the  estimated  error  bracket  is  given  by  (N  ±  \f N  )/2 02 . 

The  results  are  shown  in  Fig.  5.  The  solid  line  is  the  analytic  result  and 
the  data  symbols  represent  the  Monte  Carlo  histogram.  While  the  integrated 
energy  reflected  is  the  same  in  both  cases,  there  is  a  slight  redistribution 
of  the  energy.  The  Monte  Carlo  shows  that  because  of  partial  redistribution 
less  photons  are  reflected  in  the  line  center.  Of  course,  our  conclusions 
are  limited  to  the  case  calculated  where  the  optical  depth  at  the  line  center 
is  z  =100  and  photons  arrive  at  normal  incidence,  0q  =  0.  The  reduction 
is  about  15%  in  the  spectral  range  out  to  about  1.  5  times  the  Doppler  half 
width.  The  difference  appears  as  an  excess  in  the  region  between  1.  5  and 
2.  5  times  the  Doppler  half  width.  Qualitatively,  the  spectrum  for  both 
methods  is  described  by  strong  intensity,  comparable  to  the  incident  flux, 
in  the  neighborhood  where  the  optical  depth  is  2  or  greater  and  a  rapid 
decrease  in  intensity  occurs  as  the  optical  depth  goes  below  the  value  2. 

3.  Reflected  Angular  Distribution 

The  Monte  Carlo  results  give  an  angular  distribution  for  reflected 
photons  which  is  consistent  with  Lambert's  law.  In  our  case  the  photon 
distribution  [  I  (©  )  d©  ]  is  in  proportion  to  the  product  of  the  Lambert 
factor,  cos  0  and  the  solid  angle  factor  sin  0  d0  .  This  product  is  pro¬ 
portional  to  sin  20  . 

1(0)  dG  a  cos  0  sin  0  d0  =  j  sin  20  d0  (18) 

For  comparison  we  divide  the  0  range  into  six  15°  intervals. 

The  expression  sin  20  is  analytically  integrated  to  give  the  fraction  expected 
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NORMALIZED  REFLECTED  ENERGY  [x  =  A£/(v0£/c)] 


Fig.  5  The  reflected  spectrum  for  a  uniform  incident  spectrum  as  a 

function  of  the  final  energy  of  the  reflected  photons.  The  solid 
line  is  the  analytic  solution  and  the  data  symbols  give  the  Monte 
Carlo  histogram.  The  total  reflected  energy  is  the  same  for  the 
two  cases.  However,  the  Monte  Carlo  solution  is  between  10 
and  20%  lower  at  energies  less  than  1.  5  times  the  Dcopler  half 
width  and  correspondingly  higher  in  the  region  1.  5  to  2.  5  times 
the  Doppler  half  width. 


in  each  interval.  The  results  are  given  in  Table  2.  The  first  column 
gives  the  angular  range,  the  second  column  the  expected  fraction  for  a 
Lambertian  distribution,  the  third  column  is  the  number  of  reflected  events 
expected  for  the  769  observed  reflections,  and  the  last  column  gives  the 
actual  number  given  by  the  Monte  Carlo  calculation.  It  is  seen  that  the 
observed  distribution  is  well  within  the  uncertainty  for  the  finite  number  of 
events. 

In  Fig.  6  we  compare  the  results  of  the  observed  number  of  events 
(N~  )  and  the  analytic  curve  for  sin  29  .  The  number  of  events  is  normal- 
ized  by  dividing  by  201  which  is  the  number  of  events  (769)  multiplied  by 
the  angular  range  in  radians  (7r/l2).  The  error  bars  are  ±  vNq/201. 

The  distribution  function  for  diffuse  reflection  neglecting  redistri¬ 
bution  is  approximately  (4-6  cos  9  )/7  (see  Ref.  7,  p.  187  and  also  Ref.  9). 
Multiplying  th:s  distribution  by  sin  9  cos  9  for  the  solid  angle  and  projected 
area  effect  we  obtain  a  distribution  which  is  given  in  Table  2  labeled  No 
Redistribution  and  which  is  inconsistent  with  the  Monte  Carlo  results.  The 
difference  is  attributed  to  the  few  scattering  events  for  the  incident  light 
problem  treated  by  Monte  Carlo.  About  half  of  the  reflected  photons  experi¬ 
ence  10  or  less  scattering  events  and  hence  are  approximately  isotropic. 

In  an  example  of  100  photons  started  at  zq  =  50,  xq  =  0,  9q  =  90°  for 
-  100  the  photons  were  emitted  consistent  with  the  non- Lambertian 
expression  given  ac  the  start  of  this  paragraph.  In  the  central  source 
problem  about  l/2  of  the  photons  experienced  between  100  and  300  scatter¬ 
ing  events  before  emerging  from  the  slab  and  hence  diffusion  is  a  reason¬ 


able  description. 
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The  angular  distribution  of  reflected  light  which  is  incident 
normal  to  the  face  of  the  slab  with  optical  depth  100.  The 
solid  line  is  the  analytic  expression  for  a  Lambert  law  reflec¬ 
tion.  The  data  symbols  give  the  Monte  Carlo  histogram. 


TABLE  2 


REFLECTED  ANGULAR  DISTRIBUTION 


Angular  Rang 

e  Fraction 

Number  Expected 

Number  Observed 

90°  - 

105° 

0.  067 

Lambert 

51 

No 

Redistribution 

32 

64 

105°  - 

120° 

0.  183 

141 

112 

1  37 

120°  - 

1350 

0.  250 

192 

184 

185 

135°  - 

150° 

0.  250 

192 

210 

192 

150°  - 

165° 

0.  183 

141 

167 

144 

165°  - 

180° 

0.  067 

51 

64 

47 

The  re 

were  129  diffuse 

transmitted 

events  which  is 

0.  17  of  the 

amount  reflected.  This  is  close  to  the  analytic  ratio  of  about  0.  18.  The 
histogram  of  the  spectrum  of  the  diffuse  transmitted  photons  is  in  statistical 
agreement  with  the  analytic  method  and  has  two  peaks  near  x  =  ±  2.  The 
histogram  for  the  angular  distribution  does  not  have  enough  accuracy  to 
distinguish  between  the  Lambert  law  and  the  diffusion  distribution  given  in 
the  preceding  paragraph. 


F.  DISCUSSION 


A  word  is  in  order  on  the  effect  of  complete  redistribution.  Horak, 
Kerr  and  Tierney  assumed  complete  redistribution  in  their  numerical  cal¬ 
culation  for  the  scattering  of  sunlight  by  an  artificial  cloud  of  strontium 
atoms.  In  their  model  1)  the  neutral  cloud  is  spherically  symmetric  with 
a  radial  Gaussian  density  dependence,  2)  the  two  level  atom  has  a  Doppler 
(Gaussian)  resonance  transition,  and  3)  the  solar  flux  is  treated  as  con¬ 
stant  in  the  scattering  region.  Their  results  for  this  model  can  be  des¬ 
cribed  by  an  effective  line  width.  That  is  all  of  the  sunlight  in  the  spectral 
range  where  the  scattering  line  is  optically  thick  is  scattered  isotropically. 
Actually  there  is  somewhat  more  backscattered  toward  the  sun  than  away 
and  the  isotropic  approximation  is  better  at  90°  from  the  direction  of  inci¬ 
dent  solar  flux.  The  intensity  (power  per  unit  area  per  solid  angle)  at  the 
line  center  is  approximately  the  solar  flux  divided  by  4n.  The  effective 
width  of  the  line  is  the  spectral  region  in  which  the  optical  depth  through 
the  center  of  the  cloud  exceeds  1.  This  result  is  similar  to  that  obtained 
analytically  without  redistribution  and  with  the  correct  redistribution  of  the 
Monte  Carlo  method.  Physically  this  is  expected  since,  1)  the  cloud  is 
non- absorbing  and  all  of  the  incident  light  which  enters  and  is  scattered 
by  the  cloud  must  eventually  diffuse  out  and  2)  the  spectral  region  in  which 
the  cloud  is  optically  thick  is  an  effective  measure  of  the  part  of  the  solar 
fiUx  which  is  Scattered, 

In  summary,  the  total  light  reflected  by  a  purely  scattering  slab 


increases  with  optical  depth  because  of  the  increase  in  effective  width  of 
the  scattering  resonance.  This  implies  that  a  measurement  of  reflected 
intensity  may  be  used  as  a  measure  of  the  column  density  of  the  scattering 
material.  On  the  other  hand,  the  total  light  transmitted  diffusely  through 
the  siab  is  not  sensitive  to  column  density.  In  this  case,  two  bands  occur 
in  the  line  wings.  With  spectral  information  the  separation  of  these  bands 
may  be  measured  and  may  give  a  measure  of  the  effective  line  width  and 
hence  of  column  density. 

The  spectrum  determined  using  the  Monte  Carlo  method  with  the 
correct  redistribution  shows  a  small  effect  for  optical  depth  100;  1)  the 

*  » 

reflected  spectral  intensity  near  the  line  center  is  reduced  by  about  15% 

;  and  2)  this  energy  is  displaced  to  a  region  between  2.  0  and  2.  3  times  the 

Doppler  half  width  from  the  line  center. 

The  angular  distribution  of  the  reflected  radiation  is  found  to  be 
Lambertian. 
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3.  SPRUCE  CORRELATION 
S.  J.  Powers  Jr.  and  B.  Kivel 

A.  INTRODUCTION 

Considerable  attention  has  been  given  to  event  Spruce  at  R  +  14  min 

1  2 

because  of  the  coincidence  of  radar  and  optical  measurements.  ’  This 
note  addresses  two  problems  in  the  data  correlation  of  the  event.  First  the 
optical  analysis  depends  on  the  angle  at  which  the  cloud  is  viewed.  In  order 
to  take  account  of  the  angle  dependence  estimates  have  been  made  of  the 
orientation  of  the  ion  cloud  face  and  the  angle  of  solar  illumination  as  well 
as  the  view  angles.  It  is  found  that  Tyndall  views  the  face  opposite  to  the 
solar  illuminated  one,  while  the  solar  illuminated  face  is  observed  from 
Barin.  This  explains  the  higher  intensity  seen  from  Barin.  With  these 
angles  known  an  improved  analysis  of  the  optical  intensity  data  leads  to  a 
consistent  ion  column  density  from  Barin  and  Tyndall  and  agreement  with 
the  radar  determined  density.  The  second  problem  is  the  difference  in  the 
e  *  length  along  the  earth's  magnetic  field  direction  of  the  ion  cloud  given 

by  Thomson  scatter  data  and  that  given  by  both  optical  data  and  an  ambipolar 

3  - 1 

diffusion  estimate.  According  to  Thomson  the  observed  length  to  e  is 

~25  km  compared  to  the  ambipolar  diffusion  estimate  ~  40  km.  The  differ¬ 
ence  is  attributed  to  the  skew  angle  at  which  the  radar  intersects  the 
magnetic  field  line  (the  angle  between  these  two  directions  is  3.  5°).  This 
impli-  s,  for  an  elliptical  density  profile ,  a  width  of  ~  2  km.  The  difference 
between  the  2  km  result  and  the  beacon  determined  distance  of  ~5  km  may  in 
part  be  due  to  the  east-west  extent  of  the  several  sheets  that  comprise  the 
trailing  edge.  Also,  the  Thomson  measurement  is  west  of  the  trailing  edge 


* 


and  does  not  see  the  four  density  peaks  observed  in  the  beacon  data. 

Section  B  describes  the  method  and  results  of  the  triangulation  study 
which  locates  and  orients  the  ion  cloud  with  respect  to  the  sun  and  the 
measurement  sites.  Sections  C  &  D  consider  the  peak  density  and  the 
density  profile  along  the  magnetic  field.  Density  comparisons  are  given 
here.  Section  E  treats  the  problem  of  the  ion  cloud  length  along  the  magnetic 
field  and  estimates  the  width  of  the  ion  cloud  just  west  of  the  structure  in  the 
trailing  edge. 
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B.  TRIANGULATION 


In  order  to  locate  the  clouds,  and  to  a  somewhat  lesser  extent,  to 
determine  the  ion  cloud  shape  and  orientation,  a  system  of  triangulation 
was  employed.  The  application  of  triangulation  requires  definition  of  the 
same  point  in  space  as  seen  from  two  or  more  locations.  The  symmetry  of 
the  neutral  cloud  lends  itself  to  the  application  of  rigorous  triangulation 
since  the  center  of  the  cloud  can  usually  be  estimated  unambiguously  from 
any  location.  Identification  of  a  common  point  on  the  ion  cloud  as  seen 
from  two  different  locations  is  more  difficult  since  many  prominent  features 
as  seen  from  one  sight  become  confusingly  dissimi'ar  when  viewed  from 
another  site  with  different  viewing  aspect. 

Instead  of  a  rigorous  triangulation  method,  estimates  of  cloud  loca¬ 
tion  and  shape  were  made  by  projecting  photographic  point  images  from  one 
location  into  the  other  photographic  field  of  view. 

Photographic  records  of  the  Spruce  release  taken  from  Tyndall  and 
C-6  were  available  at  various  common  times  after  release.  Included  on 
each  record  were  the  neutral  cloud,  ion  cloud  and  background  stars.  Star 
coordinates  were  read  from  an  overlay  grid  on  each  photograph.  Approxi¬ 
mately  six  sets  of  coordinates  on  each  photograph  were  measured.  These 
measurements,  along  with  the  known  right  ascension  and  declination  of  each 
star,  and  refraction  corrections  were  combined  in  a  least  squares  solution 
for  determination  of  three  angles  of  photograph  orientation  with  respect  to 
the  local  vertical- north  coordinate  system  and  an  equivalent  focal  length 
{plate  scale'. 
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Having  oriented  the  photograph  we  obtain  a  vector  from  the  center 
of  the  earth  to  a  point  in  the  photograph  as  follows.  The  transformation 
matrix  M  between  the  photograph  and  the  observers  local  vertical  coordi¬ 
nate  system  can  be  formed  from  trigonometric  functions  of  the  orientation 
angles.  A  vector  from  the  observer  in  the  direction  of  a  photographic 
image  can  be  expressed  in  the  local  system  as: 


v  =  M_c 


whe  re 


M 


(c(3c  x  ~  sa  s(3s  x 
-  s(3c  x  —  sa  cps^ 
ca  sX 


-cPsx_sasPcx  ca  sp\ 

sPsy  -sacPcx  ca  cP  J 

cacx  sa  / 


a  =  elevation  angle  of  the  optical  axis 

P  =  azimuth  angle  of  the  optical  axis 

X  =  rotation  of  the  photographic  coordinate  system 

about  the  optical  axis 

s  =  sin 

c  =  cos 


x  and  y  are  photograph  measurements  and  f]  is  the  equivalent  focal 
length. 


The  position  vector  q  of  the  observation  site  in  the  earth  centered 
coordinate  system  was  formed  using  the  method  suggested  by  Chang  and 
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Oetzeh  The  rotational  transformation  R  between  the  observation  site 
loca)  vertical  system  and  .  earth  centered  coordinate  system  was  formed 
based  on  the  longitude  and  latitude  of  the  site.  The  line  of  sight  unit  vector 
in  the  direction  of  a  photographic  image  can  then  be  expressed  in  the  earth 


centered  system  as: 


e  =  R  v 


-  Srhcd 

-  s<hs9 
c<h 


C(hc9 
C(f)C  9 

Sr 4 


where 


longitude  measured  positive  east 
geodetic  latitude 


The  position  vector  _t  =  £  +  k  e  represents  the  earth  centered  position 
vector  of  a  target  that  is  k  units  in  range  from  the  observation  site. 

The  magnitude,  k,  of  the  target  vector  was  varied  and  the  corresponding 
altitude,  longitude,  and  latitude  was  computed  in  the  iterative  procedure 
suggested  by  Chang  and  Oetzel.  In  this  manner,  ground  plots  of  the  line 
of  sight  vector  at  different  altitudes  were  generated. 

The  line  of  sight  can  now  be  located  in  a  photograph  taken  at  the 
same  time  from  another  site.  Knowledge  of  the  latitude,  longitude,  and 
altitude  for  the  second  observation  site  enables  one  to  compute  its  earth 
centered  position  vector  c^.  The  line  of  sight  vector  to  the  target  as  seen 
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from  the  second  location  is  given  by: 


—2  =  -^  "  -S2 


where  the  subscript  2  is  for  the  second  site  values.  Expressed  in  the 
local  vertical  system  at  the  second  location: 


_  -1 

—2  ~  R2  -2 


1 


where  R  is  the  inverse  matrix  of  R.  This  is  expressed  in  the  second 
photographic  system  as: 


'  _  ..  -1 
—  2  ^2  —2 


When  the  third  component  of  c^  is  scaled  to  the  previously  computed 
effective  focal  length  of  the  second  photograph,  the  photographic  coordinates 
of  the  target  result: 

X2 


-2  I  y2 


,  fl„ 


The  results  of  projecting  points,  C,  G,  and  E,  noted  in  the  Tyndall 
photograph,  Fig.  1,  at  R  +  13  l/4  min,  into  the  photographic  field  of 
view  at  the  C-6  site  are  shown  in  Fig.  2.  Also  noted  is  the  projection  of 
point  P  from  the  Fort  Meyers  site  as  observed  at  R  +  14  min  into  both 
fields  of  view. 

Computation  of  Normal  Angles 

The  approximate  midpoint  of  the  sharp  edge  of  the  ion  clc.id,  desig¬ 
nated  point  C  in  Fig.  1,  was  determined  to  be  at  182  km  altitude  based 
on  the  line  of  sight  intersection  seen  in  Fig.  2.  This  agrees  with  the 


■  76- 


\ 


NEUTRAL 

CLOUD 


POINT  C 


-  ;  o 


POINT  P 


D257I 


VIEW  FROM  TYNDALL 


D2572 


earlier  results  (see  Ref.  5,  p.  91).  The  corresponding  latitude  and 

longitude  for  point  C  at  that  altitude  is  shown  in  the  ground  plot  in  Fig.  3. 

A  preliminary  estimate  of  the  altitude  for  point  C  was  185  km  and  initial 

calculations  were  carried  out  based  on  this  value. 

In  a  similar  manner,  approximate  positions  for  points  G,  E,  and 

P  were  determined.  The  earth  centered  position  vectors  t  ,  t  ,  t  , 

*  — c  — e  — g 

and  t  were  calculated  and  the  differences  formed;  t  -t  ,  t  -t  ,  and 
— p  — e  — c  -g  — c 

t_  -t_  .  Since  all  points  are  assumed  to  lie  in  the  plane  of  the  ion  cloud 
p  c 

sheet,  the  above  difference  vectors  were  assumed  in  that  plane.  These 
vectors  were  expressed  in  the  local  vertical  system  at  point  C  as; 

E  =  R_i  (t  -t  ) 

—  c  — e  —  c 

G  =  R-l  ( t  -t  ) 

—  c  '-g  — c 

P  =  R"1  (t  -t  ) 

—  c  — p  — c 

The  magnetic  field  vector  direction,  H,  was  defined  in  the  local 
vertical  system  at  point  C  by  means  of  the  dip  angle  and  declination.  The 
sharp  edge  of  the  ion  cloud  is  felt  to  be  parallel  to  H.  The  vector  product 
of  H  with  each  of  the  difference  vectors  represents  a  vector  normal  to 
the  plane  since  all  vectors  are  assumed  coplanar.  This  vector  is  designated 
V  .  For  H  we  have  assumed  dip  =  6l.  ?.  and  dec.  =  2.  6°E. 

The  line  of  sight  vector  from  point  C  to  the  sites  at  Tyndall  and 
Barin  is  expressed  in  the  local  vertical  system  at  point  C  as: 


-T  =  Rc  *  &T  -ic» 

cb--r;‘  (aB-  lct 
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Fig.  3  Ground  plot  locating  Spruce  at  R  +  13.25  min. 
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A  vector,  S,  in  the  direction  of  the  sun  from  point  C  was  formed 
from  the  computed  azimuth  and  elevation  angles  of  the  sun.  After  normaliz¬ 
ing  all  vector  magnitudes  to  unity,  the  scalar  product  of  the  line  of  sight 
vector  to  the  sun,  S,  and  the  vector  normal  to  the  plane,  V^,  was  formed. 
This  quantity  represents  the  cosine  of  the  angle,  0  g,  between  the  two 
vectors. 


cos  9g  =  S 


We  have  used  for  the  location  of  the  sun  at  R  +  14  min,  right  ascension  - 
20h  59.  83  and  declination  =  -  17°  3.  6'. 

In  a  similar  manner,  the  angles  to  Tyndall  and  Barin  (0^  and 
0  g)  were  computed.  The  positions  of  E,  G,  and  P  were  adjusted  until 
the  angle  between  the  resultant  vector  normal  to  the  plane  and  the  line  of 
sight  vectors  to  Tyndall  and  Barin  and  the  sun  gave  consistent  results. 

These  results  are  given  in  Table  1.  Angles  larger  than  90°  are  on  the  face 
opposite  to  those  with  angles  less  than  90°.  The  values  obtained  assuming 
an  altitude  for  point  C  of  185  km  are  given  in  Table  2,  which  may  be 
compared  with  the  results  in  Table  1  to  estimate  the  sensitivity  to  the  altitude 
of  point  C. 

The  approximate  orientation  of  the  ion  cloud  was  determined  in  the 
above  manner.  In  a  coordinate  system  centered  at  the  position  of  point  C, 
the  vectors  H,  along  the  magnetic  field;  ,  normal  to  the  plane  of  the 
sheet;  and  Vp  HX  are  used  to  describe  the  orientation  of  the  ion  cloud. 

The  R  14  min  line  of  sight  of  the  radar  located  at  the  C-6  site  is 
indicated  in  the  photographs  from  Tyndall  and  C-6  (Figs.  1  and  2). 


We  use 


the  values  provided  by  Roland  Tsunuda  of  SRI  of  61.  3°  elevation  and  189.  9° 
azimuth.  The  scalar  product  of  the  radar  line  of  sight  vector  with  each  of 
the  coordinate  axes  indicated  that  the  radar  line  of  sight  makes  an  angle  of 
3.  5°  with  the  magnetic  field  line.  The  intersection  of  the  radar  line  of  sight 
and  the  ion  cloud  appears  to  occur  at  a  point  ~  2  km  below  C  in  altitude  and 
~  5  km  in  the  direction. 


TABLE  1 


POINT  C  AT  1  82  km'1' 

dec 

d0B 

d9_ 

Point 

Alt.  (H) 

9S 

D 

dH 

0B 

dH 

e  T 

1 

dH 

(km) 

(°) 

(o/km) 

(°) 

(°/km) 

(°) 

(7km) 

P 

171 

118 

6.  6 

108 

3.  3 

59 

0.  7 

G 

180 

117 

1.  5 

108 

0.  8 

59 

0.  2 

E 

164 

118 

2.  0 

108 

0.  9 

59 

0.  06 

longitude  =  86° 

22.  4' 

W,  latitude  =  29° 

45.  1’  N 

TABLE  2 

POINT  C  AT  1  85  km* 

Ct0o 

d© 

d0 

s 

B 

T 

Point 

Alt.  (H) 

©c 

dH 

dH 

0  rr 

dH 

(km) 

b 

(°) 

(7km) 

(°)B 

(°/km) 

(°) 

(7km) 

P 

171 

121 

6.  4 

109 

3.  0 

58 

0.  3 

G 

180 

121 

1.  6 

109 

0.  8 

58 

0.  1 

E 

165 

121 

1.  9 

i  09 

0.  9 

59 

0.  15 

longitude  =  86° 

2  3.  1' 

W,  latitude  =  2  9° 

44.  9’  N 
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C.  PEAK  DENSITY 


With  information  on  the  angles  made  by  the  solar  flux  and  the  views 

from  Tyndall  and  Barin  to  the  face  of  the  cloud  we  can  make  a  correction  to 

the  analysis  of  column  density  in  event  Spruce  at  14  min  after  the  release. 

As  a  first  approximation  we  have  been  using  a  model  developed  for  a  spheri- 

1  2 

cal  cloud.  Modeling  for  a  slab  has  also  been  reported.  ’  In  this  section 

we  apply  the  slab  model  to  obtain  improved  estimates  of  the  peak  ion  column 

density.  The  peak  unfiltered  intensity  as  seen  from  Barin  is  6.  3  x  10  w/cm 

(see  Ref.  5,  p.  29).  For  the  spherical  model  this  gave  a  column  density  of 
1 3  ~  2  6 

1.2  x  10  cm  .  Assuming  a  cloud  width  of  5  km  the  corresponding  peak 

7  -  3 

ion  density  is  2.  4  x  10  cm  .  This  is  about  three  times  higher  than  indicated 
by  the  radar  data  (see  Ref.  1  ,  p.  34  and  Ref.  2,  Fig.  5).  At  normal  inci  ¬ 
dence  for  the  solar  flux  the  slab  case  is  about  4  times  that  estimated  for  the 
sphere  for  large  optical  depth.  The  solar  flux  of  ~  60°  to  the  normal  re¬ 
duces  the  slab  intensity  by  about  a  factor  2.  There  is  not  much  correction 
for  the  view  angle  in  the  neighborhood  of  60°  and  optical  depths  >  100.  As 
a  result  it  requires  less  column  density  to  provide  the  observed  intensity. 

For  the  unfiltered  light  at  densities  of  interest  the  intensity  increases  approxi¬ 
mately  as  the  square  root  of  the  column  density.  Accordingly,  we  estimate 

12  -2 

the  column  density  to  be  reduced  to  '  3  x  10  cm  .  Again  taking  a  slab 

6  **  3 

width  of  5  km  this  implies  a  peak  density  of  about  6x10  cm  in  agreement 
with  the  other  data. 

In  view  of  the  ' nsensitivity  of  cloud  intensity  to  column  density  on  the 
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dark  side  we  cannot  reliably  estimate  column  densities  from  the  observed 

intensity  at  Tyndall.  We  do,  however,  try  to  show  for  the  ion  density 
12  -2 

“3x10  cm  tnat  the  observed  intensity  is  consistent.  The  data  from 
Tyndall  is  for  the  spectral  line  at  4934A.  The  ratio  of  the  filtered  light  at 
4934A  to  the  total  unfiltered  light  for  a  column  density  in  the  neighborhood 
of  3x10  cm  ,  is  about  l/lO  (e.  g.  ,  Ref.  5,  p.  140).  At  the  correspond¬ 
ing  optical  depth  r  ^  30  the  intensity  on  the  back  face  should  be  about  l/2 
of  that  on  the  front  face  (Ref.  1,  p.10).  Combining  these  two  factors  we 
expect  the  observed  intensity  from  Tyndall  in  the  filtered  photograph  to  be 

about  l/20th  of  that  observed  from  barin  in  the  unfiltered  photograph.  Re- 

“9  2 

ducing  the  Barin  intensity  by  20  we  obtain  3x10  w/crri  Q  .  This  is  slightly 

-  9  2 

larger  than  the  observed  value  from  Tyndall  of  2  x  10  w/cm  Q  (Ref.  1, 
p.  27);  however,  not  inconsistent  in  view  of  the  approximate  method  used  to 
estimate  the  expected  intensity.  By  comparing  Tyndall  to  Barin  we  have 
avoided  making  the  correction  for  the  solar  and  view  angles  since  the  solar 
flux  angle  is  the  same  for  both  cases  and  the  result  is  insensitive  to  view 
angle.  The  view  angle  from  Barin  is  ~  70°  to  the  normal  compared  to 
~  60°  for  Tyndall.  While  the  reciprocal  of  the  cosine  of  the  view  angle 
differs  in  the  two  cases  by  a  factor  1.  5,  in  the  optically  thick  limit  the 
effect  of  view  angle  is  small  and  the  difference  can  be  neglected  (Ref.  1, 
p.  19). 
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D.  DENSITY  PROFILE 


In  order  to  estimate  the  ion  density  along  the  magnetic  field  at  posi¬ 
tions  away  from  the  peak  density,  we  have  prepared  the  correlation  between 
intensity  and  column  density  shown  in  Fig.  4.  The  dashed  lines  are  the  re¬ 
sults  for  spherical  clouds  for  unfiltered  and  4934A  intensities  as  a  function 
of  column  density.  The  solid  lines  are  the  slab  estimates  appropriate  to  the 
analysis  of  Spruce  data  at  14  min  after  release.  The  data  from  Darin  is  an 
unfiltered  photograph.  In  the  optically  thin  region  the  intensity  for  the  same 
column  density  is  expected  to  be  about  three  times  larger,  corresponding  to 

a  view  angle  of  ~  70°  from  the  normal  to  the  cloud  face.  The  effect  of  view 

12  -2 

angle  becomes  negligible  when  the  column  density  is  3  x  1  0  cm  .  The 
factor  4  for  the  geometric  difference  between  the  slab  and  sphere  reflection 
is  reduced  by  a  factor  2  for  the  solar  flux  angle  of  ~  60°.  In  a  similar 
manner  we  estimate  the  correlation  as  seen  on  the  backside  of  the  cloud 
from  Tyndall.  In  the  optically  thin  region  thu  view  angle  of  ~  60°  leads  to 
an  expected  factor  of  2  increase  in  intensity.  In  the  optically  thick  region 
the  intensity  is  insensitive  to  column  density.  Our  correlation  in  this  region 
is  only  indicative  since  it  is  not  calculated  with  precision. 

Using  the  correlation  in  Fig.  4  between  column  density  and  observed 
intensity  appropriate  for  the  Spruce  at  R  f  14  min  the  column  density  as  a 
function  of  location  along  the  magnetic  field  in  the  neighborhood  of  the  trailing 
edge  is  obtained.  The  results  are  shown  in  Fig.  5.  The  symbols  x  and  o 
indicate  respectively  the  data  from  Barin  in  the  unfiltered  photograph  and 
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Fig.  4  Correlation  of  cloud  intensity  vs  ion  column  density  comparing 
spherical  and  slab  models.  The  slab  model  is  appropriate  to 
Spruce  at  R  +  14  min. 


DISTANCE  ALONG  MAGNETIC  FIELD  (KM) 


Data  correlation  of  Spruce  at  R  +  14  min  giving  ion  column 
density  near  trailing  edge  as  a  function  of  relative  distance 
along  the  earth's  magnetic  field. 


from  Tyndall  in  the  filtered  photograph  at  4934A.  The  distance  along  the 
magnetic  field  is  the  abscissa  and  the  log  to  the  base  10  of  the  column  den¬ 
sity  is  the  ordinant.  The  high  altitude  region  is  to  the  rea.der's  left  and  the 
low  to  the  right  as  marked  on  the  figure.  In  the  figure  we  compare  an  ex¬ 
pression  tor  the  expected  distribution  for  ambipolar  diffusion  and  results 
from  radar  observations.  The  ambipolar  diffusion  has  a  Gaussian  distribu¬ 
tion  with  a  half  width  of  ^ (4Dt  +  r ) '  where  D  is  the  ambipolar  diffusion 

2 

coefficient  taken  to  be  0.  1  km  /s,  rQ  is  the  initial  radius  taken  to  be  4  km 

3 

and  t  is  the  time.  This  distribution  is  consistent  with  the  optical  data. 

The  beacon  data  is  also  consistent  with  the  Gaussian  distribution  in  the  high 
altitude  region  where  data  was  obtained/*  The  beacon  data  also  gives  an 
estimate  of  the  cloud  thickness  appropriate  to  the  trailing  edge  of  the  cloud 
of  about  5  km.  We  use  this  width  in  converting  the  Thomson  scatter  and  HF 
radar  data  electron  density  to  column  density.  Taking  this  width  the  radar 
data  is  in  agreement  with  the  optical  result. 

The  density  profile  from  the  Thomson  scatter  data  do«--s  not  agree 
with  the  optical  value  having  a  half  width  of  about  12  km  compared  to  19  km 
from  the  optical  data  and  ambipolar  diffusion  model.  We  postulate  that  the 
shorter  extent  of  the  Thomson  data  is  a  result  of  the  angular  offset  cf  the 
radar  line  of  sight  from  the  direction  of  the  earth's  magnetic  field.  In  the 
next  section  we  use  this  difference  to  estimate  the  width  of  the  Spruce  cloud 
in  the  region  of  the  Thomson  scatter  measurement. 
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E.  CLOUD  WIDTH  FROM  THOMSON  SCATTER  DATA 


The  Thomson  scatter  results  may  be  used  to  estimate  cloud  width. 
We  have  noted  in  comparing  the  profile  of  the  Thomson  scatter  with  the 
optical  density  that  the  Thomson  determined  length  along  the  field  direction 
is  significantly  less  than:  1)  that  determined  optically,  and  2)  that  expected 
from  Linson's  correlation.  We  attribute  this  to  the  skew  angle  at  which 
the  cloud  is  observed  by  the  Thomson  scatter  measurement. 

Using  the  center  beam  orientation  as  provided  by  R.  Tsunoda  of 
SRI  we  have  oriented  the  beam  with  respect  to  the  ion  cloud  at  R  +  14  min. 
The  ion  cloud  is  defined  by  the  following  three  vectors.  The  magnetic  field 
vector  pointing  from  the  cloud  down  toward  the  north  (H).  A  vector  per¬ 
pendicular  to  the  magnetic  field  and  pointing  from  the  center  of  the  trailing 
edge  in  the  east  toward  the  center  of  the  leading  edge  in  the  west  (V  ^ )  and 
a  mutually  perpendicular  vector  which  is  also  normal  to  the  face  of  the  ion 
slab  pointing  down  (V^).  In  this  coordinate  system  we  define  the  polar 
angle,  9  ,  with  respect  to  the  H  axis  and  the  azimuthal  angle,  A ,  in  the 
plane  defined  by  and  V  ,  wit^  the  A  =  0  value  in  the  VN  direction. 

In  terms  of  these  angles  the  radar  intersects  the  cloud  at  9=3.  5°  and 

A  =  -15.  9°.  W„  i  these  values  the  radar  beam  is  more  to  the  east  when 

* 

close  to  the  ra  ,.tr  site  than  a  magnetic  field  line  through  the  center  cf  the 
trailing  edge  of  the  ion  cloud,  while  further  west  when  away  from  the  radar 
toward  the  southern  conjugate  region.  This  is  consistent  with  the  azimuth 
of  the  radar  and  the  declination  of  the  earth's  magnetic  field.  The  radar 


-89" 


beam  is  below  the  magnetic  field  line  through  the  center  of  the  trailing 
edge  when  close  to  the  radar  in  the  north  and  above  the  magnetic  field  line 
when  away  from  the  radar  towards  the  south. 

To  estimate  the  cloud  width  we  assume  a  Gaussian  ellipsoidal  shape 
which  is  independent  of  the  location  in  the  east- west  direction  along  . 

The  major  axis  of  the  ellipse  is  taken  along  the  magnetic  field  line,  H,  and 
is  estimated  from  the  optical  data  and  from  Linson's  correlation  to  have  a 
half  length  of  about  19  km.  The  breadth  along  the  normal  to  the  cloud  face 
is  estimated  from  ihe  observed  Thomson  radar  profile  which  shows 
a  half  width  of  about  12  km  at  R  +  14  min.  The  breadth  b  is  obtained  from 

2  2 
,  /  r  cos  9  \  .  I  r  sin  9  cos  6  \ 

1-{--  )  \  5  ) 

where  r  is  the  Thomson  scatter  e  ^  range,  a  is  the  e  *  range  along 
the  earth's  magnetic  field  (19  km),  9  ,  (b,  are  the  polar  and  azimuthal 
angles  as  defined  above,  and  b  is  the  e  *  half-width  of  the  cloud  breadth 
normal  to  the  face  of  the  cloud.  This  evaluation  gives  a  value  for  the  width 
of  2  b  =*■  2  km. 

There  are  two  obvious  differences  between  Thomson  scatter  and 
beacon  experiment  results.  Thomson  scatter  measures  only  one  ion  den¬ 
sity  peak  and  a  2  km  width  while  the  beacon  experiment  measures  several 
peaks  and  a  5  km  width.  We  assume  that  this  difference  is  explained  by  the 
fact  that  the  Thomson  scatter  measurement  is  made  west  of  the  trailing  edge; 
that  is,  west  of  the  region  where  the  structure  is  observed  in  the  optical 
photographs.  (Note  that  in  Figs.  1  and  2  the  radar  line  of  sight  at  14  min 
is  used  with  photographs  at  13^  min.  At  14  min  the  radar  is  about  7  km 
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west  of  the  trailing  edge  and  slightly  further  west  than  indicated  in  the 

figures.  )  The  lack  of  structure  and  smaller  width  west  of  the  trailing  edge 

7 

region  is  consistent  with  current  models  for  ion  cloud  structure. 
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4.  USE  OF  PHOTOGRAPHIC  DATA  FOR  ESTIMATING 
ION  DENSIT  »  AND  INVENTORY 

B.  Kivel 


A.  INTRODUCTION 

Photographic  data  give  important  information  on  cloud  shape,  growth 
and  location;  in  addition,  when  calibrated  to  give  absolute  intensity,  photo¬ 
graphic  analysis  can  be  used  in  a  more  quantitative  way  co  obtain  estimates 
of  ion  density  and  ion  column  density.  Unfortunately  because  of  many  un¬ 
knowns,  which  we  will  discuss,  it  is  not  at  pres  3nt  possible  to  make  an  ab¬ 
solute  density  determination  independent  of  other  measurements.  The  in¬ 
tensity  must  be  normalized  by  radar  data.  The  radar  normalization  of  the 
relation  between  ion  column  density  and  photographic  intensity  is  less 
complicated  at  short  times  after  the  chemical  release  when  the  ion  cloud 
has  a  symmetric  shape.  Once  a  semi-empirical  relation  has  been  established, 
it  is  possible  to  use  the  result  at  other  times  and  places  in  the  same  event 
and,  even  more  important,  for  other  events. 

One  aspect  of  the  quantitative  use  of  calibrated  film  data  that  has  not 
been  pursued  is  its  use  in  regard  to  the  triangulation  problem.  It  is  diffi¬ 
cult  to  find  a  cloud  location  by  triangulation  because  the  cloud  is  transparent, 
not  symmetric,  and  data  from  different  sites  often  have  different  sensitivity 
and  as  a  result  it  is  difficult  to  find  points  of  coincidence  in  diverse  pictures. 
This  problem  may  be  ameliorated  using  data  calibrated  in  terms  of  ion 
column  density. 


The  remainder  of  this  chapter  is  divided  into  four  parts.  The  un¬ 
certainties  in  calibrating  and  interpreting  film  data  are  discussed  in  Sec. 

B.  Section  C  deals  with  data  in  the  optically  thin  limit.  In  this  limit  effects 
of  optical  depth  are  not  important  and  the  column  density  is  linearly  related 
to  the  measured  intensity.  This  relationship  is  particularly  useful  for 
application  at  the  cloud  edge  and  at  late  times.  Section  D  deals  with  the 
optically  thick  limit  where  an  effective  line  width  method  for  semi-empirical 
correlation  of  ion  column  density  and  film  intensity  is  presented.  Also  in 
Sec.  D  the  effects  of  cloud  shape,  the  angle  at  which  the  cloud  is  viewed 
and  the  angle  at  which  the  solar  flux  impinges  on  the  cloud  are  discussed. 
Section  E  deals  with  applications  to  data  from  releases  in  SECEDE  I  and 
Spruce  in  SECEDE  II. 


B.  DIFFICULTIES  WITH  OPTICAL  ANALYSIS 


The  difficulties  can  be  divided  into  two  classes:  1)  data  reduction, 
and  2)  interpretation  of  the  reduced  data. 

Film  calibration  to  obtain  the  absolute  intensity  of  *Le  source  is 
generally  a  difficult  problem;  however,  with  care,  consistent  results  can 
be  obtained  with  accuracy  generally  better  than  a  factor  2. 

SECEDE  observation  has.  several  specific  problems.  These  include 
reciprocity,  film  temperature,  and  atmospheric  transmission.  The  recip¬ 
rocity  effect  is  the  sensitivity  of  the  film  to  the  exposure  time.  To  correct 
for  this  effect  it  is  necessary  that  the  calibration  be  applied  at  the  same 
exposure  time.  Film  sensitivity  also  depends  on  film  temperature.  In 
SECEDE  II,  although  the  experiments  were  done  in  Florida,  the  tempera¬ 
ture  of  the  ooserving  outdoor  cameras  was  ~  30°F  while  the  calibration 
standard  was  put  on  the  film  at  an  indoor  temperature.  The  colder  film  is 
apparently  more  sensitive  and  a  correction  for  this  effect  is  required. 
Atmospheric  transmission  reduces  the  signal  from  a  barium  cloud  but  does 
not  change  the  signal  from  a  standard  lamp.  Because  the  ground  haze  in 
the  Eglin  area  was  particularly  heavy  at  times,  transmission  corrections 
should  be  made.  One  method,  which  corrects  all  three  of  these  effects 
simultaneously,  is  the  use  of  known  standard  star  intensities  within  the 
photographic  data  to  normalize  the  intensity  calibration.  The  stars  are 
exposed  on  the  film  with  the  same  exposure  time,  the  same  film  tempera¬ 
ture,  and  having  the  same  elevation  angle  as  the  observed  ion  cloud  and 
hence  the  same  attenuation  correction  for  atmospheric  haze.  There  are 
two  limitations  on  the  use  of  star  calibration  techniques:  1)  stars  are  not 
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evident  in  filtered  photographs  which  look  at  selected  ion  lines,  and  2)  the 
unfiltered  photographs  which  show  stars  are  not  useful  at  early  time  after 
release  when  the  neutral  and  ion  clouds  overlap.  These  two  problems 
may  be  alleviated  by  combining  several  data  films.  Unfiltered  late  time 
photographs  when  stars  are  available  and  the  ion  cloud  is  separated  from 
the  neutral  cloud  may  provide  a  standard  to  calibrate  the  intensity  in  terms 
of  the  star  intensities.  With  this  information  filtered  photographs  at  the 
same  time  can  be  normalized  by  comparing  ion  cloud  intensities.  This 
calibration  information  can  also  be  used  at  earlier  times.  The  reciprocity 
correction  can  be  achieved  by  comparing  photographs  at  essentially  the 
same  time  with  different  exposure  times. 

Calibrated  intensity  information  does  not  immediately  give  ion  cloud 
density.  There  are  several  unknowns  which  prevent  this  correlation  being 
made  on  an  absolute  basis.  A  basic  requirement  is  the  solar  flux  spectral 
distribution.  Because  of  the  strong  Fraunhofer  effect  the  spectrum  is  far 
from  uniform  in  the  neighborhood  of  the  barium  ion  lines.  Since  it  is 
possible  to  measure  the  spectral  flux  from  the  sun  with  good  resolution, 
it  is  recommended  that  these  measurements  be  rr.aae.  One  must  also 
consider  sunlight  reflected  from  the  earth.  Another  fundamental  limita¬ 
tion,  not  as  easily  removed,  is  the  difficult  problem  of  radiation  transport. 
Barium  has  odd  and  even  mass  isotopes.  The  odd  mass  isotopes  have 
resonant  scattering  transitions  split  into  several  spectral  components, 
ihis  complicates  the  spectrum  and  allows  photons  to  travel  with  longer 
mean-free-paths  than  in  adjacent  spectral  regions.  Another  complication 
of  the  radiation  transport  problem  is  the  occurrence  of  metastable  levels 
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in  the  barium  ion.  The  population  of  the  metastable  levels  depends  on  the 
radiation  density  in  the  neighborhood  of  the  scattering  ion,  which  in  turn 
is  determined  by  the  metastable  level  population.  Therefore,  the  atomic 
population  and  radiation  transport  are  not  uncoupled  problems.  Also  this 
problem  is  not  a  local  one  since  the  radiation  density  at  the  ion  is  deter¬ 
mined  by  radiation  emission  throughout  the  cloud.  Radiation  transport 
approximations  have  been  made  for  simple  geometric  shapes  and  will  be 
discussed  in  Sec.  D.  The  correct  solution  requires  a  complicated  itera¬ 
tion  between  guessing  cloud  geometry  and  ion  density  distribution,  solving 
the  radiation  transport  problem,  comparing  the  predicted  intensity  from 
the  radiation  transport  solution  with  the  observed  photographic  intensities, 
then  correcting  the  geometric  model  in  order  to  correct  for  the  discrepancies 
between  the  predictions  and  observations  and  then  repeating  the  process  until 
a  satisfactory  agreement  between  prediction  and  observation  is  obtained. 

One  aspect  of  the  geometry  that  has  to  be  handled  even  for  simple  methods 
of  analysis  is  the  dependence  of  intensity  on  view  angle  and  solar  flux 
angle.  Another  complication  is  the  effect  of  coilisions  which  can  reduce 
the  metastable  level  population.  This  latter  complication  becomes  more 
important  at  lower  altitudes. 

As  a  result  of  the  difficulties  enumerated  above,  it  is  concluded 
that  optics  is  an  uncertain  source  for  absolute  ion  density  determination. 

Its  basic  value  is  for  interpolating  and  extrapolating  from  radar  deter¬ 
mined  values  of  peak  ion  density  or  column  density. 
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C.  OPTICALLY  THIN  LIMIT 


In  the  optically  thin  limit  the  intensity  is  proportional  to  the  column 
density.  In  this  region  the  cloud  is  essentially  transparent  and  only  single 
scattering  events  need  be  considered.  After  the  solar  flux  has  been  deter¬ 
mined,  this  method  can  be  used  on  an  absolute  basis  independent  of  radar 
measurements.  The  difficulty  of  determining  the  metastable  populations 
is  simpler  in  the  thin  than  the  thick  limit  because  the  transition  rates  are 
determined  by  the  solar  flux  independent  of  the  cloud  properties. 

The  optically  thin  limit  is  useful  in  analyzing  the  outer  edge  of  ion 
clouds  at  early  times  even  though  the  center  is  optically  thick.  The  ap¬ 
proach  is  to  determine  the  column  density  at  the  edge  where  the  cloud  is 
optically  thin  and  then  extrapolate  to  the  cloud  center  using  a  Gaussian 
density  distribution.  Another  application  of  the  optically  thin  limit  is  at 
late  times  when  the  density  has  gotten  sufficiently  low  because  of  cloud 
growth.  This  limit  may  be  determined  by  taking  the  integral  of  the  inten¬ 
sity  over  the  area  of  the  cloud.  At  early  times  this  integral  increases 
with  time  because  of  the  increasing  cloud  area;  however,  when  the  cloud 
becomes  optically  thin  the  integral  should  remain  constant  corresponding 
to  the  constant  total  ion  inventory  in  the  cloud. 
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D.  OPTICALLY  THICK  LIMIT 


In  the  optically  thick  limit  the  cloud  is  not  transparent  and  multiple 
scattering  is  important.  In  this  case  atoms  on  the  observer's  side  of  the 
cloud  block  radiation  from  atoms  on  the  far  side.  The  effect  of  the  distant 
atoms  is  only  seen  in  the  wings  of  the  scattering  line  where  the  cloud  be¬ 
comes  transparent.  In  this  case  the  total  intensity  does  not  increase 
linearly  with  column  density  but  only  as  fast  as  the  effective  spectral  line 
width  increases.  The  effective  spectral  line  width  is  determined  by  that 
spectral  region  of  the  scattering  line  which  is  optically  thick. 

In  analyzing  intensities  for  optically  thick  clouds,  two  cloud  shapes 
should  be  distinguished.  At  early  times  the  cloud  is  spherical  or  ellip¬ 
soidal  and  at  later  time  after  release  a  slab  shape  may  be  more  appro¬ 
priate.  The  elliptical  cloud  radiation  is  essentially  isotropic  although 
slightly  more  radiation  is  reflected  toward  the  sun  than  away  from  it  and 
this  anisotropy  increases  with  optical  depth.  The  flat  cloud  gives  a  large 
reflection  from  the  solar  illuminated  side  while  the  total  transmitted 
intensity  is  only  ~20%  of  that  reflected.  The  reflected  intensity  goes 
up  with  the  effective  spectral  widrh  and  is  about  four  times  the  intensity 
from  an  ellipsoidal  cloud.  The  transmitted  intensity  occurs  in  the  wings 
of  the  scattering  line  where  the  optical  depth  is  about  2;  i.  e.  ,  the  cloud 
thickness  is  about  twice  the  radiation  mean-free-path.  The  total  trans¬ 
mitted  light  intensity  is  insensitive  to  column  density;  however,  the 
spectral  shape  does  change  and  the  separation  between  the  two  bands  in 
the  wings  of  the  line  increases  with  optical  depth.  This  separation  cannot 
be  measured  photographically  but  may  be  observed  with  an  interferometer. 


Another  aspect  of  the  slab  type  cloud  is  the  effect  of  the  angle  at 
which  the  cloud  is  viewed  and  the  angle  at  which  the  sun  illuminates  the 
cloud.  In  the  optically  thin  limit  the  intensity  increases  inversely  as  the 
cosine  of  the  angle  of  view  made  with  the  normal.  This  increase  occurs 
because  more  atoms  are  viewed  through  the  slab  at  a  slant  angle  that  at 
the  normal  direction.  On  the  other  hand,  in  the  thin  limit  the  intensity  is 
independent  of  the  solar  angle  since  all  atoms  receive  the  same  solar  flux. 
In  the  optically  thick  limit  on  the  reflected  side  the  intensity  is  independent 
of  angle  since  the  surface  is  approximately  a  Lambert  law  scatterer.  In 
1  this  case,  however,  the  solar  flux  angle  enters  as  the  reciprocal  of  the 

*  cosine  to  the  normal  since  the  solar  flux  is  spread  over  a  larger  area  of 

the  slab  face. 
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E.  APPLICATIONS 


In  this  section  we  briefly  comment  on  data  from  SECEDE  I  and  II. 

The  SECEDE  I  releases  are  used  to  illustrate  the  method  of  adjusting  the 
correlation  between  ion  density  and  intensity  using  the  HF  radar  determined 
peak  electron  density.  The  Spruce  cloud  is  of  interest  because  it  has  a 
significant  body  of  radar,  photographic  and  optical  interferometer  data. 

In  addition,  the  Spruce  cloud  at  14  minutes  after  release  may  be  described 
as  a  slab  and  the  effects  of  solar  and  view  angle  are  important. 

In  Fig.  1  we  show  a  correlation  of  intensity  as  a  function  of  column 
density  for  Apple  and  Dogwood  data  for  X  =  4554A  Ball  resonance  line. 

The  numbeis  under  the  data  points  indicate  the  time  after  release  in  minutes. 
For  these  early  time  data,  the  elliptical  cloud  model  is  reasonable.  The 
solid  line  in  the  figure  is  a  semi-empirical  model  which  is  linear  at  lower 
column  density  and  levels  off  following  an  effective  width  at  high  optical 
depth.  The  magnitude  of  the  line  is  l  ojusted  to  agree  with  the  data  points. 

The  intensity  of  the  data  points  is  determined  photographically  and  the  column 
density  is  the  product  of  the  peak  radar  determined  electron  density  and  the 
cloud  width  for  a  Gaussian  shaped  cloud  based  on  the  optical  data. 

A  more  complete  correlation  is  given  in  Fig.  2.  Correlations  for 
sphere  and  slab  models  are  given  for  both  unfiltered  photographic  data  and 
for  light  at  X  =  4934A  Ball  resonance  line.  The  dashed  lines  are  for  the 
spherical  model.  The  solid  lines  are  for  observations  from  Barin  and 
Tyndall  as  indicated.  Barin  was  on  the  solar  illuminated  side  and  views 
the  cloud  at  70  degrees  to  the  normal.  The  solar  ilux  is  incident  on  the 
cloud  making  an  angle  of  60  degrees  to  the  normal.  Tyndall  is  on  the  back 
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Fig. 


Correlation  of  photographic  data  from  SECEDE  I. 
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face  of  the  cloud  and  has  a  view  angle  of  60  degrees  to  the  normal.  Here 
again  we  see  the  effect  of  optical  depth,  the  intensity  starts  linearly  with 
density  for  the  thin  cloud  and  goes  up  more  slowly  when  the  cloud  is  optic¬ 
ally  thick.  On  the  backside  of  the  cloud  the  intensity  is  very  insensitive  to 
column  density. 

In  Fig.  3  we  give  a  correlation  for  Spruce  at  R  +  14  min  comparing 
column  density  as  a  function  of  distance  along  the  magnetic  field.  The 
comparison  is  made  in  the  striated  region  on  the  eastern  part  of  the  ion 
cloud.  The  photographic  data  are  compared  with  results  obtained  on  the 
Beacon  experiment  and  fi  om  HF  and  Thomson  scatter  radar  measurements. 
Also  shown  is  a  Gaussian  correlation  given  as  a  solid  line  in  the  figure 
based  on  a  diffusion  coefficient  of  D  =  0.  1  km  /s  and  an  initial  radius 
rQ  =  4  km.  The  photographic  data  are  given  by  the  symbol  x  for  data 
from  Barin  in  the  unfiltered  photograph  and  o  for  data  from  Tyndall  at 
X  =  4934A.  The  HF  radar  and  Thomson  scatter  radar  data  give  peak 
densities.  These  are  converted  to  column  density  using  an  effective  cloud 
width  of  5  km  determined  from  the  Beacon  experiment.  The  data  from 
Barin  being  on  the  solar  illuminated  side  is  the  more  reliable  photographic 
data.  The  rgreement  of  the  data  from  Tyndall  in  the  low  density  extremi¬ 
ties  is  expected  because  the  cloud  approaches  the  optically  thin  limit; 
hcwever,  the  agreement  in  the  central  region  of  the  cloud  is  surprising  in 
view  of  the  insensitivity  of  intensity  to  column  density.  The  interferometric 
intensity  determination  made  near  the  peak  of  the  cloud  from  the  same  side 
as  Tyndall  gives  comparable  intensity  to  the  photographically  determined 
value  and  hence  agree  with  the  data  shown  in  Fig.  3. 
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Fig.  3 


Correlation  of  data  obtained  on  Spruce  at  14  min  after 
release. 


In  conclusion,  when  treated  wi^h  care  and  used  in  a  semi-empirical 
manner,  photographic  data  measures  ion  cloud  column  density. 
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5.  THEORY  OF  1 ON  CLOUT)  DYNAMICS 
AND  MORPHOLOGY 


Lewis  M.  Linson  and  Gerald  Meltz  ’ 

A.  INTRODUCTION 

There  are  three  principal  reasons  for  acquiring  a  theoretical 
understanding  of  the  complex  phenomena  that  occur  when  barium  ion  clouds 
are  released  in  the  ionosphere.  First,  so  that  a  proper  interpretation  of 
the  data  can  be  made;  second,  so  that  future  test  planning  can  be  successfully 
carried  out;  and  third,  so  that  reliable  theoretical  extrapolations  can  be  made 
to  other  environments  where  the  conditions  are  different  from  those  under 
which  barium  ion  clouds  have  been  produced. 

We  will  shov  that  a  theoretical  understanding  coupled  with  an  empirical 
scaling  of  the  Secede  data  provides  a  basis  for  making  reliable  predictions  of 
the  motion  of  barium  ion  clouds  and  the  expected  onse..  time  of  striations 
for  future  test  planning.  To  extend  our  predictions  to  other  environments, 
it  is  necessary  to  comprehend  the  physical  mechanisms  that  determine  the 
onset  time  of  striatiors,  the  spectrum  of  the  spatial  irregularities  produced 
by  the  ionization,  and  the  time  during  which  one  could  reliably  expect  the 
structured  ionization  to  exist.  A  knowledge  of  these  factors  will  also  be 
useful  in  arriving  at  a  model  of  ion-cloud/radio-wave  interactions. 
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For  these  reasons,  we  shall  present  a  rather  complete  but  somewhat 
cursory  review  of  the  status  of  our  theoretical  understanding  of  barium  ion 
cloud  phenomenology.  This  section  will  not  review  all  the  work  that  has 
been  accomplished  during  project  Secede  nor  will  it  present  lengthy  derivations 
of  the  theoretical  results.  Instead  we  will  give  references  to  the  original 
papers  and  reports  that  contain  complete  discussions  of  the  assumptions 
and  restrictions.  Our  intent  is  to  limit  the  scope  to  contributions  that  are 
relevant  to  the  Secede  barium  release  test  series.  We  underscore  the  fact 
that  our  present  theoretical  understanding  has  grown  out  of  the  contributions 
of  a  number  of  people  and  has  been  sharpened  by  the  critical  comments  of  our 
colleagues. 

In  Section  B  we  will  briefly  review  barium-ion- cloud  phenomenology 
and  in  Section  C  present  the  model  that  has  been  the  basis  of  much  of  the 
theoretical  work.  The  principal  large-scale  effect  of  the  three-dimensional 
equations  is  to  lead  to  ambipolar  diffusion  parallel  to  the  magnetic  field. 

The  assumptions  necessary  to  make  the  transition  to  a  tv/o-dimensional 
model  are  reviewed.  The  significance  of  the  distinction  between  low  and 
high  cloud  conductivity,  and  the  current  flow  patterns,  due  to  the  presence 
of  the  cloud  and  connecting  to  other  ionospheric  layers  by  field-aligned 
currents,  are  discussed. 

Section  D  discusses  two  aspects  of  the  large  scale  behavior  of 
barium  ion  clouds:  their  motion  perpendicular  to  the  magnetic  field  and 
the  several  stages  of  deformation  that  ion  clouds  undergo.  These  stages 
include  steepening  of  the  "backside",  elongation  in  the  direction  of  motion, 
and  narrowing  transverse  to  the  direction  of  motion.  The  configuration 
of  ionization  at  this  stage  is  best  approximated  as  a  thin  flat  sheet  elongated 
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Fig.  1 


Event  Titmouse  at  Three  Stages  of  Development 


along  the  magnetic  field  and  in  the  direction  of  the  neutral  wind.  In 
Section  E  we  review  the  various  aspects  of  the  fine  scale  structuring  of  the 
ion- cloud  density  that  occurs  on  the  unstable  backside  of  the  ion  cloud  and  has 
the  appearance  of  striations.  These  aspects  include  the  onset  time  of 
striations,  the  spectrum  of  scale  sizes  produced,  the  nonlinear  evolution 
of  the  field-aligned  sheets,  and  the  eventual  dissipation  of  ionization.  Section 
F  concludes  with  a  brief  summary  of  the  principal  points  established  in 
this  chapter  and  a  discussion  of  areas  where  future  work  is  needed.  Several 
key  questions  remain  for  which  present  theoretical  understanding  has  not 
been  developed  sufficiently  to  provide  definitive  answers.  We  stress  that 
the  reliability  of  predictions  depends  on  how  well  some  of  the  basic  back¬ 
ground  parameters  are  known,  e.g.  ,  neutral  wind  velocity,  ambient  electric 
field,  and  the  ionospheric  conductivity. 

B.  BARIUM- iON-CLOUD  PHENOMENOLOGY 

Figure  1  shows  a  sequence  of  views  of  Event  Titmouse  from  a 
direction  nearly  perpendicular  to  the  magnetic  field.  *  The  picture  on  the 
left  shows  that  the  ion  cloud  i\  more  elongated  parallel  to  the  magnetic  field 
in  comparison  with  the  spherical  neutral  cloud.  The  fact  that  the  ions  extend 
further  along  the  field  than  the  neutrals  is  due  to  the  electron  pressure  adding 
to  the  ion  pressure.  The  result  is  tne  familiar  ambipolar  diffusion  of  the 
ion  cloud  in  the  direction  of  the  magnetic  field.  The  dimension  of  the  ion 
cloud  transverse  to  the  magnetic  field  is  smaller  than  the  neutral  cloud 
because  the  perpendicular  diffusion  of  electrons  and  ions  is  negligibly 
small.  In  the  middle  picture  we  see  that  the  ion  cloud  has  separated  from 
the  neutral  cloud  and  that  the  ion  cloud  has  developed  a  large  density  gradient 
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on  its  "backside";  that  is, the  side  closest  to  the  neutral  cloud  out  of  which 
the  neutrals  flow.  On  the  right  side  of  Fig.  1  we  see  that  the  backside 
has  structured  and  exhibits  magnetic-field-aligned  striations. 

It  is  important  to  recognize,  when  viewing  photographs  of  barium  ion 
clouds,  that  one  is  seeing  only  a  projection  of  i  diffuse  distribution  of  ioni¬ 
zation,  very  much  elongated  in  the  direction  of  the  magnetic  field.  The 
extent  of  this  elongation  can  be  typically  as  much  as  l/3  or  l/4  of  the  distance 
from  the  viewer.  No  photograph  of  a  barium  ion  cloud  can  give  a  true  picture 
of  the  shape  of  the  cloud  in  the  plane  transverse  to  the  magnetic  t lei d  because 
only  one  line  of  sight  is  directly  up  the  magnetic  field  line.  A  striation  that 
is  only  200  meters  wide  will  appear  twice  as  wide  on  a  photograph  if  it  is 
located  as  little  as  only  l/5  of  a  degree,  or  600  meters,,  away  from  the 
magnetic  field  zenith.  We  are  fortunate,  as  we  shall  see  later,  that  three  of 
the  barium  ion  clouds  released  during  the  Secede  II  test  series  passed  through 
the  magnetic  zenith  at  a  time  when  they  were  fully  striated.  These  unique 
views  have  greatly  enhanced  our  confidence  in  the  validity  of  our  theoretical 
understanding  of  the  complex  ion  cloud  phenomena. 

We  have  stressed  the  fact  that  much  of  what  we  discuss  depends 
on  a  number  of  ambient  parameters  and  on  the  size  of  the  barium  ion 
cloud.  Of  primary  importance  are  the  ambient  electric  field,  the  neutral 
wind  velocity,  and  the  conductivity  of  the  ambient  ionosphere.  The  Secede 
barium  ion  clouds  were  predominantly  very  large  ion  clouds  in  the  sense 
that  their  conductivity  far  exceeded  the  height-integrated  Pedersen  con¬ 
ductivity  of  the  ambient  ionosphere.  Indeed,  several  of  the  ion  clouds 
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contained  more  electrons  on  magnetic  field  lines  than  the  ambient  ionosphere 
contained. 

Figure  2  illustrates  the  different  appearance  an  ion  cloud  may  have 
as  seen  from  two  different  sites.  It  presents  a  pair  of  views  of  the  Spruce 
ion  cloud  as  seen  from  Site  C-6  on  the  left  and  from  Tyndall  on  the  right. 

At  this  time,  the  barium  ion  cloud  has  elongated  approximately  60  km 
parallel  to  the  magnetic  field  and  40  km  in  the  direction  of  the  neutral  wind, 
but  has  narrowed  in  the  direction  transverse  to  its  motion  as  seen  in  the  view 
from  Site  C-6  on  the  left.  The  bulk  of  the  unstructured  ionization  is  indeed 
representable  as  a  flat  sheet.  There  is  a  complication  due  to  the  late  time 
ionization  deposited  by  the  neutral  cloud.  The  structured  ionization  represents 
a  series  of  smaller  sheets  and  rods  also  very  much  elongated  in  the  direction 
of  the  magnetic  field. 

We  now  turn  to  a  description  of  the  model  that  has  been  the  basis 
for  most  of  the  theoretical  work. 

C.  BASIC  MODEL 

1 .  Formulation  of  the  Problem 

Ion  cloud  drift,  deformation  and  structuring  are  all  determined  by 
the  c.elicate  balance  that  exists  between  drag  forces,  due  to  ion-neutral 
collisions,  and  magnetic  forces  that  are  generated  by  Pedersen  currents 
driven  by  the  dynamo  action  of  the  v'ind-blown  barium  plasma.  These 
currents  which  flow  across  the  cloud,  normal  to  the  magnetic  field,  depend 
on  the  ratio  of  the  height-integrated  perpendicular  (Pedersen)  conductivity 
within  the  cloud  to  the  background  height- integrated  ionospheric  conductivity. 
This  ratio  determines  to  what  extent  the  dynamo  fields  will  be  reduced,  or 
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screened,  by  polarization  of  the  cloud.  A  highly  conducting  background 
will  drain  any  charge  accumulations  through  parallel  electron  currents; 
a  poorly  conducting  ionosphe?  e  will  permit  polarization  charges  to  accumulate 
anc  thus  greatly  reduce  the  difference  between  the  velocity  of  neutrals,  V  , 
and  barium  ions,  V^.  Variations  in  the  current  flow  within  the  inhomogeneous 
cloud  will  result  in  relative  motion  of  the  dense  regions  and  hence  deformation 
of  the  initial  ion  density  distribution. 

A  low  density  ion  cloud  will  be  unscreened  by  polarization  fields  and 
will  stand  still  or  drift  with  the  ionosphere  as  viewed  in  the  earth's  frame, 
and  essentially  lnaintain  its  initial  distribution  by  expanding  in  a  self¬ 
similar  manner  th-ough  ambipolar  diffusion  both  parallel  and  perpendicular 

to  the  magnetic  field.  Regardless  of  the  density  in  the  ion  cloud,  steady 

2 

state  never  is  achieved  nor  is  it  theoretically  realizable.  The  time  scales 
for  most  particle  interactions,  or  collisions,  are  generally  much  shorter 
than  the  time  scales  for  deformation  and  stria tion  formation  and  development. 
As  a  result,  inertial  effects  are  seldom  of  importance  in  describing  ion 
cloud  behavior. 

Drift  and  deformation  are  large  scale  phenomena  that  can  be 
described  in  relatively  simple  terms  as  compared  to  striation  formation 
and  morphology  which  depend  on  the  details  of  the  barium  cloud-ionospheric 
interaction.  Although  the  same  elements  enter  into  an  understanding  of 
why  striations  develop,  steepen,  and  bifurcate,  these  processes  are 
complicated  by  the  need  to  consider  the  finite,  albeit  large,  parallel 
conductivity  and  the  small  perpendicular  electron  diffusion. 

We  develop  in  this  section  the  fundamental  fluid  equations  that 
govern  ion  cloud  behavior.  Subject  to  the  assumptions  discussed  below, 
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these  equations  can  be  reduced  to  two  basic  relationships:  1)  An  equation 
stating  the  zero  divergence  of  the  total  current  consisting  principally  of  the 
ion  Pedersen  current  (including  the  contributions  of  gravitational  and 
pressure  gradient  drifts)  perpendicular  to  the  magnetic  field  and  electron 
current  along  the  field;  and  2)  the  electron  continuity  equation  which  shows 
that  the  electron  density  changes  because  of  E  x  B  drifts,  perpendicular 
and  parallel  diffusion,  and  the  divergence  of  electron  current  flow  along  B  . 

We  further  show  that  to  lowest  order  in  (k  k.)  *  ( k  and  k.  are  the 

G  1  G  X 

ratios  of  angular  gyrofrequency  to  collision  frequency  for  electrons  and  ions 
respectively),  the  problem  can  be  reduced  to  a  determination  of  parallel 
drift  (due  to  gravity)  and  ambipolar  diffusion  along  the  field  lines,  plus  a 
description  of  the  change  in  column  density  brought  about  by  incompressible 
E  x  B  convection.  The  effective  electric  field  is  composed  of  the 
amhient  field  (in  the  frair  e  of  the  neutrals),  Eq,  a  contribution  due  to  gravity, 
r  g  x  B  ,  and  a  polarization  field,  E^,  generated  by  the  presence  of  the  cloud, 
and  r  is  the  barium  ion-neutral  collision  time. 

a)  Assumptions 

We  consider  an  isothermal,  inhomogeneous  plasma  cloud  in  a 
finite,  uniform  ionosphere.  The  ionosphere  is  terminated  by  insulating 
boundaries  that  model  the  effect  on  : ^e  current  flow  of  the  unionized  lower 
atmospheie  and  the  collisionless  topside  protonosphere.  A  uniform  rnagnel  '.c 
field  threads  the  cloud  and  its  surroundings  and  is  perpendicular  to  the 
boundaries.  Gravity,  which  acts  only  on  the  ions,  has  both  a  perpendicular 
and  parallel  component  as  indicated  in  Fig.  3. 

We  assume  that  the  ambient  electric  field  and  the  neutral  wind  are 
constant  and  uniform.  All  changes  in  ion  cloud  dynamics  due  to  image 


motions  in  the  ionosphere  subtended  by  the  cloud  and  the  generation  of  these 
ionospheric  perturbations  will  be  neglected.  Large  clouds  are  probably 
not  affected  by  the  secondary  fields  generated  in  this  fashion  although  the 
lower  ionosphere  may  be  greatly  modified  on  the  field  line  below  the  cloud. 

Our  concern  will  be  with  low  p,  weakly  ionized  plasmas  that  change 
slowly  compared  to  particle  collision  times.  Inertial  terms  in  the  equations 
of  motion  will  be  small  under  these  circumstances  and  can  be  neglected 
except  at  very  high  altitudes.  Generally,  finite  Larmor  radius  effects  can 

also  be  neglected  although  they  too  may  be  important  at  very  high  altitudes. 

3  4 

Both  inertial  effects  and  the  effective  "kinematic  viscosity"  associated 

with  the  complete  ion  pressure  tensor  modify  the  dispersion  relation  for 

the  gradient-drift  instability  and  set  a  threshold  for  growth  in  the  limit  of 

collisionless  ions. 

At  altitudes  where  large  ion  yields  are  attainable,  collisions  between 
electrons  and  neutrals  or  ions  are  extremely  infrequent  on  the  time  scale 
of  an  electron  gyro-  eriod.  Ion  collisions  are  more  numerous  on  a  correspond¬ 
ing  time  scale  although  they,  too,  are  infrequent.  As  a  result  parallel 
diffusion  coefficie:  ts  and  conductivities  far  exceed  perpendicular  transport 
properties  and  Pedersen  cum  ents  always  dominate  Hall  currents.  If  we 
denote  the  ratio  of  an  angular  cyclotron  frequency  to  a  collision  frequency 
by  K  =  QT  ,  then  for  case  ;  of  inter  ^st 


K  „  >  K  •  »  K.  >  1  . 
en  —  ei  i 


Although  these  quantities  are  functions  of  altitude,  we  will  in  general 
neglect  these  variations,  taking  the  neutral  atmosphere  to  be  of  constant 
density.  A  relaxation  of  this  assumption  will  be  made  in  the  discussion  of 
ambipolar  diffusion  along  the  magnetic  field. 


b)  Three-Dimensional  Equations  of  Motion 


We  choose  a  rectangular  coordinate  system  in  the  frame  of  the 

neutrals  with  the  z-axis  pointed  along  the  magnetic  field  IT  and  the  y-axis 

aligned  with  Eq,  the  ambient  electric  field, as  shown  in  Fig.  3.  This 

system  translates,  as  viewed  in  the  earth's  frame,  with  the  velocity  of 

the  neutral  wind,  V  .  In  this  frame  E  is  the  sum  of  the  ambient  ionospheric 

no 

field  E&  (in  the  earth's  frame)  and  Vn  x  B  associated  with  the  moving 
reference  system, 


E  =  E  +  V  x  B, 
o  a  n 


Electron  and  ion  drag  are  balanced  in  the  quasistatic  moment 

equations  by  the  conservative  forces  F  .  and  the  magnetic  component 

6  ,  1 

of  the  Lorentz  force: 


F  -Vx/b=K‘1V+/c.(V-  V.).  (1) 

e  e  en  e  ex  '  e  v  '  ' 
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The  "forces"  are  given  by 
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where  b  =  B/B  is  a  unit  vector  in  the  direction  of  the  magnetic  field  and  E 

is  the  sum  of  the  ambient  and  cloud  generated  fields,  E  +  E"  .  The  common 

op 

electron  and  ion  temperature  is  expressed  in  energy  units. 

In  addition,  we  require  quasineutrality 


n.  ~  n  = 
l  e 


n 


(5) 


and  solenoidal  current  flow 


V-  T=  e  V-  (f.  -  T  )  =  0  (6) 

1  e 

whore  r  and  are  the  ion  and  electron  fluxes,  respectively.  Two 
other  equations  are  needed,  electron  continuity 
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and 

V  x  E  =  0.  (8) 

Equation  (8)  implies  the  existence  of  a  scalar  potential  <f>  such  that 


E  -  -V  (f>  (9) 

Note  thal  $  includes  an  ambient  contribution  -Eq  •  x  .  Before  proceeding 
further  it  is  useful  to  separate  §  into  components  due  to  ambipolar  parallel 
diffusion  and  gravitational  fall  and  a  potential  <j>  that  drives  currents 
along  and  normal  to  B  .  Thus 


<b  -  <j)  + 


T(l-  A)  n 
e  (1  f  A  )  n 
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where  A  ~  K{/ Ken<<  *  an<^  na  *s  ambient  ionospheric  density,  Zabusky 

5 

et  al  have  derived  the  expressions  for  the  parallel  and  perpendicular  fluxes 
in  terms;  of  <b  and  Vn/n.  By  substituting  into  Eq.  (6)  and  then  deleting 
small  terms  on  the  order  of  k ^  *,  A,  and  K^/Kq,  they  have  derived  the 
following  simplified  equation  for  the  potential  <j>  : 
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This  equation  states  that  the  divergence  of  the  ion  Pedersen  current  and 
transverse  gravitational  and  pressure  gradient  drifts  must  be  balanced 
by  electron  current  flow  along  B  . 

The  electron  continuity  equation  (7)  can  also  be  simp’ified  by  eliminating 
the  V|  •  nV  d>  term  by  using  Eq.  (11).  The  result  is 


Vi"'  xb  +  8,1  V  n  -  7b  Vm  Ki  vn  " 


2T  - 1  1  ~ 

7B  VlKei  Vin  +  T  Vnn*e  VM»  =  °- 


This  equation  states  that  changes  in  the  electron  density  ire  due  to  E  x  B 
drifts,  gravitational  fall,  parallel  and  perpendicular  diffusion  and  the 
divergence  of  the  electron  current  along  the  magnetic  field. 

These  equations  can  be  reduced  to  a  particularly  simple  and  useful 
form  by  neglecting  diffusion  (low  temperature  limit)  and  gravity  and  assuming 
that  the  potential  is  constant  along  field  lines.  With  these  restrictions, 
which  will  be  discussed  in  the  next  section,  Eqs.  (11)  and  (12)  become 


V  ,•  n  V  i<j)  -  0 
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and 


3  n  l  _  „  —  /\ 

-  p-V.n.V.^xb, 
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If  we  introduce  a  new  potential  v|j  defined  by 


(14) 


$  =  i \)  -  Eq  .  x 


and  integrate  Eq.  (13)  along  B,  then 


V  .  (NV.  4>)  =  E  .  V  .  N 
11  o  1 


(15) 


This  Poisson  type  equation  determines  the  polarization  fields  in  terms  of 

r  l 

the  ratio  N/Nq  where  N  =  J  ^  n  dz  and  Nq  is  the  ambient  ionospheric 
column  density.  The  integral  of  Eq.  (14)  yields 


3  N 


/Tn 


—  +  (Uo  -  V  +  x  b/B)  .  VN  =  0 


(16) 


...  _♦  .  2 
which  describes  an  incompressible  plasma  motion  with  the  local  E  x  B/B 

velocity  and  Uq  is  simply  the  ambient  Eq  x  B/B^  drift.  Alternatively, 

we  could  have  written  Eqs.  (15)  and  (16)  as 


V  •  NE  =  0 


(1?) 


3  N  ,  E  x  B  __T  „ 

+  - o—  .  V  N  =  0. 


3  t 


B 


/  >  r,  * 

Vi  o; 


a  form  used  extensively  in  later  sections. 

Parallel  diffusion  can  be  introduced  by  assuming  a  separable  density 
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function  and  noting  that 

n  (x,  y,  z)  =  f  (z)  N  (x,  y). 


Then  one  finds  from  Eq.  (12)  that  f  must  satisfy  the  generalized  one¬ 
dimensional  ambipolar  diffusion  equation 


3  f 
3  t 


+ 


2T  _3_  _3_ 

eB  3  z  K  i  3  z 


0. 


(19) 


c)  Simplifying  the  Model 

5 

An  important  simplification  is  achieved  by  recognizing  that  the 
large  parallel  conductivity  implies  nearly  constant  potentials  along  the 
field  lines.  We  assume  that  cj>  can  be  expanded  in  powers  of  (/cg /c)  S 
thus 


where  A  =  KeKi  *s  treated  as  a  constant.  From  the  potential  Eq.  (11) 
in  lowest  order,  we  find 


3 

3  z 


n 


3  z 


0  . 


(20) 


Since  the  insulating  boundary  conditions  require  I  =  0,  Eq.  (20) 

E.  -  L 

_ I Q\  _ 

implies  Vm  <t>  '  1  =  0.  Thus  the  first  order  potential  equation  becomes 


Vi’  (nVl>  (0>)  "  Vln'  !TgxB>+7  V±2n+ 


(4.1 


=  0. 


(21) 
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If  we  integrate  Eq.  (21)  along  B  and  note  that  V,|  <j)  ^  vanishes  at  the 
boundaries  then  we  obtain 


Vx  •  (N 


«V 


—  —  2T  2 

VN-rgxB  +  —  V  N  =  0  (22) 

6  _L 


where  the  electron  column  density  N  =  /  n  dz  is  proportional  to  the 

•  J'L  V 

height-integrated  Pedersen  conductivity  p  since  and  K g  are 

assumed  to  be  constant. 

We  introduce  a  polarization  potential  ^  through 


T(°)  .  2T  .  N  -  - 

<4  -  ib  -  -  In  -r——  -  E  •  x  . 

v  ^  e  NQ  o  i 


(23) 


where  Nq  is  the  ambient  electron  column  density  which  is  proportional  to 

Si 

L,  the  ambient  integrated  Pedersen  conductivity.  Combining  Eqs.  (22) 

p 

and  (23)  we  obtain  the  basic  two-dimensional  potential  equation 


Vx  •  (N  i |i)  =  E  *  .  Vjl  N  (24) 

where  _  _  _ 

E  =E  +  V  xB+TgxB 
an  ° 

is  an  effective  background  field - the  source  of  transverse  Pedersen 

currents.  The  potential  ijj  extends  beyond  the  cloud  but  becomes  vanishingly 
small  at  large  distances.  Note  that  the  perpendicular  component  of  gravity 
can  also  drive  currents.  In  fact  it  acts  like  an  electric  field  in  the  potential 
equation  of  magnitude  r  gB  which  for  barium  ions  is  about  0.  3  mV/m 
at  midlatitudes  and  altitudes  of  185  km. 
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We  can  proceed  the  same  way  with  the  electron  continuity  equation, 
L 

and  n  vanish  at  the  boundaries,  to  obtain  the  second 
-L 


noting  that  V||  n 
basic  relationship 


8  N  (E  -  V  i}i)  x  b 

-  +  — 2 -  •  VN  =  0 


8  t 


(25) 


B 


or  re-expressed  in  terms  of  E  =  Eq  - V 4*, 


8N,  ETxET  ^ 

-  +  - —  •  V  N  =  0 


(26) 


8 1 


B 


Since 

dN  8  N  —  M 

the  preceding  equations  imply  that  the  column  density  is  incompressible 
and  hence  that  deformations  are  caused  by  the  relative  E  x  B  convection 
of  various  parts  of  the  cloud. 

To  summarize,  we  find  that  ion  cloud  drift  and  deformation  can  be 
described  by  parallel  diffusion  and  incompr es sib'  onvection.  It  can  be 

shown  that  if  we  assume  that  the  plasma  density  is  separable,  then  <5  ^ 
is  also  vanishingly  small.  Rigorously,  our  model  does  not  lead  to  a 
separable  z-dependence.  However  if  the  Pedersen  conductivity  of 
the  cloud  dominates  the  background  then  the  z-dependence  of  n  will  be 
nearly  cepuruolc  wxtnxn  the  clcuc  unez  tm-  auiutctit-  piao  in  a  wiu  pi'OuuLc 
the  parallel-resistivity  effects.  Thus,  the  density  variation  along  the  mag¬ 
netic  field  follows  from  the  one-dimensional  parallel  diffusion  equation  (19) 
with  gravitational  fall  which  has  been  studied  by  Bannister  and  Davis^  and 

7 

applied  to  barium  plasmr  Touds  by  Pendya’a  as  discussed  in  Chapter  1 
of  this  report. 
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An  extensively  studied  feature  of  barium  ion  cloud  behavior  is  the 
ambipolar  diffusion  of  the  ions  parallel  to  the  magnetic  field.  The  parallel 
diffusion  coefficient  (cf.  Eq.  (19))  is  given  by 


so  that  its  measurement  yields  an  estimate  for  the  barium  ion- nitrogen 
collision  time,  T  .  Measurements  have  yielded  collision  times  in  agree¬ 
ment  with  earlier  classical  measurements  of  barium  ion  mobility  in  nitrogen. 
At  altitudes  around  185  km,  the  effective  collision  time  is  approximately 
0.7  s  and  we  obtain  a  value  of  25  for  the  product  of  the  ion  gyrofrequency  and 
the  ion-neutral  collision  time,  =  (Qt)^  The  perpendicular  diffusion 
coefficient  based  on  the  measured  ic  .  is  several  orders  of  magnitude  smaller 
than  the  parallel  diffusion  coefficient  and  can  be  taken  co  be  zero  as  far  as 
large  scale  motion  of  ion  clouds  is  concerned.  We  shall  see  later  that  the 
value  of  the  perpendicular  diffusion  coefficient  is  exceedingly  important  in 
determining  the  striation  scale  size  and  lifetime. 

a  number  of  optical  analyses  have  yielded  a  parallel  diffusion  constant 
2  . 

of  0.  1  km  /s  at  altitudes  around  185  km.  Figure  4,  taken  from  Ref.  20, 
shows  the  agreement  between  the  observed  ion  cloud  thickness  along  B 
as  inferred  from  Thomson  scatter  measurements  at  Site  C- 6  and  the  calcu¬ 
lated  value  based  on  a  solution  of  the  one-dimensional  diffusion  cauation.  The 
agreement  is  excellent  for  the  first  10  or  15  minutes  for  both  Event  Spruce 
and  Event  Olive  I.  At  later  times,  the  Thomson  scatter  radar  line  of  sight 
through  the  ion  cloud  was  not  parallel  to  B. 


The  lack  of  diffusion  perpendicular  to  the  magnetic  field  implies  that 

in  a  uniform  atmosphere  the  maximum  electron  density  in  the  ion  cloud 

_  i  h 

shouLd  decrease  as  t  '  in  accord  with  one-dimensional  parallel  diffusion. 
The  general  behavior  of  the  peak  electron  density  as  a  function  of  time  is  in 
broad  agreement  with  this  prediction.  In  an  exponential  atmosphere,  the  de¬ 
crease  in  the  peak  density  is  limited  because  the  cloud  comes  into  diffusive 
equilibrium.  The  bottom  side  of  the  ion  cloud  will  become  very  flat  while 
the  topside  becomes  broader  than  a  gaussian.  In  addition,  gravity,  a  down¬ 
ward  or  northward  component  of  the  neutral  wind,  or  a  westward  component 
of  the  electric  field  will  cause  the  location  of  the  peak  electron  density  to 
descend.  The  change  in  altitude  will  result  in  a  compression  of  the  ion  density 
profile  since  the  neutral- atmosphere  scale-height  is  smaller  at  lower  altitudes. 
These  effects  are  illustrated  in  Fig.  5  taken  from  Ref.  7.  Not  shown  in 
this  figure,  but  discussed  in  detail  in  Chapter  1  of  this  report,  is  the  fact 
that  the  peak  electron  density  can  slowly  increase  as  the  ion  cloud  descends 
in  altitude.  Figure  5  was  drawn  under  the  assumption  that  the  vertical  drift 
due  to  the  neutral  wind  and  electric  field  effects  was  constant  in  altitude. 

Let  us  turn  to  the  question  of  the  value  of  the  diffusion  coefficient 
perpendicular  to  the  magnetic  field.  The  cloud  generated  electric  field, 
which  generally  extends  beyond  the  cloud  boundaries,  determines  the  dom¬ 
inant  diffusion.  It  has  been  shown  to  first  order  in  (K  eK  J  *  and  for  large 
racios  of  cloud  to  background  height- integrated  Pedersen  conductivities 
(see  Eq.  (12))  that  perpendicular  diffusion  proceeds  at  a  very  slow  rate  con- 

4 

trolled  by  electron  collisions.  A  low  cloud  density  expansion  due  to  Simon 
leads  to  the  result  that  the  perpendicular  diffusion  is  controlled  by  ions. 


-129- 


DISTANCE  ALONG  THE  MAGNETIC  FIELD  (KM 


Simon's  solution  is  valid  as  long  as  the  changes  in  electron  density  which 
occur  as  a  result  of  the  higher-order  electric  field  do  not  give  rise  to  a  large 
change  in  the  zero-order  electric  field. 

2.  Two-Dimensional  Model 


We  turn  our  attention  to  the  dynamics  of  barium  ion  clouds  in  the  plane 
transverse  to  the  magnetic  field.  These  dynamics  are  best  studied  by  the 
system  of  coupled  nonlinear  two-dimensional  equations  (15)  and  (16)  derived 
from  the  moment  equations  by  integratin  g  parallel  to  the  magnetic  field.  These 
equations  can  also  be  obtained  by  starting  with  simplified  equations  of  motion 
for  electrons  and  ions  in  the  limit  of  zero  electron  collision  frequency.  The 
resulting  two-dimensional  equations  express  the  conservation  of  electron 
flow  and  zero  divergence  of  the  total  current.  In  a  purely  two-dimensional 
model  everything  is  uniform  parallel  to  the  magnetic,  field  and  there  is  no 
flov  of  either  species  in  the  vertical  direction.  However,  the  same  two- 
dimensional  equations  also  represent  the  limit  of  vanishingly  small  parallel 
resistivity  in  which  the  electrons  are  allowed  to  flow  freely  up  and  down 
the  magnetic  field  lines  but  ions  are  restrained  to  flow  perpendicular  to  the 
magnetic  field  in  response  to  the  transverse  electric  fields. 

Figure  6  shows  a  schematic  of  the  cross  section  of  the  ion  cloud  in 

the  plane  perpendicular  to  the  magnetic  field.  The  ion  cloud  is  moving  to  the 

right  with  respect  to  the  neutrals.  In  the  two-dimensional  model,  the  ratio 

of  height-integrated  conductivities;  Z~/Zj*  ,  can  be  represented  as  the 

ratio  of  the  cloud  center  dens.ty,  Nc>  to  the  ambient  ionospheric  density, 

N  .  Due  to  the  slow  rate  of  development  of  the  ion  cloud  and  the 
o 
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small  plasma  pressure  relative  to  the  magnetic  field  pressure  (low  |3),  the 
ambient  electric  fields  are  electrostatic  as  expressed  by  the  condition 
V  x  E  =  0,  Since  the  electrons  do  not  experience  any  collisions  they  flow 
at  the  local  E  x  B/B2  velocity.  The  flow  of  the  ionization  is  incompressible 
as  shown  in  Section  C.  1.  The  fact  that  the  flow  is  incompressible  is  ex¬ 
tremely  important  in  theoretically  determining  the  behavior  of  baiium  ion 
clouds.  This  assumption  precludes  the  possibility  of  an  electron  density 
increase.  The  ion  motion  closely  follows  the  electron  motion  but  the  ions 
do  have  a  component  of  velocity  in  the  direction  of  the  electric  field  as  a 
result  of  collisions  with  the  neutrals.  At  high  altitudes  the  ion  velocity 
with  respect  to  the  neutrals  is  given  approximately  by 


E  x  B  ,  J_  _E 
B2  B  ' 


(27) 


The  simple  set  of  coupled  nonlinear  differential  equations  that  des¬ 
cribe  the  motion  of  barium  ion  clouds  transverse  to  the  magnetic  field  are 
repeated  here 

+  -E'XZ—  •  V  N  =  0;  (28) 

B 

V  '  7  a  V  •  NE=  0;  (29) 

P 

V  x  E  =  0.  (30) 
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Equation  (28)  is  the  continuity  equation  for  the  incompressible  flow 
of  the  electrons  at  the  local  E  x  B  /B2  velocity.  Equation  (29)  describes 
the  conservation  of  the  current  flow  perpendicular  to  the  magnetic  field  in 
the  ionosphere  and  Eq.  (30)  is  the  quasistatic  assumption  of  the  low- (3  approxi¬ 
mation.  The  Pedersen  current,  ,  is  proportional  to  the  product  of  the 
conductivity,  N  (proportional  to  the  ion  density)  and  the  electric  field,  E  , 
all  measured  in  the  frame  of  the  ambient  neutrals. 


Three  (Afferent  approaches  have  been  used  successfully  to  obtain 
information  about  the  nature  of  the  solutions.  A  number  of  analytical 
studies  have  been  carried  out2,<*’  **  and  numerical  calculations  have 

been  used.  ^  A  third  method  involves  the  use  of  a  hydrodynamic  analog*  * 
discussed  below.  This  experiment  has  been  particularly  instructive  in  re¬ 
gard  to  the  detailed  structure  of  striation  development, 
b)  The  Hele-Shaw  Cell 


A  variable-density  but  incompressible  fluid  confined  between  two 
closely  spaced  plates  and  moving  under  the  influence  of  gravity,  Fig.  7 
obeys  the  following  momentum  equation  in  the  quasi- steady  limit: 


-  J  Vpl  qV2  U  +7  =  0 

where  p  is  the  pressure,  p  the  fluid  density,  q  the  kinematic  viscosity  and 
u~  the  fluid  velocity.  For  small  plate  separation  the  viscous  drag  may  be 
estimated  by 

T)  V2  u  ^  -  8r)  ^ /D2 
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where  is  the  midplane  velocity  and  D  is  the  plate  spacing.  Since  the 
fluid  is  incompressible  (V  •  u  =  0)  there  exists  a  two-dimensional  stream 
function  £  (x,  y)  such  that 

u  =  -  V  x  ((•  ia  ) 

where  b  is  a  unit  vector  perpendicular  to  the  plate  It  follows  from  this 
relation  that 

V,  •  (pv,  5)  =  ia-V  "b 

and 


rv 

-*$■  +  U  •  V  p  =  o 

If  we  neglect  gravity  and  write  Eq.  (24)  in  the  frame  where  the  ambient  E 

field  vanishes  with  E  =  V  x  B.  we  obtain 

n 

V  .  (NV  * )  =  (V,NxV  ).b 

where  is  the  uniform  neutral  wind  and  $  -  ^/B.  This  equation  is  in 

exactly  the  same  form  ?.s  the  fluid  motion  relationship  above.  Thus  the 

1  3 

correspondence  shown  in  Table  1  can  be  drawn, 
c)  Ionospheric  Currents 

When  the  barium  ion  cloud  is  considered  to  be  of  finite  extent 
parallel  to  B,  currents  flowing  in  the  ionosphere  below  the  cloud  are  affected 
by  the  presence  of  the  ion  cloud  and  give  rise  to  currents  flowing  parallel 
to  the  magnetic  field. 


TABLE  1 


DUALITY  BETWEEN  HELE-SHAW  CELL  AND  SIMPLIFIED 
TWO-DIMENSIONAL  MODEL  OF  AN  ION  CLOUD 


^vi  '  >=  VPX 

cr  =  -  v  x  a 
•  viP  =  0 

CO 

Boundary  condition  at 
u*=  0;  £  =  constant 


bxVi  ■  (NVj_  $  )  =  V  N  X  Vn 
r  =  -  V  X  ( 

If  +  u  .  VXN=  0 

00 

Boundary  condition  at 
u"  =  0;  $  =  constant 


Fluid 


Plasma  f  in  drift  frame  (E  =  0)1 

•  CO  1 


Density  (p  ) 

Fluid  velocitv  (  u  ) 
"Free-fall" 

2  * 

velocity  (  D  g/8r)  ) 
Stream  function  (  £  ) 


Electron  Column  Density 
Plasma  velocity  (u) 

Neutral  wind  velocity  in 

the  orift  frame  (V  ) 

n 

Electric  poten  .al  $/B 


(N) 
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In  Fig.  8  we  indicate  qualitatively  the  effect  that  the  presence  of 
the  highly  conducting  barium  ion  cloud  has  on  the  current  flow  pattern  in 
the  ionosphere.  Figure  8a  is  a  schematic  diagram  in  the  plane  perpen¬ 
dicular  to  the  magnetic  field  showing  the  distortion  of  the  current  flow  due 
to  the  presence  of  the  more  conducting  ’on  cloud.  In  a  completely  two- 
dimensional  model  this  same  current  flow  pattern  would  exist  at  every 

altitude  parallel  to  the  magnetic  field  and  drift  as  a  "current  cell"  with  the 
—  2 

E  x  B/B  velocity.  That  situation  is  shown  schematically  in  Fig.  8b  which 
represents  a  cut  in  the  yz  plane  through  the  center  of  the  ion  cloud.  Note 
that  the  current  density  in  the  ion  cloud  is  larger  than  in  the  less  dense 
ionosphere.  From  Fig.  8a  we  see  that  this  increase  results  from  the 
concentration  of  the  current  in  the  ion  cloud. 

In  Fig.  8c  the  horizontal  currents  are  the  ion  Pedersen  currents 
whereas  the  parallel  currents  are  due  to  electrons  flowing  along  the 
magnetic  field.  Figure  8d  represents  a  generalization  to  a  three- layered 
ionospheric  model.  The  more  conducting  layer  at  the  top  of  this  figure 
represents  a  conjugate  ionosphere  which  may  be  highly  conducting  and 
therefore  car.  contribute  a  large  perpendicular  current  density.  Again  the 
presence  of  a  conducting  ion  cloud  in  the  local  ionosphere  can  effect  the 
current  flow  in  the  conjugate  ionosphere  since  electron  currents  flow  freely 
along  the  magnetic  field  through  the  magnetosphere. 

While  it  is  recognized  that  the  diagrams  given  in  Fig.  8c  and  8d 
are  qualitatively  correct  insofar  as  they  indicate  current  flow  in  different 
layers  of  the  ionosphere,  they  are  strictly  correct  only  if  (1)  there  is  no 
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Changes  in  the  Ionospheric  Current  Flow  Due  to  the  Presence 
of  a  Highly  Conducting  Barium  Ion  Cloud;  a)  Height-Integrated 
Current  Density  Flowing  Perpendicular  to  the  Magnetic  Field; 
This  Same  Pattern  Results  in  the  Cases  of  b)  a  One- Layered 
Ionosphere,  c)  a  Two- Layered  Ionosphere,  and  d)  Local  Plus 
Conjugate  Ionosphere 
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distinction  between  the  barium  ions  and  the  ambient  ionospheric  ions; 

(2)  if  the  neutral  density  is  the  same  in  all  three  layers;  and  (3)  if  the 
electric  field  which  drives  the  currents  is  externally  imposed  on  the  com¬ 
plete  system.  The  actual  interaction  of  several  different  layers  in  the 

14 

ionosphere  with  a  barium  ion  cloud  is  only  now  beginning  to  be  considered 
and  is  a  very  complex  problem.  It  is  clear  that  the  plasma  density  will  be 
increased  where  the  electrons  flow  into  the  lower  ionosphere  and  decreased 
where  electrons  leave  to  neutralize  the  charge  accumulation  in  the  cloud.  Thus 
the  ionosphere  subtended  by  the  cloud  will  be  locally  perturbed  and  this 
perturbation  will  move  with  the  cloud. 

We  can  see  qualitatively  that  the  effect  of  an  ion  cloud  on  the  current 
flow  pattern  depends  sensitively  on  the  cloud  conductivity  compared  to  the 
conductivity  of  the  ionosphere.  The  appropriate  parameter  that  characterizes 
the  cloud-ionosphere  interaction  is  the  ratio  of  the  height- integrated  Pedersen 
conductivity  of  the  barium  ion  cloud  to  the  height- integrated  Pedersen  con¬ 
ductivity  of  the  ambient  ionosphere, 

D.  LARGE  SCALE  BEHAVIOR 
1 .  Motion  of  Ion  Clouds 

Equation  (29)  indicates  that  at  any  instant  in  time  the  electric  field 
depends  on  the  distribution  of  ionization  in  the  ion  cloud.  In  particular  the 
r  olarization  electric  field  due  to  the  presence  of  the  ion  clov.d  reduces  the 
total  field  at  the  center  of  the  cloud  as  indicated  at  the  top  of  Fig.  9, 

As  a  result,  the  center  of  the  ion  cloud  moves  with  a  slower  velocity  with 
respect  to  the  neutrals  than  does  a  background  ion  that  is  flowing  with  the 
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velocity  U  .  The  value  of  the  electric  field  it  the  center  of  the  ion  cloud 
o 

depends  on  the  ratio  of  the  height- integrated  Pedersen  conductivity  of  the 

ion  cloud,  Z^  ,  to  the  height- integrated  Pedersen  conductivity  of  the  back- 

ground  ionosphere,  Z  ,  as  well  as  on  the  distribution  of  ionization  in  the  ion 

P 

cloud.  For  a  uniform  ionosphere  and  a  constant  neutral  density,  Z^/Z^' 
can  also  be  expressed  as  the  ratio  of  height- integrated  cloud  to  ambient 
density,  Nc/Nq.  The  electric  field  at  the  center  of  the  ion  cloud,  and  thus 
its  motion,  can  be  expressed  in  terms  of  a  coupling  parameter,  %  ,  defined 
as  the  ratio  of  the  electric  field  at  the  center  of  the  ion  cloud  (in  the  frame 
of  the  ambient  neutrals)  to  the  external  electric  field, 


E  ~  £  r  E  =  £,  E  , 

c  o  pc  o 


(31) 


where  the  subscript  c  represents  the  value  at  the  center  of  the  ion  cloud. 
Since  the  difference  in  velocity  between  the  ion  cloud  and  the  neutral  wind 
is  frame  independent  and  equal  to  E  x  B /B  ,  Eq.  (31)  implies 


Y.  -  V 
l  n 


uc  =  « 


(E  x  B  \ 


(32) 


This  equation  also  defines  the  frame  invariant  velocity  of  the  ion  cloud, 

U  .  It  follows  that  £,  ranges  between  one  and  zero  for  very  small  and  very 


large  ion  clouds  respectively; 


_  E  x  B 

v - a - 

V  .  ?  , 

1  B^ 


(N  «N  ); 
'  c  o 


4  -o, 


V.  —  V 
x  n 


(N  »  N  )  . 

c  o 


(33) 
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Very  small  clouds  or  low  density  perturbations  provide  a  direct  measure  of 

the  ambient  electric  field,  =  -  V.  x  B  ,  while  very  large  clouds  are 

completely  tied  to  the  neutral  wind. 

Theoretical  models  have  been  explored  that  give  analytic  expressions 

for  the  dependence  of  £,  on  N^/N^.  ^  A  well-known  example  corresponds  to 

a  cloud  of  constant  density  N  with  an  ellipsoidal  cross  section  as  illustrated 

in  Fig.  9.  Two  interesting  limits,  a  circular  cloud  with  a  =  b  and  an 

elongated  cloud  with  b«a  are  also  shown.  The  former  is  the  shape  used  by 

15 

Haerendal,  Lust,  and  Reiger  in  deriving  their  now  classic  equation 

for  determining  the  ambient  electric  field  from  the  motion  of  barium  ion 

clouds.  In  the  high  altitude  limit  (k^»  1),  their  equation  can  be  brought  into 

the  form  of  Eq.  (32)  if  the  identification  £  =  2/(1  +  X*)  is  made  where 

X"'  =  (N  +  NJ/N  .  The  elongated  shape  is  most  tied  to  the  neutral  wind 

and  its  motion  corresponds  to  that  of  a  field-aligned  sheet  of  ionization. 

The  degret  of  coupling  also  depends  on  the  density  profile  of  the  ionization 

for  large  density  clouds.  A  cylindrically  symmetric  cloud  with  a  density 

2-4/3 

profile  N  «  (1  +  r  )  '  is  less  coupled  to  the  neutral  wind  for  the  same 

value  of  peak  column  concentration  as  indicated  on  the  right-hand  side  of 
Fig.  9. 

Figure  10  is  a  plot  of  the  coupling  parameter  as  a  function  of  the 
ratio  of  cloud  conductivity  to  ionospheric  conductivity  for  several  of  the 
Secede  releases.  The  estimated  uncertainties  in  the  values  for  £,  are  indicated 
by  the  vertical  error  bars,  and  the  horizontal  error  bars  indicate  a  factor 
of  two  uncertainty  in  the  estimates  for  the  ratio  of  integrated  conductivities. 

The  two  curves  illustrate  the  two  extremes  in  coupling*'*  that  depend  on  the 
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Fig.  10  Coupling  Neglecting  the  Influence  of  the  Conjugate  Ionosphere. 

The  open  circle  is  event  Quince  while  the  half-filled  circles 
are  events  Apple  and  Dogwood.  In  order  from  largest  coupling 
parameter  to  smallest,  the  solid  dots  are  the  Plum  Puff, 
Redwood,  Plum,  Spruce,  and  Nutmeg  clouds  respectively. 
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configuration  of  the  ionization.  The  sheet- like  shape  is  strongly  coupled  to 
the  neutral  wind,  whereas  the  least  coupled  is  the  cylindrically  symmetric 
cloud  with  a  monotonic,  radial  density  profile.  A  first  glance  at  the  data 
would  suggest  that  the  simple  theory  is  inadequate  to  describe  the  degree 
of  coupling  of  ion  clouds  to  the  neutral  wind.  We  see  that  only  the  Quince, 
Apple,  and  Dogwood  releases  come  near  to  agreeing  with  the  theory.  Note 
the  large  disagreement  for  the  Secede  II  evening  releases  since  the  scale  on 
the  bottom  is  logarithimic  and  covers  several  decades.  Thus  the  ion  clouds 
are  observed  to  separate  from  the  corresponding  neutral  clouds  much  faster 
than  the  simple  theory  would  allow.  They  behave  as  though  the  background 
ionosphere  were  far  more  conducting  than  the  measurements  of  the  ambient 
ionosphere  indicate. 

The  observation  that  the  winter  evening  (Secede  II)  releases  are  those 
that  are  in  most  serious  disagreement  with  the  theory  led  to  the  suggestion 
that  the  conjugate  ionosphere  in  the  southern  hemisphere  was  contributing 
to  the  background  ionospheric  conductivity.  This  suggestion  seemed  plausible 
because  the  ionosphere  conjugate  to  Eglin  (at  54  S,  103°W)  is  further 
west  and  was  experiencing  summer  an!  hence  was  sunlit.  On  the  other  hand, 
the  ionosphere  conjugate  to  the  Secede  I  morning  releases  (at  49°S,  66°W) 
was  experiencing  winter  and  thus  was  still  in  darkness.  More  quantitatively, 
when  the  sun  was  9°  below  the  horizon  after  sunset  at  Eglin  in  the  last  half 
of  January,  1971,  and  9°  below  the  horizon  before  sunrise  at  Puerto  Rico 
in  the  first  half  of  May.  1969,  the  sun  was  26°  above  the  horizon  and  23° 
below  the  horizon,  respectively,  at  the  corresponding  conjugate  ionospheres. 
Hence  it  is  reasonable  to  assume  that  the  height- integrated  Pedersen  con¬ 
ductivity  of  the  conjugate  ionosphere  for  Secede  II  had  a  typical  daytime  value 


of  5.  0  mho.  By  adding  these  conductivities  to  the  local  values,  we  obtain 
new  estimates  for  the  ratio  of  cloud  to  ionospheric  conductivity. 

Figure  1 1  is  a  replotting  of  the  data  assuming  that  the  conjugate 
ionosphere  does  contribute  to  the  background  ionospheric  conductivity.  By 
making  this  assumption  we  see  that  there  is  no  longer  any  disagreement  be¬ 
tween  the  simple  theory  and  the  data.  It  is  clear  that  these  clouds  do  not 
move  as  though  they  were  shaped  like  sheets.  However,  the  data  do  not 

clearly  distinguish  between  the  two  other  cylindrical  models  whose  theoretical 

1 2 

curves  are  shown.  Zabusky,  Perkins,  and  Doles  have  solved  Eqs.  (2  8) 
through  (30)  by  finite  difference  techniques  for  a  cloud  with  a  cylindrically 
symmetric  initial  density  distribution,  N=Nq  [1  +  4  (1-  r^/4)^]  for 
r<  2.  The  velocity  of  the  peak  density  in  the  E  x  B  frame  is  0.  63  Vn  . 

The  corresponding  value  for  the  coupling  parameter,  £,  ,  is  0.37.  This  value  1 
for  Nc/N0  =  4  is  shown  as  a  solid  square  in  Fig.  11.  We  note  that 
this  point  is  in  excellent  agreement  with  the  analytical  models,  lying 
closer  to  the  curve  with  a  density  profile  (which  should  provide  an  upper 
limit)  than  to  the  curve  for  a  cylinder  with  constant  density. 

We  note  that  the  agreement  with  the  theoretical  model  includes  data 
in  which  the  ratio  of  conductivities  ranges  over  two  and  one-half  orders  of 
magnitude.  This  agreement  is  highly  suggestive  that  the  conjugate  ionosphere 
must  participate  in  the  dynamics  of  the  motion  of  barium  ion  clouds.  We 
estimate  that  an  Alfven  wave  (the  propagation  velocity  of  an  electromagnetic 
disturbance)  can  travel  from  the  northern  to  the  southern  hemisphere  and 
back  in  half  a  minute,  so  there  is  sufficient  time  to  communicate  the 
polarization  electric  fields  between  the  two  hemispheres.  The  most  sensitive 
test  of  the  theory  would  be  to  make  an  Apple  size  release  (12kg)  at  dawn  and 
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COUPLING  INCLUDING  CONJUGATE  IONOSPHERE 


•'ig.  1  I  Coupling  Including  the  Influence  of  the  Conjugate  Ionosphere. 

The  dots  represent  the  same  clouds  shown  in  the  preceding 
figure.  The  solid  square  represents  the  point  obtained  from 
numerical  modeling. 

-147- 


uiaMtA 


dusk.  The  degree  of  coupling  with  the  neutral  wind  should  change  from  about 
65%  (4  =  0.  35)  to  about  30%  (4  =  0.  7)  for  the  assumed  change  of  a  factor  of 
five  in  the  ratio  of  conductivities. 

There  are  two  other  possible  explanations  for  the  lack  of  agreement 
between  data  and  theory  shown  in  Fig.  10;  either  the  data  are  plotted 
incorrectly,  or  the  theoretical  treatment  is  wrong.  The  largest  uncertainty 
in  the  velocities  listed  is  the  ambient  ionospheric  velocity,  E&  x  B/B^,  which 
has  been  assumed  to  be  less  than  20  m/s  where  better  estimates  were  rot 
available.  The  electric  field  probe  indicates  that  this  velocity  may  be  us 
much  as  60  _+  20  m/s  westward  during  Spruce.  Granting  a  velocity  of  this 
magnitude,  the  effect  on  the  data  points  shown  in  Fig.  10  would  be  to 
reduce  4  hy  a  factor  of  two.  It  is  easy  to  see  by  inspection  that  the  data 
points  would  still  be  in  disagreement  with  theory.  It  is  unlikely  that  this 
large  a  discrepancy  in  E  x  B/B^  will  allow  a  reconciliation  with  theory. 

The  remaining  possibility  is  that  the  theory  does  not  correctly 
describe  the  coupling  of  the  barium  cloud  to  the  rest  of  the  ionosphere.  It 
is  true  that  the  complications  of  this  interaction  are  hidden  by  ascribing 
to  the  ionosphere  a  height- integrated  conductivity,  a  procedure  that  would 
be  strictly  valid  only  if  the  value  of  (S2r)^  outside  the  cloud  were  constant 
and  the  same  as  for  barium.  There  is  no  assessment  yet  as  to  the  effect 
to  be  expected  from  a  more  rigorous  treatment  of  the  dynamics  of  the  rest 
of  the  ionosphere.  On  the  other  hand,  the  assumption  that  the  reconciliation 
between  data  and  theory  is  achieved  by  including  the  effect  of  the  conjugate 
ionosphere  is  a  simple  and  plausible  one. 

There  are  several  important  implications  that  follow  from  assuming 
that  the  conjugate  ionosphere  plays  a  role  in  the  ion  cloud  dynamics.  First 
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is  that  the  cloud-produced  polarization  electric  fields  must  project  into  the 
southern  hemisphere  which  suggests  that  the  assumption  that  the  earth's 
magnetic  field  lines  are  equipotentials  for  a  large  distance  is  reasonable. 

A  second  conclusion  is  that  there  must  be  magnetic- field- aligned  currents 
in  the  magnetosphere,  and  a  third  conclusion  is  that  the  presence  of  the  ion 
cloud  and  its  resulting  motion  in  the  northern  hemisphere  must  cause  a 
disturbance  in  the  ion  density  in  the  conjugate  ionosphere  in  the  southern 
hemisphere.  An  independent  verification  of  any  of  these  effects  would  be 
very  valuable. 

Let  us  adopt  a  simplified  model  of  an  ion  cloud  in  the  ionospheric- 
magnetospheric  system  as  shown  in  Fig.  12  in  order  to  make  a  quantitative 
estimate  of  some  of  these  effects.  We  find  that  the  magnetic  field  lines  are 
highly  conducting  so  that  the  electric  fields  do  not  suffer  any  appreciable 
resistive  losses  in  being  projected  into  the  conjugate  ionosphere.  The 
most  important  effect  is  the  disturbance  in  electron  density  that  the  current 
flow  in  the  background  ionosphere  is  likely  to  make.  The  parallel  currents 
must  be  closed  by  Pedersen  currents  which  are  due  to  the  flow  of  the  back¬ 
ground  ions.  This  ion  flow  will  lead  to  significant  changes  in  the  local 

14 

ionospheric  density.  These  changes  should  be  observable  by  ionosondes. 
In  particular,  for  the  Secede  II  releases,  the  motion  of  the  ion  clouds  should 
have  produced  large  changes  in  electron  density  in  the  conjugate  ionosphere 
in  the  southern  hemisphere.  These  changes  should  have  been  detectable 
had  ionosondes  been  deployed  there, 

2.  Deformations  of  Large  Ion  Clouds 

The  electric  field  varies  over  the  distribution  of  density  in  the  ion 
c'oud  as  discussed  above  and  as  shown  schematically  in  Fig.  9.  Since 
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the  ionization  moves  with  the  local  E  x  B/B^  velocity,  different  parts  of  the 
ion  cloud  move  with  different  velocities.  Thus,  as  the  ion  cloud  moves,  its 
shape  distorts.  Linson^  has  analytically  investigated  the  changes  in  shape 
of  equal  density  cortours  by  expanding  the  solution  of  the  nonlinear  equa¬ 
tions  for  the  equipotentials  (streamlines)  as  a  power  series  in  time.  If  only 
the  first,  time- independent  term  is  retained,  then  the  deformation  of  an 
initially  cylindrically  symmetric  concentration  is  shown  in  Fig.  13a. 
However,  the  streamlines  are  not  time-independent.  By  retaining  the  next 
time-dependent  term,  there  is  a  qualitative  difference  in  the  distortion  of 
the  equal  density  contours  as  shown  in  Fig.  13b.  Note  that  the  contours 
are  convex  outward  on  the  righthand  or  front  side  of  the  barium  ion  cloud 
when  the  more  accurate  expression  is  used.  We  also  note  that  the  backside 
has  become  steeper  during  the  same  time  interval. 

Equations  (29)  and  (30)  have  also  been  solved  by  a  finite  difference 

12 

technique  on  a  digital  computer1  for  the  case  of  an  initially  cylindrically 
symmetric  ion  cloud  with  a  maximum  ratio  of  integrated -conductivity 
Nc/Nq  =  4  as  shown  in  Fig.  14.  The  steepening  of  the  backside--that  is, 
the  side  on  the  left  where  the  neutral  cloud  would  be  --is  prominent.  The 
steepening  shown  in  this  and  the  proceeding  figure  provides  an  explanation 
of  the  fact  that  the  trailing  edge  of  ion  clouds  develop  a  large  density  gradient 
when  seen  from  a  variety  of  views.  We  also  note  that  the  ion  cloud  elongates 
in  the  direction  of  its  motion.  Since  the  flow  as  described  by  Eqs.  (29)  and 
(30)  is  incompressiole,  the  area  enclosed  in  each  equal  density  contour 
remains  constant. 

Analytical  models  yield  an  estimate  of  a  steepening  time.  However, 
such  estimates  are  sensitively  dependent  upon  both  the  assumed  density 


COMPUTER  SOLUTION  SHOWING  STEEPENING 
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profile  in  the  ion  cloud  and  the  definition  of  steepening.  Figure  15  shows 

steepening  times  normalized  by  the  ion  cloud  radius  at  half-maximum  density 

and  the  velocity  UQ  as  obtained  from  two  different  models.  The  curve  (a) 

results  from  a  model^  in  which  the  ion  cloud  density  has  a  slowly  decreasing 
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profile  varying  like  N  «:  (1  4  r  )  '  .  The  steepening  time  in  this  case  is 

defined  in  terms  of  the  time  for  the  peak  density  to  move  a  typical  distance 

over  which  the  streamlines  vary  significantly.  On  the  other  hand,  the  curve 

(b)  represents  the  relative  motion  of  a  cylindrical  ion  cloud  with  constant 

density  and  a  sheet- like  configuration  with  the  same  density.  We  see  that  for 

ion  clouds  with  low  conductivity  the  steepening  time  varies  inversely  with 

cloud  conductivity.  We  also  see  that  for  ion  clouds  with  large  conductivity 

the  steepening  is  very  model  dependent. 

An  estimate  of  the  steepening  time  has  ai»o  been  obtained  from  the 

numerical  solution  shown  in  Fig.  14.  The  gradient  in  density  between 

the  third  and  fifth  density  contours  from  the  center  increase  linearly  with 

time  up  to  t  =  6.  0,  apparently  doubling  at  t  =  4.  5.  If  we  take  this  time  as 

the  steepening  time,  and  normalize  it  with  respect  to  the  half-density  radius, 

a  =  0.  8,  we  find  t  =  5.  63  a/U  which  is  shown  as  a  black  square  in 
'  s  '  o 

Fig.  15,  in  reasonable  agreement  with  the  analytical  models.  An  analytic 
calculation  corresponding  to  Event  Holly^  (  a  2.4  kg  release  at  194  km 
altitude)  is  given  in  Ref.  10  and  is  shown  in  Fig.  15  as  a  black  dot. 

5 

Zabusky,  Perkins,  and  Doles  have  presented  a  simple  elongated  model 
of  an  ion  cloud  that  is  useful  not  only  for  showing  steepening,  but  also 
for  illustrating,  the  narrowing  of  the  middle  portion  of  an  ion  cloud 
after  it  has  become  elongated.  Assume  a  cloud  elongated  in  the 
direction  of  its  motion  (the  x-direction)  such  that  3N/3y  »  3N/3x. 


Steepening  Time  According  to  Two  Different  Models;  a)  Broad 
Density  Profile;  b)  Constant  Density  within  a  Cylinder.  The 
filled  square  represents  a  point  obtained  from  numerical 
modeling.  The  filled  circle  represents  the  steepening  time 
for  Kvent  Holly. 


The  y- component  of  the  electric  field  as  determined  by  Eq.  (29)  is  given 
roughly  by 


and  the  equation  V  x  E  =  0  gives 


SN  f 
6x  o' 


Hence,  the  velocity  components  in  the  drift  frame  (which  best  shows  the 
cloud  distortions)  are 


V 

x 


No  9N  Eo 


and  are  illustrated  in  Fig.  16a.  The  isodensity  contours  then  deform 
into  a  "dumbbell"  shape  narrowing  significantly  in  the  middle  as  illustrated 
in  Fig.  16b. 

These  distortion  features  are  evident  in  the  many  optical  observa¬ 
tions  of  barium  ion  clouds.  For  instance, the  Olive  ion  cloud  shown  in 
Fig.  17  at  17  minutes  after  release  has  elongated  approximately  55  km 
in  the  direction  of  its  motion  and  its  midpoint  had  narrowed  to  a  dimension 
of  the  order  of  1  km  transverse  to  the  direction  of  its  motion.  By  the  same 
time,  the  ion  cloud  had  diffused  parallel  to  the  magnetic  field,  resulting  in 
an  extension  of  over  60  km  in  the  direction  of  the  magnetic  field.  This 
configuration  is  best  described  as  c  sheet,  and  it  is  important  to  take  this 
configuration  of  ionization  into  account  when  attempting  to  describe  radar 
effects.  The  regions  of  highest  densities  do  not  occur  in  this  thin  sheet 
region  but  are  confined  towards  the  back  edge  of  the  cloud  which  becomes 
unstable  and  breaks  up  into  a  number  of  sheets  and  rods  as  we  shall  discuss 
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below. 


(b) 


Fig.  16  Narrowing  of  an  Elongated  Ion  Cloud 
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Fig.  17  Elongated  Olive  Cloud  Stretched  into  a  Thin  Sheet  (Technology- 
International  Corporation) 
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Figure  18  is  a  photograph  of  the  Spruce  ion  cloud  from  Site  C-6 
looking  nearly  up  the  magnetic  field  line.  The  features  described  above 
are  evident  in  this  photograph.  The  ion  cloud  is  highly  elongated  in 
the  direction  of  the  neutral  wind.  The  middle  portion  of  the  ionization  is 
very  thin. 

We  conclude  this  section  by  summarizing  the  principal  features  of 

the  large  scale  behavior  of  barium  ion  clouds.  The  ion  cloud  lengthens  in 

the  direction  of  the  magnetic  field  in  accordance  with  ambipolar  diffusion 

theory  and  distorts  into  a  thin  flat  sheet  of  ionization  with  the  highest  density 

portion  developing  into  smaller  sheets  of  ionization.  On  the  scale  of  the 

whole  cloud,  there  is  negligible  diffusion  of  ionization  across  the  magnetic 

field  in  the  time  scale  of  l/2  hour.  The  electron  density  decreases  as 
_  1  h 

t  in  accordance  with  one  dimensional  diffusion  in  a  uniform  medium. 

At  late  times  this  decrease  slows  as  the  ion  cloud  achieves  diffusive  equi¬ 
librium  in  the  variable  scale  height  neutral  atmosphere  as  discussed  in 
Chapter  1.  The  dependence  of  the  ion  cloud  motion  on  the  ambient  electric 

field,  E  ,  the  neutral  wind  velocity,  V  ,  and  the  ratio  of  height- integrated 
ci  n 

Pedersen  conductivity  of  the  cloud  to  the  ambient  ionosphere,  Nc/No>  is 
well  understood.  An  important  point  is  that  the  conjugate  ionosphere 
participates  in  the  dynamics  of  barium  ion  clouds. 

E.  STRIATIO  NS 


One  of  the  most  interesting  features  of  barium  ion  clouds  is  their 

tendency  to  break  up  into  magnetic-field-aligned  structures  that  have  the 

appearance  of  striations.  The  process  that  leads  to  the  formation  of  these 

3 

striations  is  the  gradient- drift  instability.  This  structuring  occurs  as  a 
result  of  the  tendency  of  the  peak  densities  in  the  ion  cloud  to  move  more  with  the 
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Fig.  18  The  Spruce  Ion  Cloud  Showing  Elongation  in  the  Direction  of  the 
Neutral  Wind  and  Narrowing  in  the  Middle  (Technology  Internati 
Corporation) 
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neutral  wind  velocity.  Thus  the  peak  electron  densities  that  originally  are 
created  near  the  center  of  barium  ion  clouds  eventually  end  up  in  the  fine 
scale  striations  at  late  times.  This  behavior  has  been  show.'  -vith  the  use 
of  numerical  simulations  .  A  great  deal  has  been  learned  about  the 
detailed  development  of  structure  in  the  ion  cloud  by  using  the  fluid-dynamic 
analog  or  Hele-Shaw  cell.  The  behavior  exhibited  by  this  experiment  is 
believed  to  be  relevant  to  the  high-altitude  barium  cloud  case  not  only  be¬ 
cause  the  dynamics  in  the  appropriate  limits  are  described  by  equivalent 
sets  of  equations  but  also  because  we  are  fortunate  to  have  observations  of 
barium  ion  clouds  as  they  passed  through  the  magnetic  zenith  as  seen  from 
Site  C-6  at  a  time  when  they  were  fully  striated.  These  unique  observations 
of  three  different  barium  clouds  allow  us  to  confirm  that  many  of  the  struc¬ 
tures  seen  in  the  Hele-Shaw  cell  also  occur  in  barium  ion  clouds  which 
suggests  that  the  simple  set  of  equations  (2  8)  through  (30)  adequately  describe 
the  features  that  develop  in  high  altitude  barium  ion  clouds.  Before  pro¬ 
ceeding  with  some  of  the  quantitative  results  obtained  from  the  theory,  let 
us  first  give  an  overall  description  of  striation  moiphology. 

Figure  19  is  a  schematic  drawing  of  the  several  stages  of  develop¬ 
ment  of  striations  that  form  on  the  backside  of  ion  clouds.  Initially,  a 
broad  spectrum  of  density  perturbations  can  grow.  There  does  tend  to  be 
a  preferred  scale  size  evident  in  photographs  of  barium  ion  clouds.  The 
perturbations  grow  into  magnetic-field- aligned  sheets  of  ionization.  In  the 
third  stage,  these  sheets  undergo  a  series  of  complex  evolutions  that  include 
pinching  of  sheets  to  form  rods  and  bifurcations  of  the  sheets.  At  late  times 
the  ionization  pinches  off  forming  a  series  of  a  large  number  of  rods.  Small 
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Fig.  19  Morphology  of  Nonlinear  Striations  Deduced  from  Photographs, 
Simulations,  and  Theory 
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rods  are  seen  to  slowly  fade  away  apparently  into  very  thin  sheets.  The 
fact  that  small  scale  sizes  in  the  rods  and  sheets  and  large  density  gradients 
exist  over  times  of  minutes  indicates  that  the  diffusion  of  ionization  trans¬ 
verse  to  the  magnetic  field  in  these  large  density  ion  clouds  is  controlled  by 
electron  collisions. 

1 .  Onset  Time 

Barium  ion  clouds,  while  elongating,  appear  to  remain  rather 

regularly  shaped  lor  significant  periods  of  time  before  finer  scale  structure 

is  seen.  There  apparently  is  a  delay  between  the  formation  of  the  ion  cloud 

and  the  time  at  which  structure  is  first  seen.  The  ability  to  understand 

and  predict  this  onset  time  of  striations  is  very  important.  The  parameters 

discussed  in  the  previous  section  contain  only  two  independent  units  of  time. 

They  are  the  ion-neutral  collision  time,  T  ,  and  the  time  proportional  to  the 

reciprocal  of  the  background  electric  field  measured  in  the  frame  of  the 

neutrals,  U  \  Effects  associated  with  the  first  time  scale  lead  to  diffusion 
o 

which  does  not  generate  structure.  The  ambient  electric  field,  on  the  other 
hand,  is  the  mechanism  that  drives  the  gradient- drift  instability.  It  is 
reasonable  to  expect  that  the  onset  time  of  striations  should  vary  inversely 
as  U  ,  all  other  conditions  remaining  constant. 

In  the  Secede  II  test  series,  both  the  1  kg  clouds  and  the  48  kg 
Spruce  and  Plum  clouds  were  released  approximately  at  the  same  altitude. 
Figure  20  shows  a  plot  of  the  onset  time  of  these  clouds  as  a  function  of 
the  estimated  UQ  *  for  each  of  these  releases.  The  large  horizontal  error 
bars  represent  the  present  uncertainty  in  the  value  of  the  ambient  electric 
field.  We  see  that  there  is  agreement  with  the  expected  behavior  indicated 
by  the  solid  lines.  We  should  point  out  that  the  neutral  wind  can  affect  the 
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peak  electron  density  produced  due  to  the  finite  time  needed  to  ionize 
barium  atoms.  As  a  result,  the  Plum  cloud  did  not  perturb  the  conductivity 
of  the  ambient  ionosphere  to  the  same  extent  that  the  Spruce  cloud  did  and 
hence  they  were  not  identical  clouds.  Figure  20  can  be  used  for  making 
predictions  as  to  the  onset  time  of  striations  for  48kg  releases  at  185  km 
altitude  at  Eglin  AFB  during  the  winter  evening  once  the  values  of  the  neutral 
wind  and  the  ambient  electric  field  have  been  given.  All  of  these  caveats  have 
been  made  because,  although  the  linear  relationship  between  onset  time  and 
Uo  1  can  be  expected,  the  value  of  the  slope  does  depend  on  the  altitude, 
yield,  release  location  on  the  earth,  the  time  of  dav  (morning  dawn  or  evening 
dusk)  and  the  season  of  the  year.  Our  theoretical  understanding  combined 
with  empirical  scaling  of  the  Secede  data  would  allow  us  to  make  fairly  reliable 
estimates  of  what  the  slope  should  be  once  the  future  test  location  and  time 
have  been  specified. 

Figure  21  is  a  purely  empirical  scaling  of  the  onset  time  of 

striations  normalized  to  U  =75  m/s,  as  a  function  of  altitude  between 

o 

150  and  250  km  for  48kg  payloads  released  at  midlatitude  regions.  While 
there  is  no  theoretical  basis  for  the  curve  drawn  in  Fig.  21,  it  should 
provide  the  basis  for  an  accurate  prediction  of  onset  time  once  the 
appropriate  value  of  Uo  =  j  +  Ea  x  B/B  J  is  known. 

Figure  22  is  a  comparison  of  the  normalized  observed  onset  time 
of  striations  as  a  function  of  the  estimated  ratio  of  height- integrated  cloud 
conductivity  to  the  height- integrated  background  ionospheric  conductivity 
with  a  theoretical  curve.  The  scale  to  which  the  onset  times  are  normalized 
is  a/UQ,  the  time  for  the  background  ionization  to  flow  a  distance  equal  to 


NORMALIZED  ONSET  TIME  vs  ALTITUDE 


the  estimated  initial  radius  at  half-maximum  density  of  the  ion  cloud.  The 
error  bars  represent  the  uncertainty  in  the  estimate  of  the  data.  The  pri¬ 


mary  uncertainty  lies  in  the  value  of  the  ambient  electric  field,  E^,  and 

in  the  ratio  of  conductivities,  N  /N  .  The  theoretical  curve  represents  a 

c  o 

steepening  time,  that  is  the  time  it  takes  a  sheet-like  configuration  of 
ionization  to  move  a  distance  equal  to  the  initial  cloud  radius  at  half-maxi¬ 
mum  density  relative  to  the  motion  of  the  center  of  the  ion  cloud.  The 
solid  curve  results  from  a  constant  density,  cylindrical  cloud;  by  using 
the  expressions  given  in  Fig.  9,  we  obtain 


T 


U 


s _ o 

a 


U 

_ o_ 

AU 


+  3  + 


(34) 


c  o 

The  shaded  area  is  an  indication  of  the  variation  that  results  from  using 
different  initial  density  models.  The  similarity  in  the  shape  of  the  theore¬ 
tical  curve  and  the  distribution  of  the  data  points  is  striking.  The  agreement 
in  magnitude  is  fortuitous  because  there  is  no  reason  why  the  numerical 
coefficient  relating  the  onset  time  to  the  dimensional  time,  a/AU,  should 
have  been  one  in  this  nonlinear  problem.  It  appears  that  a  numerical 
coefficient  a  little  greater  than  one  would  give  even  better  agreement.  Note 
that  the  daca  extend  over  two  orders  of  magnitude  in  the  ratio  of  cloud  to 
ionospheric  conductivity. 

The  data  plotted  in  Fig.  22  include  all  of  the  variations  in  yield 
and  altitude.  It  should  further  be  emphasized  that  the  estimates  for  the 
height- integrated  background  conductivity  were  based  on  inclusion  of 
estimated  effect  of  the  conjugate  ionosphere  as  was  done  in  Fig.  11.  If 
only  the  local  ionospheric  conductivity  had  been  included,  all  of  the  Secede  II 


data  points  shown  by  the  solid  circles  would  have  been  displaced  an  order 
of  magnitude  to  the  right.  This  agreement  of  the  observed  onset  times  of 
these  barium  ion  clouds  with  the  theoretical  curves  is  an  independent  indi¬ 
cation  that  the  conjugate  ionosphere  participates  in  the  dynamics  of  barium 
ion  clouds.  We  do  note  that  there  are  two  data  points  that  do  not  agree  with 
the  theoretical  curve.  The  Apple  ion  cloud  striated  later  than  it  should 
have  while  the  Nutmeg  ion  cloud  produced  striations  much  more  rapidly 
than  expected.  At  present  there  is  no  definite  explanation  for  this  discrep¬ 
ancy.  Perhaps  tne  error  bars  on  Apple  should  be  larger;  perhaps  the  lower 
altitude  Nutmeg  cloud  had  large  initial  perturbations. 

In  general,  we  would  conclude  that  we  have  a  solid  basis  for  making 
predictions  of  the  onset  time  of  striations  both  empirically  by  using  the 
data  base  and  by  a  theoretical  model  that  can  be  used  to  obtain  predictions 
under  a  wide  variety  of  conditions. 

2.  Scale  Size  of  Striations 

We  now  turn  our  attention  to  the  question  of  tne  scale  size  of  the 
striations.  Much  of  our  information  has  been  obtained  from  photographs 
that  were  taken  looking  directly  up  the  magnetic  field  line  into  the  highly 
striated  portions  of  the  Redwood,  Spruce,  and  Nutmeg  ion  clouds  as  shown 
in  Figs.  23,  24,  and  25  respectively.  We  note  that  there  is  a  wide 
variety  of  scale  sizes  evident  in  the  sheet  and  rod-like  structures  shown. 

We  shall  examine  some  of  the  details  of  the  fine  scale  structures  below. 

We  do  note  that  in  spite  of  the  variety  of  detailed  structure  present,  there 
is  an  overall  sheet-like  appearance,  particularly  in  the  Redwood  and  Spruce 
ion  clouds.  We  conclude  that  there  is  apparently  a  predominant  overall 
wavelength. 
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Redwood  Striations  Viewed  Directly  Up  the  Magnetic  Field 
Line  (Technology  International  Corporation) 
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Nutmeg  Looking  Directly  Up  the  Magnetic  Field  Line 
(Technology  International  Corporation) 
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We  now  point  out  the  overall  similarity  between  the  photographs  of 

the  barium  ion  cloud  taken  up  the  magnetic  field  line  and  the  hydrodynamic 

17 

analog  experiments  conducted  in  a  Hele-Shaw  cell  by  Thomson.  In  Fig.  26 
we  see  overall  sheet- like  structures  in  addition  to  finer  structures  that  we 
shall  describe  below.  The  fluid  motion  exhibited  in  Fig.  26  is  described 
by  a  set  of  equations  identical  to  Eqs.  (28)  through  (30)  that  apply  to  high 
altitude  barium  ion  clouds  as  discussed  in  Section  C.  The  structures  seen 
in  the  Hele-Shaw  cell  are  identical  in  many  respects  to  those  that  can  be 
seen  in  barium  ion  clouds.  The  Hele-Shaw  cell  has  the  advantage  that  one 
can  follow  the  nonlinear  evolution  of  a  given  sheet  or  striation  without 
becoming  confused  by  its  changing  aspect.  When  viewing  a  barium  cloud 
striation,  it  is  located  near  the  magnetic  zenith  only  for  a  short  time,  and 
as  it  moves  out  of  the  region  close  to  the  zenith  its  appearance  changes  as 
more  of  its  extension  parallel  to  the  magnetic  field  can  be  seen  projected 
into  the  field  of  view. 

In  looking  at  the  detailed  evolution  of  these  sheets,  it  is  seen  that 
typically  a  large  sheet  grows  and  then  bifurcates,  breaking  up  into  two 
smaller  sheets.  Figures  27  arid  28  are  tracings  of  the  characteristic 
shapes  that  evolve  at  the  tip  of  sheets  in  the  Hele-Shaw  experiment.  Figure 
27  shows  a  typical  fiburcation  process  in  both  a  symmetric  and  asymmetric 
form.  Figure  28  shows  a  more  complex  evolution  into  three  structures  on 
the  left-hand  side;  the  right-hand  side  represents  J  .e  growth  of  a  sheet  of 
less  dense  ionization  into  denser  ionization.  It  seems  apparent  that  this 
process  of  a  sheet  growing  and  then  breaking  up  into  smaller  structures 
that  develop  into  larger  sheets  keeps  repeating.  This  process  appears  to 
continue  until  a  sheet  becomes  pinched  off  at  some  location.  The  isolated 
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Fig.  26  Striation  Development  in  the  Hele-Shaw  Cell  (Physical 
Dynamics,  Incorporated) 


-174- 


INITIAL 


Pinch-Off  and  Formation  of  Isolated  Rods  (Physical  Dynam 
Incorporated) 


pie.e  of  ionization  continues  to  split  and  elongate  until  it  reaches  what  is 

apparently  a  stable  size. 

1 8 

Thomson  has  carried  out  a  detailed  study  of  the  evolution  of  the 
tips  of  sheets  as  illustrated  by  the  symmetric  shapes  shown  in  Fig.  27. 

It  was  concluded  that  the  tips  grow  in  size  linearly  with  time.  When  a 
characteristic  size  is  reached,  the  tip  bifurcates  again.  There  seems  to 
be  a  factor  of  three  between  the  maximum  stable  size  and  the  smallest 
aew  size  that  is  formed.  The  time  scale  for  a  tip  to  bifurcate  also  appears 
to  scale  with  the  size  of  the  tip. 

Figure  29  shows  a  computer  simulation  of  this  bifurcation  process. 

In  order  to  get  the  initially  cylindrically  symmetric  rod  or  ion  cloud  to 
bifurcate,  ar>  artificial  density  perturbation  had  to  be  introduced  at  t  =  0. 

We  note  that  the  scale  length  of  the  nonlinear  sheets  far  exceeds  the  scale 
associated  with  the  original  perturbation.  We  also  see  that  the  peak  electron 
density  at  late  times  is  at  the  very  tip  of  the  developing  sheets.  The  sheets 
have  begun  to  pinch  off  at  t  =  18  due  to  the  artificial  diffusion  included  in  the 
computer  routine.  Note  also  that  by  t  =  18  the  background  ionization  between 
the  two  sheets  has  also  bifurcated. 

Figures  30  through  34  show  a  sequence  of  views  of  the  Spruce 
ion  cloud  taken  up  the  magnetic  field  line.  The  magnetic  zenith  is  easily 
identifiable  in  these  photographs  as  the  point  of  convergence  of  the  rayed 
structures.  In  Figs.  30,  32,  and  33  we  see  evidence  of  bifurcated 
sheets.  In  Figs.  30,  32,  and  34  we  see  structures  quite  similar  to  the 
tracings  taken  from  the  Hele-Shaw  cell  shown  in  Fig.  28.  The  similarity 
in  appearance  of  these  phenomena  indicates  that  the  Hele-Shaw  cell  is  a 
useful  simulation  of  the  nonlinear  sheet  de  lopment  process. 
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COMPUTER  SIMULATION  OF  BIFURCATION 


Fig.  29  Deformation  of  Isodensity  Contours  for  an  Initially  Circular 
Plasma  Cloud  with  a  5%  Initial  Perturbation  Showing  a  Simu¬ 
lation  of  Bifurcation  (Bell  Telephone  Laboratories) 
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Fig,  30  Spruce  From  C-6  at  R  +  22:25  (Technology  International  Corp. ) 
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Fig.  32  Spruce  From  C-6  at  R  +  24:05  (Technology  International  Corp.  ) 
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Fig.  34  Spruce  From  C-6  at  R  +  27:08  (Technology  International  Corp. ) 
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Some  of  the  scale  sizes  evident  in  this  series  of  photographs  can  be 
estimated  by  noting  that  the  separation  between  the  pair  of  bright  stars  near 
the  bottom  of  Fig,  31  corresponds  to  approximately  5  km  at  an  altitude  of 
185  km,  the  height  of  the  middle  of  the  rod  along  the  magnetic  field  line. 

In  Fig.  32  note  the  very  thin  sheet  that  is  attached  to  the  newly  formed 
rod  near  the  magnetic  zenith.  This  sheet  appears  to  be  approximately  30  m 
thick  in  the  original  photographs.  Other  thin  sheets  can  be  identified  through 
out  Figs.  31  through  34.  We  will  discuss  the  lifetime  of  these  fine 
structures  below. 

We  call  the  readers  attention  to  a  small-,  medium-,  and  large- size 
rod  in  this  series  of  Spruce  photographs.  A  small  60  m  wide  rod  can  be 
seen  in  the  original  of  Fig.  30  on  the  southernmost  edge  of  the  Spruce  ion 
cloud.  This  rod  is  not  evident  in  Fig.  31  and  appears  to  have  faded  away. 
Note  the  pair  of  bifurcating  sheets  seen  in  Fig.  30.  In  Fig.  31  the 
newly  formed  sheets  have  pinched  off  to  form  rods.  The  smallest  of  these 
appears  to  be  about  120  m  in  width.  It  can  still  be  identified  faintly  in 
Fig.  34,  indicating  a  lifetime  of  at  least  4  minutes.  Directly  below  the 
pair  of  bifurcating  sheets  seen  in  Fig.  30  is  a  larger  rod  t^at  has  a  well- 
developed  sheet  attached  to  it.  This  rod  can  still  be  clearly  seen  in  Fig.  34 
and  appears  to  be  approximately  the  larg  t- sized  rod  that  does  not  bifurcate 

No\ '  let  us  see  what  theoretical  '  nvestigations  have  to  say  about  the 

question  of  scale  size.  The  linearized  stability  analysis  applied  to  the 

3 

backside  of  a  barium  ion  cloud  has  been  extended  in  several  ways,  pri¬ 
marily  by  taking  into  account  the  finite  size  of  the  ion  cloud  and  the  finite 
ratio  of  the  ion  cloud  conductivity  to  the  background  conductivity. 
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The  stability  of  a  slab  with  a  density  profile 
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has  been  analyzed  using  the  two-dimensional  equations  including  electron- 
controlled  diffusion  perpendicular  to  the  magnetic  field.  In  the  respective 
limits  the  growth  rate  assumes  the  form 
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where  k  =  2tt/X  is  the  wavenumber  of  the  perturbivtion  in  the  y-direction 
and  a  =  N  /Nq.  U  represents  the  relative  velocity  between  the  back  edge 
of  an  ion  cloud  and  the  velocity  of  the  neutral  atmosphere. 

A  typical  linear- wave  spectrum  is  shown  in  Fig.  35.  There  is  a 
wavelength  for  maximum  growth,  but  it  should  be  pointed  out  that  the 
wavenumber  scale  is  logarithmic  and  there  is  a  significant  growth  rate 
for  a  wide  range  of  different  scale  sizes.  It  is  reasonable  to  expect  a 
variety  of  wavelengths  attempting  to  grow  initially.  It  should  be  emphasized 
that  >.he  value  of  the  wavelength  for  maximum  growth  depends  sensitively 
on  the  parameter  d  =  |vj£rj  N  j  *  which  so  far  can  only  be  estimated  for  large 
density  barium  ion  clouds  that  have  steepened.  When  striations  have 
formed,  they  exhibit  very  large  gradients  in  density  corresponding  to  small 
values  for  d.  Hence,  the  striations  themselves  should  be  unstable  to  very 
small  wavelengths  as  we  shall  see  below. 
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Table  2  summarizes  a  preliminary  estimate  of  some  of  the  scale 
sizes  identified  in  photographs  of  the  three  events  Redwood,  Spruce,  and 
Nutmeg.  The  dominant  width  tabulated  represents  the  overall  sheet 
size  when  the  ionization  is  fully  striated.  It  is  basically  estimated 
by  dividing  the  distance  between  the  two  outermost  sheets  by  one  less  than 
the  number  of  sheets  that  can  be  counted.  This  distance  is  compared  in  the 
table  with  the  peak  wavelength  given  by  the  linear  wave  spectrum  in  Fig.  35. 
The  appropriate  values  for  U  and  D^were  used  for  each  event.  We  see  that 
the  uniform  assumption  that  d  =  500  m  does  not  give  good  agreement  with 
the  observed  dominant  wavelength.  However,  if  d  =  500  m  is  assumed 
appropriate  for  Spruce  and  d  is  scaled  to  the  initial  transverse  width  of 
each  of  the  other  two  clouds,  then  the  wavelength  for  peak  growth  gives 
excellent  agreement  with  the  observed  dominant  widths.  It  is  not  clear 
why  the  smaller  size  sheets  that  form  larger  sheets  are  never  smaller  than 
100  to  150  m  in  full  width.  If  one  applies  the  two-dimensional  linear  stability 
analysis  to  the  tip  of  a  sheet,  one  calculates  that  perturbations  on  a  wave¬ 
length  scale  of  meters  should  be  unstable.  Hence,  the  fact  that  a  sheet  or 
rod  with  a  transverse  width  of  500  m  does  not  appear  to  bifurcate  suggests 
that  the  two-dimensional  model  may  not  be  appropriate  for  these  small 
scale  lengths. 

5 

Zabusky,  Perkins,  and  Doles  have  pointed  out  that  the  finite  resis¬ 
tivity  parallel  to  the  magnetic  field  in  the  ambient  ionosphere  due  to  electron 
collisions  has  a  significant  effect  on  the  minimum  unstable  wavelength.  The 
finite  resistivity  means  that  magnetic  field  lines  are  no  longer  cquipotentials 
and  electrons  can  cross  field  lines.  Thus  they  no  longer  prevent  ions  from 
crossing  field  lines  and  the  minimum  scale  length  becomes  controlled  by 
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TABLE  2 

COMPARISON  OF  SCALE  SIZE  (X /l  =  2a) 


Event  and 

Scale  Altitude 

Size 

REDWOOD 

2  50  km 

SPRJCE 
185  km 

NUTMEG 

1 50  km 

Maximum  stable  (m) 

500 

500 

? 

Minimum  developed  (m) 

150  •> 

150 

100  ? 

Thinnest  sheets  (m) 

<50 

<6  0 

<50 

Dominant  Width 

Observed  (m) 

300-500 

400 

250 

Theory 

d  =  500  (m) 

170 

400 

500 

d  oc  a  (m) 

340 

400 

300 

ion  diffusion  in  the  ambient  ionosphere.  Although  they  have  not  solved  a 
dispersion  relation,  they  did  obtain  an  estimate  of  the  minimum  unstable 
wavelength  by  using  quadratic  forms; 


/  N  2T  L  1  1 

X  .  =  2?rd  (  -  —  -  - s 

mm  l  ,T  .2  c 

\  N  eEd  d  K  .  K. 
\  o  eA  i 


~ 0 .  8  km.  (37) 


Here  Nc/N  is  the  ratio  of  cloud  to  ambient  electron  content  and  L  is  the 
length  in  the  ionosphere  indicated  in  Fig.  3.  The  above  estimate  is 
arrived  at  by  taking  d  =  1  km,  Nc/Nq  =  2,  T/eEd  -  0.  05,  L/d  -  102, 

4 

-  1*3x10  and  K ^  =  25.  We  note  that  this  estimate  suggests  that 
striations  with  a  width  of  order  400  m  or  less  should  be  stable  in  agreement 
with  the  observations. 


It  seems  apparent  that  sheets  exist  that  are  as  thin  as  can  be  re¬ 
solved.  Since  the  barium  ion  has  a  gyroradius  of  approximately  8  m,  one 
might  conclude  that  16  m  is  the  minimum  size  sheet  that  can  exist.  As  a 
final  comment,  the  density  gradient  lengths  associated  with  rods  can  easily 
be  10  m  or  below.  There  is  evidence  from  the  short  pulse  radar  that 
density  gradients  on  this  short  a  scale  do  exist. 

3.  Lifetime  of  Rods  and  Sheets 

Now  we  turn  our  attention  to  a  process  that  we  suggest  determines 
the  lifetime  of  rods  and  sheets.  Figure  36  shows  a  schematic  of  a  rod 
with  the  ambient  ionospheric  flow  around  it.  A  barium  ion  that  is  trans¬ 
ferred  from  the  rod  into  the  background  ionosphere  will  immediately  acquire 
the  relative  flow  velocity  labeled  U  .  Diffusion  transverse  to  the  magnetic 
field  is  a  process  for  taking  ions  from  the  rod  and  putting  them  into  the 
region  of  the  background  flow.  Rods  are  seen  to  have  very  well  defined 
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Fig.  36  Schematic  of  a  Dissipation  Mechanism  for  Rods 
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edges  indicating  very  sharp  density  gradients.  If  we  assume  that  their 
density  distribution  is  flat  with  sharp  sides  as  indicated  schematically, 
the  transverse  diffusion  across  the  magnetic  field  line  would  broaden  the 
edge  into  a  finite  width,  f.  This  width  is  determined  by  the  value  of  the 
perpendicular  diffusion  coefficient,  Dj_,  and  the  value  of  the  characteristic 
flow  time  of  the  ambient  ionization  around  the  rod,  a/Um,  where  a  is  the 
half-width  of  the  rod.  The  lifetime  of  the  rod  is  then  given  by  the  number 
of  flow  times  it  takes  to  erode  the  whole  rod  by  removii  g  a  layer  of  thick¬ 
ness  C  each  flow  time  as  indicated  by 


T 


3/2 

a 

-n/D.U 
I  ^ 


(3  8) 


The  ionization  that  is  carried  by  the  ambient  flow  ends  up  in  the  sheet 
whose  width,  w,  apparently  can  be  from  below  30  m  up  to  500  m. 

Due  to  the  uniformity  of  the  ionization  in  the  direction  of  the  flow,  the  life¬ 
time  of  the  sheet  is  significantly  larger  than  that  of  the  rod  and  is  controlled 

by  ordinary  one-dimensional  diffusion  transverse  to  the  direction  of  the 

2  - 1  /2 

sheet  --  i.  e.  ,  N~  (w  +  4Djt)  --  which  yields 


l/2 


3w‘ 


D 


(39) 


as  an  estimate  for  the  time  ^^/z  ^or  s^eet  ^ensity  to  decrease  by  a 
factor  of  2. 

Figure  37  shows  a  preliminary  comparison  of  the  lifetime  for  rods 
and  sheets  estimated  by  Eqs.  (3  8)  and  (39)  with  some  of  the  rods  and  sheets 
discussed  in  the  Spruce  photographs  (Figs,  30  through  34).  The  agree¬ 
ment  between  the  data  and  the  theoretical  curves  is  sufficiently  good  to 
encourage  pursuing  this  idea  further. 
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F.  SUMMARY 


1 .  Conclusions 

A  solid  theoretical  basis  has  been  established  for  understanding  a 
number  of  the  complex  phenomena  that  occur  when  barium  ion  clouds  are 
released  into  the  upper  ionosphere.  We  have  indicated  that  a  simple  set 
of  coupled  nonlinear  equations  (see  (28)- (30)  are  capable  of  describing  a 
number  of  these  phenomena.  We  know  what  governs  the  motion  of  barium 
ion  clouds  and  their  elongation  transverse  to  the  magnetic  field  in  the  direc¬ 
tion  of  the  neutral  wind,  and  we  know  why  they  become  unstable.  This  de¬ 
tailed  understanding  provides  a  basis  for  reliable  predictions  of  barium 
ion  cloud  motion  provided  certain  ambient  parameters,  primarily  the 
ambient  electric  field,  E  ,  the  neutral  wind  velocity,  V  ,  and  the  ionospheric 
conductivity  are  known  at  the  location  of  the  planned  release.  Analytic  work, 
numerical  solutions,  and  hydrodynamic  simulation  prove  to  be  useful  for 
understanding  the  kinds  of  motion  and  instabilities  described  by  the  non¬ 
linear  equations  (28)  through  (30).  In  particular,  analytic  work  has  been  suc¬ 
cessful  in  correlating  the  motion  of  ion  clouds,  the  onset  time  of  striations, 
the  dominant  striation  wavelength,  and  the  lifetime  of  rods  with  data 
corresponding  to  a  wide  range  of  ambient  and  cloud  parameters. 

A  combination  of  empirical  scaling  of  the  Secede  data  base  and  our 
theoretical  understanding  also  allows  reliable  predictions  of  the  onset 
time  of  striations  for  barium  ion  clouds  at  midlatitudes  and  for  ioni¬ 
zation  produced  in  other  circumstances.  A  point  that  has  not  been 
brought  out  so  far  is  the  fact  that  now  a  configuration  of  ionization  may  be 
chosen  that  will  striate  quickly  and  yield  high  electron  densities  in  the 
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striations,  and  permit  a  more  reliable  estimate  of  where  the  striations 
will  occur  (see  Chapter  2  of  Vol.  II  of  this  report). 

It  is  important  to  recognize  that  in  large  high-density  ion  clouds  the 
transport  mechanism  perpendicular  to  the  magnetic  field  is  controlled  by 
electron  collisions.  This  fact  explains  why  small-scale  structure  can  be 
observed  for  long  times  and  why  the  ionization  persists  for  hours  in  the 
ionosphere.  We  have  also  shown  that  the  conjugate  ionosphere,  more  than 
30,  000  km  distant,  apparently  participates  in  the  dynamics  of  barium  ion 
clouds . 

2.  Future  Work 

There  are  important  areas  where  additional  theoretical  work  is 
necessary.  The  influence  of  the  initial  ion-deposition  process  on  the 
later  time  behavior  of  the  ion  cloud  is  not  well  known.  Since  the  neutral 
cloud  is  optically  thick  to  the  solar  metastable  exciting  radiation,  the 
barium  ions  are  not  created  uniformly.  They  are  created  initially  on  the 
sunlit  side  of  the  neutral  cloud  and  only  later  are  created  on  the  dark 
side  after  the  metastable  barium  ions  have  been  consumed  by  photoioni¬ 
zation  and  oxidation.  It  is  possible  that  this  process  leads  to  initial 
constant  ion-density  contours  that  are  quite  distorted  from  cylindrical 
symmetry.  Once  this  distortion  is  known,  an  evaluation  of  the  effect  on 
the  later  time  behavior  of  the  ion  cloud  can  be  made. 

The  origin  of  the  long  delay  time  before  striations  occur  in  barium 
ion  clouds  is  not  completely  understood.  The  theoretical  modeling  that 
seems  to  correlate  quite  well  with  the  data  suggests  that  the  delay  results 
from  the  competition  between  the  tendency  of  barium  ion  clouds  to  steepen 
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on  their  backside  and  their  tendency  to  striate.  Both  processes  are  a  re¬ 
flection  of  the  tendency  of  the  higher  density  portions  of  the  ion  cloud  to 
remain  more  coupled  with  the  neutral  wind.  The  model  giving  the  onset 
time  shown  in  Fig.  22  suggests  that  as  long  as  steepening  can  occur,  the 
tendency  to  striate  is  quite  reduced  and  the  ion  cloud  only  begins  to  develop 
structure  after  it  has  steepened  considerably. 

There  are  two  other  possible  explanations  for  the  long  delay  before 
striations  appear.  One  possibility  is  that  there  are  stabilizing  effects  that 
are  operating  at  early  times  due  to  potential  variations  parallel  to  the 
magnetic  field.  A  realistic  model  that  shows  such  an  effect  has  not  yet 
been  developed.  However,  mechanisms  that  can  lead  to  the  stabilization 
of  striations  at  early  times  ought  to  be  investigated.  It  would  perhaps  be 
logical  to  assume  that  the  striations  are  driven  unstable  when  the  backside 
has  steepened  sufficiently  so  that  linearized  growth  rate  increases. 

Another  possibility  is  that  the  long  delay  time  before  striations  occur 

is  just  a  reflection  of  the  extremely  small  initial  amplitudes  of  perturbations 

which  need  to  have  many  e-foldings  before  they  can  be  seen.  Indeed, 

1 2 

Zabusky,  Perkins,  and  Doles  are  lead  to  this  conclusion  as  a  result  of 
the  numerical  computations  shown  in  Fig.  14.  Their  smooth  cloud  did 
not  striate  and  in  order  to  obtain  the  bifurcation  shown  in  Fig.  29,  they 
had  to  introduce  an  initial  perturbation  in  the  electron  concentration.  Two 
arguments  against  this  as  an  explanation  can  be  given.  The  first  is  that  a 
linearized  growth  rate  appears  to  be  so  large  that  the  delay  time  corresponds 
to  tens  of  striation  growth  times.  It  is  unlikely  that  initial  perturbations 
would  need  that  many  e-folding  times  to  become  apparent.  A  second  argu¬ 
ment  derives  from  the  observation  that  once  striations  begin  to  form,  they 
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form  rapidly.  Their  rate  of  formation  is  consistent  with  a  linearized  growth 
rate  but  not  with  the  long  delay  before  they  appear. 

There  is  a  question  as  to  what  does  determine  the  scale  size  of 

striations.  The  linearized  analysis  of  the  two-dimensional  slab-like  model 

presented  earlier  seems  to  give  scale  sizes  in  rough  agreement  with  the 

overall  size  of  the  sheets  that  are  formed.  However,  the  same  linear 

analysis  predicts  that  the  tips  of  the  sheet  should  be  subjected  to  the  same 

instability,  and  should  continue  to  break  up  into  smaller  striations  down  to 

scale  lengths  or.  the  order  of  tens  of  meters.  Thus  the  fact  that  a  300-500 

meter  diameter  sheet  or  rod  appears  to  be  stable  and  does  not  bifurcate 

into  smaller  sheets  is  not  explained  by  the  linearized,  two-dimensional 

5 

analysis.  Perkins,  Zabusky  and  Doles  have  indicated  that  the  large 
conductivity  of  electrons  parallel  to  the  magnetic  field  can  allow  diffusion 
by  ions  in  the  ambient  ionosphere  to  control  the  minimum  scale  size  of 
striations.  The  resulting  estimate  for  the  maximum  stable  scale  size 
appears  to  be  in  agreement  with  the  minimum  unstable  scale  size.  The 
range  of  applicability  of  each  of  these  two  linearized  analyses  ought  to  be 
investigated  further.  An  indication  that  the  striation  scale  size  may  be 
insensitive  to  the  wavelength  of  a  small  initial  perturbation  is  provided  by 
the  numerical  computation  shown  in  Fig.  29. 

The  possibility  exists  that  the  scale  size  of  the  striations  is  deter¬ 
mined  by  the  fact  that  there  is  a  maximum  size  that  does  not  feel  the  influence 

of  the  conjugate  ionosphere.  In  support  of  this  possibility  is  the  observation 
1 9 

by  Chesnut  that  the  motion  of  the  striations  after  they  ha^e  been  formed 
in  Event  Spruce  is  not  colinear  with  the  motion  of  the  neutral  wind  and  the 
main  barium  ion  cloud.  Both  the  theory  describing  the  motion  of  barium 
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ion  clouds  and  the  modeling  of  the  onset  time  for  striations  indicate  that  the 
conjugate  ionosphere  participates  in  the  dynamics  of  barium  ion  clouds. 
Estimates  were  presented  that  show  that  the  typical  voltage  drop  across  a 
barium  ion  cloud  is  of  order  12  volts  wnereas  the  finite  resistivity  of  the 
local  ionosphere  parallel  to  the  magnetic  field  would  cause  a  much  smaller 
voltage  drop  of  order  1.4  volts.  However,  the  voltage  drop  across  a  cloud 
of  ionization  is  directly  proportional  to  its  transverse  scale  size.  Hence 
once  striations  have  been  formed,  they  should  be  decoupled  from  the  conjugate 
ionosphere. 

Each  of  the  above  theoretical  questions  require  further  study.  In 
addition,  an  investigation  into  the  interaction  of  different  ionospheric  layers 
is  needed.  No  work  considering  different  winds  at  different  altitudes  has 
yet  been  attempted.  Analytic  work  is  also  needed  on  the  nonlinear  effects 
exhibited  by  barium  ion  clouds,  for  example,  the  bifurcation  and  the  pinching 
characteristics  of  sheets  need  to  be  more  thoroughly  studied  as  well  as 
the  distortions  of  barium  ion  clouds.  Finally,  there  are  a  number  of  unre¬ 
solved  questions  regarding  the  effect  of  ion  cloud  density  disti  ibutions  on 
HF  radar  scattering.  A  resolution  of  these  questions  would  provide  an 
additional  diagnostic  of  striation  characteristics. 
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I  3  ABSTRACT 


Several  aspects  of  the  analysis  of  barium  ion  clouds  are  presented  including 
ion  cloud  modeling,  comparison  of  radar  and  optical  data,  and  correlation  of 
data  with  theory.  A  quantitative  model  has  been  developed  from  which  various 
properties  of  barium  ion  clouds,  primarily  their  size,  time  history  of  the 
peak  electron  concentration,  ,\nd  height-integrated  conductivity  can  be 
estimated.  These  estimates  are  in  excellent  agreement  with  the  observations 
of  Secede  ion  clouds.  The  modeling  of  the  radiation  transport  in  the  ion  cloud 
has  been  confirmed  by  a  Monte  Carlo  calculation.  The  optical  analysis  of 
Spruce  at  R  +  14  min  has  been  corrected  for  effects  of  cloud  geometry  and  the 
results  for  ion  density  are  found  self-consistent  and  in  agreement  with  radar 
measurements.  A  summary  review  is  given  of  the  values  and  limitations  of 
phot  ^graphic  data.  The  current  status  of  theoretical  understanding  of  the 
dynamics  of  barium  ion  clouds  is  reviewed.  Particular  attention  is  given  to 
their  motion,  deformation,  and  the  properties  of  striations  including  onset 
time,  scale  size,  and  dissipation. 
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